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DYNAMICS OF TILTING PROPROTOR AIRCRAFT IN CRUISE FLIGHT 

Wayne Johnson 


Ames Research Center 
and 

U.S. Army Air Mobility R&D Laboratory 


SUMMARY 


A theoretical model is developed for a proprotor on a cantilever wing, 
operating in high inflow axial flight. This theory is used to investigate the 
dynamic characteristics of tilting proprotor aircraft in cruise flight. The 
model, with a total of nine degrees of freedom, consists of first mode flap 
and lag blade motions of a rotor with three or more blades and the lowest fre- 
quency wing bending and torsion motions \ rotor blade pitch control and aero- 
dynamic gust excitation are included. The equations of motion for a four- 
degree-of- freedom model (lateral and longitudinal tip path plane tilt, pylon 
pitch and yaw) are obtained, primarily to introduce the methods and formula- 
tion to be used in deriving the rotor and cantilever wing equations. The 
basic characteristics of the rotor high inflow aerodynamics and the resulting 
rotor aeroelastic behavior are discussed. The problems of classical whirl 
flutter (a truly rigid propeller on a pylon) and the two-bladed rotor are 
discussed briefly. The influence of the proprotor on the stability deriva- 
tives of the aircraft is considered. The theoretical dynamic behavior of two 
full-scale proprotors is studied, and comparisons are made with the results of 
tests of these rotors in the Ames 40- by 80-Foot Wind Tunnel and with the 
results of other theories. These studies show the sensitivity of the theoreti- 
cal results to several features and parameters of the proprotor configuration 
and to various elements in the theoretical model. In particular, these studies 
demonstrate the important influence of the rotor blade lag degree of freedom on 
the dynamics of both stiff inplane and soft inplane proprotor configurations, 
the dominance of the section lift curve slope (<?£ a ) terms in the high inflow 
aerodynamics of a rotor and the importance of a good structural model of the 
rotor blade and the wing in predicting the dynamic behavior of a proprotor. The 
comparisons also establish the theoretical model developed as an adequate 
representation of the basic proprotor and wing dynamics, which then will be a 
useful tool for further investigations. 


INTRODUCTION 


The tilting proprotor aircraft is a promising concept for short-haul 
V/STOL missions. This aircraft uses low disk loading rotors located on the 
wing tips to provide lift and control in hover and low-speed flight, it uses 
the same rotors to provide propulsive force in high-speed cruise, the lift 
then being supplied by a conventional wing. Such operation requires a 90. 
change in the rotor thrust angle, which is accomplished by mechanically tilt- 
ing the rotor shaft axis. The rotor is vertical for helicopter mode 



operation - landing and takeoff, hover, and low-speed flight - and is tilted 
forward for airplane mode, high-speed cruise flight. Thus the aircraft com- 
bines the efficient VTOL capability of the helicopter with the efficient, 
high-speed cruise capability of a turboprop aircraft. With the flexible’ blades 
of low disk loading rotors, the out-of-plane and inplane (flap and lag) motions 
of the blades are significant, so the blade motion is as important an aspect 
of tilt rotor dynamics as it is for helicopters. When in the cruise mode 
(axial flight at high forward speed), the rotor is operating at high inflow 
ratio (ratio of axial velocity to rotor tip speed); this introduces aerodynamic 
phenomena not encountered with the helicopter rotor, which is characterized by 
low inflow. The combination of flapping rotors operating at a high inflow 
ratio on the tips of flexible wings leads to dynamic and aerodynamic character- 
istics unique to this configuration, and which must be considered in the design 
of the aircraft. The combination of efficient VTOL and high-speed cruise capa- 
bilities is very attractive, so it is important to establish a clear under- 
standing of the behavior of this aircraft and to formulate adequate methods for 
predicting it, to enable a confident design of the aircraft. Experimental and 
theoretical investigations have been conducted over several years to provide 
this capability (refs. 1 to 50). This report develops a model of the aero- 
elastic system for use in some initial studies of the system character and 
behavior. Of particular interest are the features specific to the configura- 
tion: high inflow aerodynamics of a flapping rotor in axial flow and the 

coupled dynamics of the rotor/pylon/wing aeroelastic system. Therefore, this 
work concentrates on the proprotor in airplane configuration: axial flow and 

high inflow ratio. In addition, rigid body degrees of freedom of the aircraft 
are not considered, only the elastic motion of a cantilevered wing. Many fea- 
tures of the coupled wing and rotor motion can be studied with such a model, 
theoretically and experimentally, with the understanding, of course, that the 
model must eventually incorporate the entire aircraft. 

An introduction to the problems characteristic of a high inflow proprotor 
is provided by the following discussion (found in the early proprotor litera- 
ture, e.g., refs. 3 and 8). Consider the behavior of the rotor in response to 
shaft pitch or yaw angular velocity, with the rotor operating in high inflow 
axial flight. A moment on the rotor disk is required to precess it to follow 
the shaft motion. With an articulated rotor (a rotor with a flap hinge at the 
center of rotation), this moment cannot be due to structural restraint between 
the shaft and the blade root, so it must be provided by aerodynamic forces on 
the blade. For example, pitch angular velocity of the shaft will require a 
yaw aerodynamic moment on the disk to precess it to follow the shaft. The 
aerodynamic moment is due to incremental lift changes on the blade sections; 
the component normal to the disk plane provides the yawing moment required. 

For high inflow flight, this incremental blade section lift also has a large 
inplane component and, as a result, the moment to precess the disk is accom- 
panied by a net inplane force on the rotor hub. This force is directed to 
increase the rotor shaft angular velocity, so it is a negative damping force 
that increases with the inflow ratio. There is also the usual rotor positive 
damping due to tip path plane tilt of the thrust vector, plus the damping due 
to the hub moment for a hingeless rotor. If the inflow is high enough, the 
negative inplane force (H force) damping can dominate. The rotor and aircraft 
can be designed so that the velocity for any instability is well above the 
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flight regime, but the high inflow aerodynamics are always important in the 
analysis and design. 

The behavior of the proprotor in high inflow (as outlined above) implies 
the following characteristics: decreased rotor/pylon/wing aeroelastic stabil- 

ity since the negative H force damping of the high inflow aerodynamics can 
reduce the dynamic stability at high forward velocity; decreased damping of 
the aircraft short period modes, again due to the negative II force damping 
contribution of the rotor; and large flapping in maneuvers and gusts. (The 
last arises because the moment to precess the rotor to follow the shaft is due 
to the flapping motion of the blades with respect to the shaft; a given shaft 
velocity requires a fixed component of the section aerodynamic force normal to 
the disk, which means then that increased incremental lift is required at high 
inflow and thus more flapping since flapping is the source of the lift.) 

These features were first delineated in the studies with the XV-3 aircraft 
(refs. 1 to 3), the first experimental tilting proprotor aircraft. Investiga- 
tions of the concept and its problems with the XV-3 provided the initial impe- 
tus for further theoretical and experimental work with the configuration, much 
of which is still in progress. The work with proprotor dynamics has its basis 
in propeller/nacelle whirl flutter investigations (refs. 4 to 7); however, the 
flapping motion of the rotor introduces many new features into the dynamics. 
Experimental and theoretical work has been done by several organizations in 
the helicopter industry on the various features of tilting proprotor aircraft 
dynamics, aerodynamics, and design (refs. 8 to 24). This work has culminated 
in tests of full-scale, flight-worthy proprotors (refs. 25 and 26) and prelim- 
inary design of prototype demonstrator vehicles (refs. 27 and 28) as part of 
the current NASA/ Army-sponsored tilt rotor research aircraft program. How- 
ever, in the literature there is little concerning the details of the analysis 
of proprotor behavior. There are some early reports on very simple analytical 
models (e.g., refs. 8 and 18), and some recent reports on the most sophisti- 
cated analyses available (refs. 29 and 30). Further exploration of the basic 
characteristics of the proprotor dynamics is therefore desirable. 

The objectives of this report are to establish a verified method to 
predict the dynamic behavior of the tilting proprotor aircraft in cruise 
flight; to develop an understanding of the dynamics of the vehicle and of the 
theory required to predict it; and to assess the applicability, validity, and 
accuracy of the model developed. The model of the wing/rotor system developed 
here will be useful for future investigations as well as for these initial 
studies. The primary application of the theory in this report is a comparison 
with tests in the Ames 40- by 80-Foot Wind Tunnel of two full-scale proprotors. 
The analysis begins with a treatment of the four-degree-of- freedom case: 
pylon pitch and yaw plus rotor longitudinal and lateral flapping (i.e., tip 
path plane pitch and yaw). With this derivation as a guide, the equations of 
motion are derived for a rotor with flap and lag degrees of freedom and a six- 
degree-of- freedom shaft motion. The high inflow aerodynamics involved are 
discussed, followed by some elementary considerations of the rotor behavior in 
high inflow. Next, the special cases of classical whirl flutter (no blade 
motion degrees of freedom) and the two-bladed rotor are considered briefly; 
the implications of the basic rotor behavior concerning the aircraft stability 
are investigated. After these preliminary discussions, the development of the 
rotor and cantilever wing model is resumed. The equations of motion for a 
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cantilever wing with the rotor at the tip are obtained and combined with the 
rotor equations of motion to produce a nine-degree-of-freedom model for tilt- 
ing proprotor aircraft wing/rotor dynamics. This model is applied to two 
proprotor designs, in order to examine the basic features of the rotor and 
wing dynamics. Finally, the results of the theory are correlated with those 
from full-scale tests of these two proprotors in the 40- by 80-Foot Wind 
Tunnel . 

The author wishes to thank Troy M. Gaffey of the Bell Helicopter Company 
and H. R. Alexander of the Boeing Vertol Company for their help in collecting 
the descriptions of the full-scale rotors given in table III and figures 14 to 


SECTION 1: BASIC THEORY FOR PROPROTOR DYNAMICS 


Four-Degree-of-Freedom Model 


Consider a flapping rotor on a pylon with pitch and yaw degrees of 
freedom operating in high inflow axial flight. Eventually, at least a few 
more degrees of freedom must be added to this model for both the rotor and the 
support. This limited model is examined first, however, to demonstrate the 
methods used to derive the equations of motion, and because this case is 
studied in the literature. 

The model is shown in figure 1. The pylon has rigid-body pitch and yaw 
motion about a pivot, with the rotor forces acting at the hub forward of the 
pivot. The pylon degrees of freedom are pitch angle ay, positive for upward 
rotation of the hub, and yaw angle a x , positive for left rotation of the hub. 
The rigid-body pitch and yaw motion has inertia, damping, and elastic restraint 
about the pivot. At the hub, a distance h forward of the pylon pivot ( h is 
the mast height) is a rotor with N blades. The rotor has clockwise rotation 
when viewed from the rear, with azimuth angle ^ measured from vertically 
upward. The azimuth position of the mth blade, m = 1, 2, . . . , N is 

= ip + Aip where A\p = 2ir/N is the angle between succeeding blades. The rotor 
degrees of freedom are the out-of-plane motion given by the flapping angles 
6 Cm) for each blade, defined positive for forward displacement of the blade 
tip from the disk plane (upward in helicopter mode, which is the usual heli- 
copter convention) . The blade out-of-plane deflection is assumed to be the 
result of rigid-body rotation of the blade about a point at the center of 
rotation (by the angle gO 7 *)). The dimensionless rotating natural frequency of 
the flap motion is allowed to be greater than 1/rev so that blades with canti- 
lever root constraint may be treated as well as articulated blades (which have 
an actual hinge at or near the center of rotation) . The mode shape for the 
flap motion is assumed proportional to the radial distance r, that is, rigid- 
body rotation. The net forces exerted by the rotor on the hub from all A? 
blades are rotor thrust T, rotor vertical force H, anc rotor side force Y. It 
is assumed in the derivation of the equations of motion that an engine governor 
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supplies the torque required to hold the rotor rotational speed Q constant 
during any perturbed motion, and that the pivot supplies the reaction to the 
rotor thrust T. The pivot also reacts the rotor vertical and side forces so 
that the only pylon motion is pitch and yaw about the pivot. With a flap 
natural frequency greater than 1/rev, as with cantilever root restraint or 
with a flap hinge offset or spring, blade flap motion results in a moment cn 
the hub. The rotor pitch moment on the hub is My and the rotor yaw moment, M x . 

The rotor is assumed to be operating in purely axial flow in the 
equilibrium, unperturbed state, at velocity V . The inflow ratio V/SIR (which 
may be written simply V, with the nondimensionalization implied) is assumed to 
be of order 1. Only rotor aerodynamics are considered; any pylon aerodynamic 
forces are neglected. 

Equilibrium of forces and moments gives the equations of motion: flap 

moment equilibrium for each blade and pylon pitch and yaw moment equilibrium 
(about the pivot). The linearized equations of motion, that is, for small 
angles of the blade and pylon displacement, are then: 

mth blade (m = 1, . . .,/!/): 

Ibi^ + - («y ' 2fta x )cos ip m + (a x + 2f2a^)sin <l> m ] = Mp m (1) 

Yaw: 

+ C-hrr + = AL - hY (2) 

vO %AJ IV tv W IV 

Pitch: 

I y a y + c y a y + K y a y = M y + hH 


where 

Ifo flapping moment of inertia of the blade 

3^ flap motion of mth blade with respect to the hub 

Mp^ aerodynamic flap moment on the blade 

v rotating natural frequency of flap motion (1/rev for an articulated 

^ blade with no hinge spring or offset; greater than 1/rev for a 

cantilever blade) 

I I x pitch and yaw moment of inertia of the pylon about the pivot, includ- 

“ ing the mass of the rotor (as a point mass at the hub) 

Cy 3 C x pitch and yaw damping 

Ky 3 K x pitch and yaw spring restraint of pylon motion about pivot 
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These equations are now made dimensionless with p, fi, and R; the inertias are 
normalized by dividing the flap equation of motion by I and the pylon equa- 
tions of motion by (N/2)I b . The normalization of the pylon inertia, damping, 
and spring constants (division by (N/2)I b ) are denoted by a superscript *; for 
example, Iy A = Iy/(N/2) I b . The rotor Lock number y = paaR h /I b and solidity 
a = Nc/vR are introduced; the Lock number represents the ratio of aerodynamic 
to inertia forces on the rotor blade, and the solidity is the ratio of total 
blade area to disk area. Also notice that the normalized and dimensionless 
hub force H may be written in terms of the rotor coefficient: 


///pfi 2 /? 1 * _ p 0£R^_ it R 2 H _ 

(N/2)I b /pR 5 T b Nc a P TTi? 2 (Q7?) 2 Y oa 


and, similarly, for the other forces and moments. The equations of motion 
then become 


S W * V. 



- C -d y - 

2& tr )C0S + 

(«x ■ 

^X*' d X 

+ ^ x *^ 2 ; 

+ K x*«x = y ( 

2C m x 

oa 

I y* a y 

+ c y**y 

+ K y* a y = y ( 

'2 C M 

-Oil 

oa 


Me T 
x rn 


* m 


2 C 


Y 


oa 


(4) 


These equations are straightforward except perhaps for the pylon acceleration 
terms in the flap moment equilibrium. Blade flap with respect to space is 
composed of gf™), flap with respect to thejuib plane, plus oq and a x , which 
give the tilt of the hub plane; hence the and a x contributions to the flap- 
wise acceleration. The remaining terms are due to Coriolis acceleration; the 
blade has a velocity Or in the hub plane, which has an angular velocity 
a x cos ^ m + a zy s i n Pm due to pylon motion, and the cross-product of these 
gives a flapwise Coriolis acceleration of the blade. In the flap equation, 
the dimensionless aerodynamic flap moment Mp^/pQ^R^ is written as Mp m for 
simplicity; that is, the nondimensional i zation is now implicit in the notation 
^Frrr This practice is followed in the following equations. 


Now introduce a coordinate transform of the Fourier type, defining the 
new degrees of freedom as 


N 


1 N 

nZ 


m= 1 

N 


s = £ v 

ns N Z-i 


.(«) 


f m) 

6 sm nit 


m 


m= 1 


J nc 



6 


(m) 

cos n\\> m 


m= l 



C5) 


6 



so that 


» w = B 0 * S f6 «= C0S ” "m * Si " * V2 t ' 1) " (61 

n 


The coning angle is 6 q> Bic an d B]_£ are tip path plane tilt coordinates; and 
6/7/2 i- s the reactionless flapping mode. The summation over n goes from 1 to 
(N - l )/2 for A 7 odd, and from 1 to (N - 2)/2 for N even; the 6/7/2 degree of 
freedom appears only if N is even. 

The quantities $ nc , 6 ns , and 6/7/2 are degrees of freedom, that is, 
functions of time (which, when dimensionless, is the rotor azimuth angle ifO 
just as the quantities 6 are. These degrees of freedom describe the rotor 
motion as seen in the nonrotating frame, while the terms describe the 

motion in the rotating frame. This coordinate transform must be accompanied 
by a conversion of the equations of motion for ftC™) from the rotating frame to 
the nonrotating frame. This is accomplished by operating on the equations of 
motion with the summation operators: 


1 _ 

N 




• )cosn v 

m 


. )sin nip , 
J r m 



C. • O(-l) 


m 


The usefulness of the Fourier coordinate transformation lies in the 
simplifications it produces in the equations of motion. The above equations 
of motion have periodic coefficients because of the nonrotating degrees of 
freedom in the rotating equations of motion and vice versa; the periodic coef- 
ficients only appear explicitly so far with the pylon inertia terms in the 
flapping equation, but there are actually many more in the aerodynamic forces 
in all the equations. Since the Fourier coordinate transform converts the 
rotor degrees of freedom and equations of motion to the nonrotating frame, the 
result is constant coefficients for the inertia terms, and also for the aero- 
dynamic terms for axial flow through the rotor (as considered here) . In 
addition, only a limited number of the rotor nonrotating degrees of freedom 
couple with the pylon degrees of freedom; in this case, only the 61 c an d 615 
degrees of freedom couple with cty and a x » The other rotor degrees of freedom 
are coupled from the pylon motion and represent only internal rotor motion. 
Thus the transformation reduced a set of N + 2 equations with periodic coef- 
ficients to four equations (considering only those influenced by the pylon 
motion) with constant coefficients. The rotor behavior for this problem is 
basically part of the nonrotating system, so the transformation which converts 
the rotor degrees of freedom and equations of motion to that frame is the 
appropriate one. 
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Operating with (l/N) n (. . .)» ( 2 /N) (. . .)cos i | > m , and 

m m 

(2 /N)^, (• • • )sin on the blade flapping equations gives the nonrotating 

YYl 

equations for coning and tip path plane tilt motion, assuming that N = 3; 


where 


Mt 


h + v b 2 ^o y ac 


Mp 

he * 2 hs * ( V' 1)6 ic - % * 2 “x ' 'I 


(7) 




Ft Q 

hs - 2 V 4 ( V ‘ 1)e is 4 V 4 2 “» = v ~ 


M n • f E % 


m 


Mf 1 C "I E % cos *m 


m 


Mp.„ - — Mp sin \p 

IS N Z-*. t m y m 


m 


The pitch and yaw moments on the rotor disk are and , respectively. 

Note that the transformation introduces Coriolis and centrifugal acceleration 
terms into the 8^ and 8^ equations. The equation for 3 q does not couple 
inertial ly with a y and ax, nor will such coupling be found in the aerodynamics 
hence it may be dropped. A set of four coupled equations remains for the 
degrees of freedom that describe the rotor tip path plane tilt and the pylon 
pitch and yaw motion: $is* a v , anc * If N > 3, the equations for 8 q 

8 yc> and 815 remain as above. To these are added equations of motion for the 
degrees of freedom $ 2C , 8 2 $> • • • > B na ,& ns , and B^ /2 as appropriate; like 
the 8o equation, these equations are not coupled with a y and a x , so they may 
also be dropped from the set, since they represent only internal rotor motion. 
The four-degree-of-freedom model then is sufficient to represent the coupled 
rotor/pylon motion for the general case of a rotor with three or more blades. 
The exception is a two-bladed rotor, N = 2, which is considered separately in 
a later section. 

The equations of motion for the four degrees of freedom (6l C> Big, a y, 
and a x ) are then 
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Mp lc /ao 

Mp is /aa 

2C M y/oa + h{2C H /aa) 
2Cy[Joa - h(2Cy/aa') 


( 8 ) 


The rotor aerodynamic forces (right-hand side) introduce much more coupling of 
the equations. 

The hub pitch and yaw moments due to the rotor. My and M x , might be found 
by integrating the forces on the blade (as is done for the other forces on the 
hub), but it is simpler to express them directly in terms of the rotor flap- 
ping motion. The source of the hub moment is the bending moment at the blade 
root due to flapping, M m = -Zfc(vg 2 -1 )b( w ). Transforming the moment into the 
nonrotating frame and summing over all N blades gives the hub pitch and yaw 
moments : 

= ]Pt-vV - 1)e(m)cos V = - 1 W ' 1)B ic 

m 

► (9) 

M x = ^ [I b (v 6 2 - l)B (m) sin U.J = j I b (v & 2 - 1)B 1S 
m 


where the definition of the tip path plane coordinates $ 1C and & 15 has been 
applied; vg is the rotating natural frequency of the flap motion. If the 
rotor blade has a flap hinge at the center of rotation, then the only spring 
restraint of the blade is due to the centrifugal forces, resulting in vg = 1; 
in that case, no moment on the hub is produced by tip path plane tilt B 1C and 
BlS (except for the torque terms), as required for a hinged blade. With hinge 
offset, hinge spring, or a cantilever root, the natural frequency is greater 
than 1/rev and so tip path plane tilt produces a hub moment. Dividing by 
y(N/2)Ib gives 



( 10 ) 


9 



Rotor aerodynamics- Consider now the rotor aerodynamics. Figure 2 shows 
the aerodynamic environment of the rotor blade section, and the definition of 
the section velocities and forces. A hub plane reference frame is used, that 
is, a coordinate frame fixed with respect to the shaft and tilting with pylon 
pitch and yaw (a y and a^) . All forces and velocities are resolved with 
respect to the hub plane coordinate system, and the blade pitch angle and flap 
angle are measured from the hub plane. The velocities seen by the blade sec- 
tion are up (in the hub plane, positive in the blade drag direction), up 
(normal to the hub plane, positive rearward through the disk), and Up (in the 
hub plane, radially outward along the blade). The resultant of up and u T in 
the blade section is U. The blade pitch angle, 0, is composed of collective 
root pitch, built-in twist, and any increment due to control of the perturbed 
blade motion. The inflow angle is <J> = tan -1 up/up, and the section angle of 
attack, a = 0 - <p. The aerodynamic forces on the blade section are lift L, 
drag D, and radial force Fyi. The latter is positive outward (in the same 
direction as positive up) and has contributions from the tilt of the lift vec- 
tor by blade flapping and from the radial drag due to up. The section lift 
and drag are resolved with respect to the hub plane into normal and inplane 
forces F s and F x . 

The section aerodynamic lift and drag forces are expressed in terms of 
the lift and drag coefficients as 


L = i P e{u 2 + u/)o z = | U* 0i 


D ~ \ 00 {U T Z + u p 2 )a (i = y u 2 c 


( 11 ) 


Working with dimensionless quantities from this point on, the air density p 
has been dropped in the last step in equations (11). The coefficients are 
functions of the section angle of attack and Mach number: 


°d = °d {a ’ m 


where 


a = 0 - tan 


- 1 


Up 


M = M, . U 
tvo 


U 2 = u 2 + u p 2 


and is the tip Mach number, ftfl divided by the speed of sound. The 

section forces resolved into the hub plane are then 
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( 12 ) 


F~ = 


- 


F„ = 


Lu t - Du p ^ 

u 

Lu p + DUp 
U 

uj) 

~TT ~ ) 


The radial force F r has terms due to radial drag and due to the tilt of F z by 
the flap angle 3. The radial drag term in F. p is derived assuming that the 
viscous drag force on the section has the same sweep angle as the local sec- 
tion velocity. Such a model for the radial drag force is only approximate, 
but is adequate for proprotors since this term is not important in high inflow 
aerodynamics. Substituting for L and Z), and dividing by ac, where a is the 
two-dimensional section lift curve slope and c is the section chord, yields 


L z 
ao 


m °A f ± 

U ^ U T 2 a ' U P 2 a 


F X 

ac 


Q n o F 

TTl £ a 

P 2 a + U T 2a, 


V 

L r 
ac 


d z 
Uu v - 6 — 
R 2a ac 


(13) 


The net rotor aerodynamic forces are obtained by integrating the section 
forces over the span of the blade and summing over all N blades. The forces 
required are thrust, rotor vertical force, rotor side force, and flap moment: 


T 



F dr 
z 


H 



Y 






rF dr 
z 


04) 
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or, in coefficient form. 


Yf F -Zar 

oa N / j J ac 

m 0 


2 C 

H 


oa 



aa 




> (15) 


M _F 

aa 



and for the flap equations of motion 


Mp = 4 V Wv 
p ° N / j - m 

m 

M F\ C C0S *m 

m 

% S sin 

m 

The net blade forces required then are, if one substitutes for F z > F x> and F'. 
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(Eqs . (16) continued on next page.) 




The expressions in equations (16) give the net blade force normal to the hub 
plane (thrust) and its moment about the hub (flap moment), the net blade force 
in the hub plane (blade drag force), and the net blade radial force. 

To evaluate the blade forces, the blade section pitch angle and the 
velocities seen by the blade section are required. Each velocity component 
has a trim component and a perturbation component, the latter due to the blade 
and pylon degrees of freedom. When the differential equations of motion are 
linearized, the perturbation components of the velocity are assumed to be 
small. The trim velocity components for operation in purely axial flow are 

Wy = 

Up = V + v 

U R = 0 

The velocity wy is due to the rotation of the blade; the rotor rotation speed 
is included here to show the source of this velocity, but it is usually 
dropped when dimensionless quantities are used. The inflow u P is composed of 
the forward velocity V plus the induced inflow y; the latter given by momentum 
theory as 


or 


V = -7/2 + / (7/2) 2 + C T / 2 


(17) 


V + y = 7/2 + / (7/2) 2 + C T f 2 
= 7 + C t /2V 

where the last approximation is valid for large inflow V (really, the inflow 
ratio V/tiR, since it is dimensionless). The induced inflow will, in fact, be 
very small, v/V « 1, for typical proprotor operation; this is due to the high 
inflow 7, and also to the low working Cy of a proprotor in cruise. 
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Consequently, induced inflow is not generally an important factor in high 
inflow proprotor aerodynamics, and the assumption of uniform induced inflow, 
or even neglecting it entirely, is reasonable for an investigation of the rotor 
aeroelastic behavior. Since the rotor in the unperturbed state is operating 
in purely axial flow, the radial velocity component Up has no trim term. The 
trim blade pitch angle is determined by the collective pitch and the blade 
built-in twist. 

The perturbation velocities are due to the rotor and pylon degrees of 
freedom (g 3 a y 3 a x here) and to the aerodynamic gusts. The convention used 
for the gust velocities is shown in figure 1. The gust velocities are normal- 
ized based on the forward speed 7, so that the vertical and lateral gusts (ot£ 
and are angles, and the longitudinal gust {uq) is a fractional change in 
the forward speed. This convention follows the usual practice for aircraft 
stability and control investigations. The gust velocities are a small pertur- 
bation to the direction and magnitude of the forward velocity 7, assumed uni- 
form over the entire flow field. The gust influence is entirely aerodynamic; 
the gust velocities do not involve a change of the aircraft velocity with 
respect to an inertial frame, but only a change with respect to the air. 
Therefore, the gust velocities do not appear in the inertia terms of the equa- 
tions of motion, but only in the aerodynamic terms. The perturbation 
velocities are 
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( 18 ) 


angular velocity of the pylon about the pivot; inplane component of the for- 
ward velocity V + V due to the tilt of the pylon; and the inplane velocity due 
to vertical and lateral gusts. In Sup there are two terms: flapwise velocity 

with respect to the air, due to both flapping with respect to the shaft and 
angular velocity of the shaft (this term is proportional to r) ; and longitudi- 
nal gusts (this term is independent of r) . If 8up in equations (18) is 
written as 


8u p 


,6u pg * Su pA 


(18a) 


14 



then Su T , Su Rs Sup B , and &up A are all independent of r, and so may be factored 
out of the integrands m the aerodynamic forces. The perturbation of the 
blade pitch motion is 

60 = 6 - Kp$ (19) 

where now 0 is just the perturbation of the blade pitch, an input variable in 
the equations of motion, also available for feedback control. Since this 
pitch perturbation is made through the control system, it is uniform over the 
blade span (independent of r ) . Also included above is pitch/flap coupling 
with K P the gain of negative feedback of blade flap angle to pitch angle. 

This feedback is usually accomplished by mechanical means inherent in the 
control-system geometry; it is then usually referred to as 6 3 coupling where 
here Kp = tan 6 3 . 

It is now possible to find the perturbations of the aerodynamic forces on 
the blade, that is, the forces due to the rotor and pylon degrees of freedom 
and gusts. The following relations are made: 


Sa SL = 9^r 6oi + a/ m 


3c , 3c , 

6c , = ~ 6a + — ~ m 

a 3a 3 M 


6a = 66 - 


u T &u p - u p Su T 
U 2 


SU = 


U T^ U T + u p^ u p 
~U 


SM = M , . &U 


( 20 ) 


where the coefficients of the perturbation quantities are evaluated at the 
trim state. Hence the net blade forces may be expended as linear combinations 
of the perturbations of the rotor blade velocity and pitch angle: 



\ 


l (21) 
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The coefficients are constants, independent of rotor azimuth since trim 
axial flight is considered; they are integrals of the blade aerodynamics over 
the span (expressions for them are obtained later). The first terms, subscript 
o y are the trim forces and moments on the blade* The second terms, subscript 
y , are forces and moments due to hub inplane velocity, the third terms, sub- 
script 3, are due to flapwise velocity of the blade; the fourth terms, sub- 
script X, are due to axial velocity of the rotor; and the last terms, subscript 
0, are due to blade pitch control. The thrust forces on the blade are T, the 

the flap moments, M ; the blade drag forces, H ; and the blade radial force R . 

In the blade radial force, the trim value of the coefficient of 3 
(Z 1 F~/ac dr) is required, which is just C^/aa\ therefore, 

c lF , 

/ — dr - R - — 3 (22) 

ac v R oa 

The last term is the radial tilt of the blade thrust vector. 

The blade forces may now be summed over all N blades to find the net 

rotor forces. If the expressions for the blade forces (eqs. (21)) are substi- 
tuted into those for the rotor forces (eqs. (15)), the aerodynamic coefficients 
are independent of m (blade index) so the summation operates only on the per- 
turbations of the blade velocities and pitch. If the definitions of the rotor 
nonrotating degrees of freedom (the Fourier coordinate transform described 
above) are used, the following is obtained: 

ft., \ 

aa o I 


= {H + R ) [-ha + [V + v)a + Va r ] 


* V ®1 s ■ B 1 c * ix) 

* V 6 1 s - w 


- 6 1C 


+ + iv + v)a x + 7 b g ] 

■ V*ic + B is ' V 

' *V e i c " *t b i Z 

■ - B .S 


(Eqs. (23) continued on next page.) 
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— — = M [ -hct + (7 + u)a + 7a ] 
aa m l y y G 
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+ V e is - W 


( 23 ) 


The perturbations of the thrust are no longer needed (they are entirely reacted 
by the pivot) so only the trim term, T Qi is retained. Again, N - 3 has been 
assumed in evaluating the sums. The fixed frame coordinates for the blade 
-pitch motion are 


0 O = J 


r £ 
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e = § W"’ cos * 
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(24) 


These coordinates represent control inputs by means of the usual rotor swash- 
pl&te mechanism: 0 q is the rotor collective control, and 0 ^ and 8 are 

rotor lateral and longitudinal cyclic control (control plane tilt) . 


The decoupling of the equations of motion (which has been seen in the 
inertia terms) is maintained by the aerodyanmics also because axial flow is 
assumed for the trim operating state. In the forces that excite the four- 
degree-of-freedom model (£#, Cy > A and Afp^) the only rotor degrees of 

freedom involved are and 615 . The aerodynamic forces introduce some input 
variables, but even with these, there is limited coupling: only lateral/ 

longitudinal Control plane tilt ( 0 ^ 7 , ©i 5 ) and lateral/longitudinal gusts 
appear in the four-degree-of-freedom set. As for the inertia terms, the aero- 
dynamic terms due to the higher rotor degrees of freedom (02 c> S> • • -j $no> 

$ns> $N/2 as a PP ro P r i ate f° r > 3) do not involve any coupling with the shaft 
motion or with the blade pitch control or gusts (assuming conventional 
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swashplate control inputs and uniform gusts); hence these degrees of freedom 
remain internal rotor dynamics. 


From helicopter rotor aerodynamics, the tilt of the tip path plane ( 6 ^ or 
$1$) is expected to tilt the rotor thrust vector and hence give an inplane 
force on the rotor hub. The tip path plane tilt terms in C # and Cy are (from 
eqs. (23)): 



aa 
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1C 


A 
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The first terms are the inplane forces due to radial tilt of the blade mean 
thrust vector by the blade flapping. They are only half that expected because 
of the tilt of the rotor thrust by tip path plane tilt, assuming that the 
thrust vector remains perpendicular to the tip path plane. The other half is 
in Hg. Rotor tip path plane tilt 6 ^ or 6 ^ 5 , steady in the fixed system, 
causes a flapping velocity in the rotating frame. This flapping velocity 
changes the blade angle of attack and so tilts the blade mean thrust vector in 
the chordwise direction (like induced drag). The inplane force due to flapping 
velocity, Hg, may be written 


C T 

//•= — + H‘* 
B aa B 


(25) 


where the first term is the tilt of the blade thrust, and Hg* is due to the 
rotor inflow. Thus the inplane hub forces due to tip path plane tilt are, 
combining that due to direct radial tilt of the blade thrust by 6 , and that 
due to chordwise tilt of the blade thrust by 0: 



The first term is the inplane component of the rotor thrust due to tip path 
plane tilt, as expected, and the second is the inplane force due to the inflow 
term of Hg acting on 6 . The inflow term Hg* is negative, so it decreases the 
inplane force due to tip path plane tilt. For low inflow, the effect of Hg* 
is small, but for large inflow (as considered here) it dominates the thrust 
vector tilt term. It is, in fact, the negative // force, already mentioned as 
an important feature in high inflow rotor aerodynamics . Notice that Hg acts 
on the blade flapwise velocity to produce an inplane force, regardless of the 
source; hence the angular velocity of the tip path plane (with respect to the 
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hub plane) or the shaft (the hub plane) also produces a hub force, with no 
corresponding term from the blade radial force. 

Substituting now for the rotor forces and moments into equations (8), one 
obtains the equations of motion for the four-degree-of- freedom model: 



The influence of the rotor aerodynamics in this set of equations is as follows: 
damping of the flap motion, Afg, which also acts on flapping velocity due to 
shaft angular velocity and tip path plane tilt; speed stability flap moments 
M due to hub velocity produced by angular velocity of the pylon about a pivot 
aft of the hub, and due to the inplane component of the forward velocity pro- 
duced by the shaft tilt; positive damping and a negative spring on the pylon 
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motion due to + 7? y (which is positive); and inplane hub forces on the pylon 
acting through H g, due to flapping velocity produced by tip path plane tilt 
(Bl c> 615 ) or angular velocity 3 1 / 7 ) or shaft angular velocity (d^, d^.) . 

The blade pitch input produces flap moments and hub forces through and 7 /q ; 
through pitch/flap coupling Kp these coefficients enter the coefficients of 
and S>ip also. Gusts produce flap moments and hub forces through the hub 
inplanc velocity coefficients, A/y and //y + 7?y. 

Aerodynam ic coeff lei ents - The blade forces required are given in 
equations (16). Substitute for the velocities and section force coefficients 
in terms of the trim plus perturbation values, linearize about the trim state, 
and compare with the expanded forms (eqs. ( 21 )) to identify the aerodynamic 
coefficients. For the moment, the effects of drag ( yv) and of compressibility 
(14) will be neglected. Moreover, only certain coefficients are required for 
the four-degree-of-freedom model, namely, T QJ Q 0i _, // 3 3 //y_, and // 0 - 

Only these coefficients will be examined now; in fact, this set is sufficient 
to describe the general behavior of all the rotor aerodynamic coefficients 
required here. 

If the drag coefficient is neglected (except for Q 0 ) , the forces 
required are 



where ci = c^(a) 9 a = 0 - tan 1 up/uy, and U 2 = wj> 2 + up 2 . With the perturba- 
tions about the trim state, the flap moment becomes 
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= ^0 * + M^Sup^ * 60 

and, similarly, the blade drag force is 
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Also required is R^ y but R, is due to blade radial drag force and since all 
drag terms are neglected for now, R^ = 0. 



The trim values of the velocities (uy, up y and U) and of the blade loading 

and ) should be substituted into these expressions. For the velocities, 
the trim values are 


Uf = v 
up = V + v 

V = Jr 2 + (V + v ) 2 
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For the lift coefficient and lift curve slope, assuming small angle of attack 
so that C£ a = a (two-dimensional lift curve slope), the trim values are 


2a 2 

Q 

= | = \ [9 - tan - 1 {V + v)/r] 

where 0 is here the trim collective value plus the built-in twist of the blade. 

The inplane force due to flapping velocity is also written 
H g = {Cy/aa) + // g*; substituting the expression for Cy/aa from II g (with uy /') 
yields 



Approximate expressions for the aerodynamic coefficients may be obtained 
by evaluating the integrands at an effective radius. Since the inflow angle 
is <p = tan' 1 up/uy , then also Uy/U = cos <j> and up/V - sin from this, it is 
possible to substitute for up and U in the integrands, in terms of uy and <£. 
Then a £ and <p are evaluated at an effective radius v Q (r e = 0.75 usually) and 
uy = y* is used in the integrand; Co /2a = 1/2 may also be used. Then, 

/ I 

r 2 a , s , 1 a 

cos <p 2 ' ll cos <p 6 
u 

where a is the rotor mean angle of attack. This expression is used to evaluate 
a in the aerodynamic coefficients. The coefficient TT is approximately 



6 aa 


The last step follows since the C T term (a term) is significant only for low 
inflow, when <f> is small; therefore, [1 + (cos <j>) -2 ] = 2. Similar approximations 
may he found for the other coefficients. 
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The aerodynamic coefficients are then approximately: 
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a *r . v S1 l n ♦ . Zz v 
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0 4 cos 


where 4 is evaluated at r py so that 


cos <p = p c /Jv c 2 + V 2 
sin 4) = 7/ 2 + V 2 - 


and usually ry., = 3/4 is satisfactory. The thrust coefficient (a^) terms in 
Mg (flap damping) are always negligible and were therefore dropped. There are 
no terms in the pitch coefficients Mq and //g . For operation in high inflow, 
V is of order 1 and so cos (p and sin <f> arc of order 1 also; C^/oa << 1 (in 
fact, for proprotor operation in cruise, the blade loading is even lower than 
usual for helicopters). Hence for high inflow operation, the thrust effects 
on all the aerodynamic coefficients terms) may reasonably be neglected 
compared with the <r?£ a terms. The reason the <?£ terms dominate the coeffi- 
cients is that with high inflow both inplane and out-of-plane velocity pertur- 
bations give large angle-of-attack changes. Therefore they give (through <?£ a ) 
large section lift perturbations, which have significant components in both 
the out-of-plane and inplane directions. The high inflow thus allows a great 
simplification of the rotor aerodynamic derivatives. 
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Retaining now only the terms in equations (27), that is, assuming 
that V is of order 1 and a « 1, and substituting for the trim velocities, one 
obtains the coefficients (writing V for V + V for convenience): 



If <?£ is assumed to be independent of r, that is, o^j2a - 1/2, the integrals 
may be evaluated exactly as 


V 1 ? V 3 1 * ^ * f2 
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The behavior of these coefficients is clearer from the expansions for small 
and large V : 
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This behavior is also shown in the approximate expressions based on the inflow 
angle <p at an effective radius. These expressions (eqs. (29)) are in fact most 
convenient for examining the general behavior, since they contain (in simple 
form) the behavior over the complete range of inflow ratio. For numerical 
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work, it is straightforward to evaluate the coefficients by use of the exact 
integrals (eqs. (31)) or even more complete expressions (as derived in a later 
chapter) . 

The flap damping coefficient Mq is negative (which is positive damping); 
the inplane force due to flapping velocity Hq has the term C T /aa as expected, 
and the inflow term H g* is negative and therefore opposes the contributions 
from tip path plane tilt of the rotor thrust. The speed stability coefficients, 
M v and H y, and pitch control power coefficients, Mq and Hq, are all positive. 

All coefficients are of order 1 for high inflow. For low inflow only, the 
flap damping and control, Mq and Mq, are of order 1; the flap moment due to 
inplane velocity is an order V smaller in low inflow, and all inplane force 
coefficients are an order V smaller than the corresponding flap moment coeffi- 
cients. Flap damping Mq and the mean blade angle of attack (for given rotor 
thrust) are decreased by high inflow, but remain the same order as for low 
inflow; the other coefficients increase with increased inflow ratio. For low 
inflow, the rotor thrust coefficient terms must be retained for My, ffy, and Hq 
(but not for Hq*; Hq, of course, has only C£ a terms), but for hign inflow, they 
may be neglected for all coefficients. 


SECTION 2: THEORETICAL MODEL FOR A ROTOR IN HIGH INFLOW 


Equations of Motion and Forces for the Rotor 

With the procedures to be followed established from the derivation and 
discussion of the simpler four-degree-of-freedom model, consider now a more 
comprehensive model for the rotor motion. The blade motion is extended to both 
flap and lag (first mode out of plane and inplane) degrees of freedom, and the 
shaft motion to all six degrees of freedom; inputs from blade pitch and aero- 
dynamic gusts complete the model. The equations of motion are derived for the 
rotor degrees of freedom, and expressions for the rotor forces and moments 
acting on the hub are obtained. In a later chapter, a wing is added to this 
model, thereby completing the equations for use in the study of proprotor 
dynamics. 

The model considered and the conventions for the hub forces and moments, 
pylon motion, and aerodynamic gust are shown in figure 3. The pylon motion is 
defined about a pivot a distance h aft of the hub. The pivot linear displace- 
ment degrees of freedom are xp, yp, and zp - vertical, lateral, and longitudi- 
nal, respectively. The angular degrees of freedom are a^., a y, and a z (yaw, 
pitch, and roll) . The forces and moments exerted by the rotor on the hub and 
the gust velocities are as defined in section 1. The torque reaction between 
the rotor and pylon is Q; following shaft-driven rotor convention, Q is the 
torque exerted by the shaft on the rotor, hence the torque moment on the hub 
due to the rotor is -Q (as indicated in fig. 3) . The rotor blade azimuth 
angle ip is measured with respect to the pylon, which is rotated by a z in roll; 
so the rotational velocity of the blade with respect to space is SI + a z (with- 
out blade flap or lag motion). The equilibrium velocity V with respect to the 
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air is assumed to be purely axial flow. The blade motion is defined by flap 

and lag (degrees of freedom) and pitch (input) motion with respect to the hub 
plane. 

The blade motion is represented by two degrees of freedom per blade: 
flap and lag motion 6 and which are pure out-of-plane and pure inplane 
deflection of the blade spar, respectively. The motion is defined with respect 
to the hub plane. The mode shape of the blade deflection is ng(r’) for flap 
and n^(f) for lag. These modes are functions of r and are normalized to 1 
at the tip. The out-of-plane deflection of the blade is then a distance 
BOIOngC**) normal to the hub plane, with 6 defined positive for deflection above 
the hub plane (forward in airplane cruise mode) . The inplane deflection is a 
distance ?(^)n^(f) from the undeflected spar line, measured in the hub plane, 
with £ defined positive for deflection opposing the rotor direction of rota- 
tion Rotating mode shapes are used, that is, natural vibration modes includ- 
ing the centrifugal spring due to blade rotation. A major influence on the 
mode shape is the root restraint, that is, either a hinged or a cantilever 
root. However, the centrifugal stiffening is so strong that the effect of the 
root restraint on the lowest flap and lag mode shapes is restricted mainly to 
the root area. The influence of the root restraint on the natural frequencies 
of the modes is of primary importance. The first (lowest frequency) flap and 
lag bending modes even for a cantilever blade are then nearly r\ = r; near the 
root of a cantilever blade, the mode shape must deviate from this, of course 
to satisfy the boundary condition of zero slope. 

form for the equations of motion is in terms of the nonrotating 
rotor degrees of freedom. It is possible to have different mode shapes for 
the various nonrotating degrees of freedom, for example, one for the coning 
mo e and one for the tip path plane tilt modes, depending on how the hub 
restraint appears during deflection of the blades in that particular rotor 
model. Two rotors are considered in applications of this theory; a canti- 
lever rotor and a gimballed rotor. For the cantilever rotor, the mode shape 
tor all nonrotating degrees of freedom of the blade is that of elastic bending 
with cantilever root restraint. For the gimballed rotor, the mode shape for 
tip path plane tilt degrees of freedom £>i q and is that of an articulated 
blade, namely, rigid-body motion about a hinge at the center of rotation, 
ng = v For all other nonrotating modes of the gimballed rotor (specifically 
for the coning and blade lag modes), the rotor blade acts as a cantilever 
blade, with corresponding blade deflection mode shapes. 


The motion of a cantilever rotor blade in elastic bending is actually 
more complex than the representation used here. The inplane and out-of-plane 
deflections are highly coupled by the collective pitch and built-in twist of 
the blade, which are large for the proprotor. Consequently, although the low- 
est bending modes are usually still identifiable as predominantly flap or lag 
motion, there is actually both inplane and out-of-plane motion in each mode. 

The neglect of this effect, by assuming that the blade flap and lag degrees of 
freedom are pure out-of-plane deflections and pure inplane deflections, 
respectively, is probably the severest limitation of the theory presented here. 
The basic features of the flap and lag motion are represented, so this model 
may be expected to predict proprotor behavior fairly well. 
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The equations of motion are derived for a constant rotor rotational speed 
n (with respect to the pylon); this is to be considered the model for powered 
operation of the rotor. The autorotation case - where the rotor rotates freely 
on the shaft, the rotor speed being determined by equilibrium of torques on 
the rotor - can also be handled with this model if the collective lag mode Cq 
is used. This mode involves the simultaneous motion of all blades in the lag 
direction (opposite the rotor rotation direction); if there is no hub restraint 
for this mode, it will be equivalent to a perturbation of the rotor azimuth or 
rotational speed. If the mode shape for rigid-body rotation is used, and if 
the natural frequency in the rotating frame is set to zero (n = P and v = 0 for 
the z„ mode), then, indeed, t 0 will be just the degree of freedom that repre- 
sents the rotor rotational speed perturbation. This is a good representation 
of the autorotation case. The other limit, a fixed rotor rotation speed 
will be considered as powered operation. With a constant rotor rotation speed 
the collective lag mode is then elastic bending of a cantilever blade with 
respect to the hub (which rotates at a constant speed). This limit is, in 
fact, the case of operation with a perfect governor on the engine or rotor 
speed. For an actual rotor in powered flight, the engine/drive train/governor 
dynamics must be included to give a complete representation of the behavior. 

The blade also ha? pitch motion about the feathering axis at the blade 
root (given by 0 ), with the actual blade pitch measured from the hub plane. 

The pitch has trim and perturbation contributions as before. The trim value 
is due to root collective and built-in twist; the perturbation value is due 
to a control input and pitch/flap coupling. Pitch/flap coupling ( 63 ) is 
included for the gimballed rotor. 

The equations of motion for flap and lag degrees of freedom are obtained 
from equilibrium of moments on the blade. For the w?th blade (m - 1 , . . • > ^0 
in the rotating frame, the equations are 


J B (B + v g 2 0) + J Bct [-ca y - 2<ycos + (d x + 2<ysin Tig + s^z p = m f 
I d + v^ 2 ?) + S ? [(* p + h'dy)sin i> m - (y p - ha^ cos = M L f 

The flap and lag aerodynamic moments on the blade are 


( 32 ) 
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and the inertia constants are integrals of the blade section nuyrs 
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The rotating natural frequencies of the flap and lag motions are vg and v r 
respectively. A subscript o will be added to the mode shape inertias and 
natural frequency for the collective modes (coning or collective lag)' sind 
these terms may not be identical to those for the cyclic modes (e g for the 

gimballed rotor or the autorotation case). 1 8 ’’ th 


and the lfl h qUatl ° n as fo ^ ced b y Pure out-of-plane aerodynamic forces (F,) 

flan and in h ^ lnp ^ me forces ™ ’ because of the assumption of decoup ed 
flap and lag bending modes. The flap equation is as before, 1-ith the addit on 
of the acceleration due to longitudinal motion of the shaft. The flap mode 
shape ng influences the effective inertias of the flap motion and the^haft 
angular acceleration; with rigid-body flap motion, n « = r, the equation reduces 
o that used for four-degree-of-freedom case. The lag motion couples with 

witl^rol ac ^ ele y a 1:1011 op the rotor hub (resolved into the rotating frame) and 

I 1 . ang ^ ar acce leration of the rotor shaft. The Coriolis inertial cou- 

P . " g t the fl ?5 w nd 138 equations has been neglected. The coefficients of 
these terms would be proportional to the rotor trim coning angle which is of 

a^for M T ' “° WeV6 ^ aerodynamic terms also contribute to fhis coupling, 

- . • lg in f^ ow these coefficients are of order 1. Hence the Coriolis 

forces! C ° UP ing may be ne S lec ted compared with the high inflow aerodynamic 

r 1 

Now let I b = jf r 2 m dr and normalize the inertias by dividing by I b ; this 
normalization is denoted by superscript *, for example, = I$/I b . The Lock 
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number is defined (as before) by y = p acR^/Ifo. Ifo is only a normalization 
factor - a representative moment of inertia of the blade. It is used in the 
blade Lock number (the ratio of blade aerodynamic to inertial forces) and to 
normalize the blade masses so they are of order 1. A convenient inertia is 
that of the blade about the shaft, that is, the rotary moment of inertia of 
the entire rotor divided by N. This inertia is a well-defined property of the 
rotor and also should be the largest possible moment of inertia of the blade. 


This normalization yields the following equations: 


V (B * v e 




« + - 2 «P COS + + 2 V sin 'V + 5 S *r = Y w 


M r 


I*(i + v 2 0 + S *[( x p + fto y )sin ^ - (y p - ha J cos - 1^ = y — J 


> (33) 


If n s v for the flap and lag modes, then the I* terms are all nearly 1 and the 
C A terms, nearly 3/2 (for constant mass distribution); with usual blade con- 
struction, the I* terms are slightly less than 1, and the S* terms around 1. 


The Fourier coordinate transformation is now applied to convert the 
equations of motion and degrees of freedom of the rotor from the rotating 
frame to the nonrotating frame. Again, the nonrotating degrees of freedom 
above 0, 1C, and 1C are not coupled with the shaft motion, so these higher 
rotor degrees of freedom involve just internal rotor motion. In studies of 
the coupled rotor and shaft motion then, the 0, 1C, 1 C set is sufficient to 
treat the general case of N > 3. The nonrotating equations of motion for the 
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where the flap aerodynamic forcing 


moments are 


M p = 4 V Mp 

0 M ‘ h m 

m 

~ ^ T COS 

10 iv < 1 m 

m 


*i8 - I E % sin *™ 


m 

and similarly, for the lag moments. The flap modes are coning (g n ) and tin 
path plane tilt (gy-., Siy). The collective lag mode ? 0 is simultaneous lagging 
mo ion of all the blades (with respect to the hub. rotating at constant speed a 
or the powered case; for the autorotation case, i Q is the rotor speed pertur-’ 
ation degree of freedom). The cyclic lag modes £ lc - and Cio produce rectilin- 
eal inplane motion of the net rotor center of gravitv, laterally for r ir (- uu 
direction) and vertically for r, • (x h direction). Note that the equations 
-separate into a lateral/vertical group (1C, IS, x, y) and a longitudinal group 
1U, i), with no inertia coupling between them. This decoupling is maintained 

b> the aerodynamics also (because of the trim axial flow); the shaft motion due 
to the actual wing degrees of freedom will, in general, couple the two groups 
of equations. 6 1 


The hub moment due to the rotor may be expressed (as before) in terms of 


tip path plane tilt 


Bjr and fci- 
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The inertia contributions to the rotor drag, side force, thrust and 
torque acting on the hub are ’ 
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( 36 ) 



The drag and side forces are the net inplane acceleration of the rotor due to 
the motion of the shaft and blade; similarly, the thrust is the net longitudi- 
nal acceleration of the rotor; and the torque, the net angular acceleration. 
The new inertia constants are 


^ = 


Mu = 


f, 

i: 


r 2 m dr 


m dr 


Therefore, NI 0 is the moment of inertia of the entire rotor about the shaft ^ 
and NMu, the mass of the entire rotor. When normalized (divided by I b ) , T 0 
is nearly 1 (exactly 1 if I 0 is used for I b ) and M b * is around 3 for a uniform 
mass distribution is greater than 3 for usual rotors]. In coefficient 

form, dividing the side and drag forces by {N/2)I b y and the thrust and torque 
by NI b y, these forces are 
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Rotor aerodynamics- The analysis follows that of the previous section, 
the section aerodynamic environment, with the conventions for forces and 
velocities, remains as shown in figure 2. With the present degrees of free- 
dom and shaft motion, the perturbation velocities are 


Su ? = r(a 2 - O - h(i y sin ^ cos * n ) 


(V * y)(a sin * a x cos 4-^) * V(.& G cos sin * m ) 
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Su R = h(-a cos < 1 ^ sin 4^) + O'’ + lO (“^ cos % " ° x Sln ^ 
+ 7(-B c sin ip n + a Q cos - Cp p sin i | > m * i p cos 4^) 
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coupHng^:' PitCh pertUrbati °" is as before . including control and pitch/flap 


66 = e - K p Q 

The trim velocities are, again, for equilibrium axial flow 
and Up = 0. ' 


Uf = r, up = V +v y 


The rotor motion contributions to the velocity perturbations (in 6ur and 
upp) assume that ng - n ? = r for the flap and lag mode shapes This aDnroxi - 

moments on the blade so that 6 sed in the aerodynamic 
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The use of this mode shape for the aerodynamic greatlv simplifies the 
aerodynamtc coefficients involved or, at least, reduces th'c number of coeffi- 
lents required With the correct n B and nr, separate coefficients are 

equired for blade motion and shaft angular motion, and for the lag moments 
and torque moments on the blade. With n. = n - « n „, , , moments 

appears in the integrands of the aerodynamic coefficients^ never^r*^ 
hence the evaluation of the coefficients is also simplif^d. " ’ 

The expressions for the section aerodynamic forces (L D and F th P ir 
^ecomposition^into the hub plane F and F x , and the rotor forces and^moments 

’ ’ - ’ u ^ in teims the net rotating forces on the blade arc the 

same as m the previous section. Again, the net blade forces may be funded 
as linear combinations of the velocity and pitch perturbations: ' 1 


fd 

l' 

fd 

r. 


Jr=M 0 * M tu T 


H 

P 1 B 


M k Su ?A 


H k SUT A 


M x 6u Pa 


+ 'V 9 


h \ &u Pa + !I b 6u Pb * /7 e 6 0 


do dV = T ° + \^B + + T A &U PA + T B S “P H + ' e 50 


v — ~ dp = 


+ Q v Su Tg + Q^u Ta * Q x 6up A +■ QgSup B + y 0 66 


~dr- = R SUr , - 6 

ii ft oa. 


34 



where M is the flap moment; H, the blade inplane force (drag direction); T, the 
blade thrust; Q, the blade torque moment; and R, the radial force. The sub- 
scripts denote the source of the force or moment; subscript 6> indicates trim 
values; subscript y, hub inplane velocity (speed); t> blade rotational velocity 
(lag damping); 3* flapwise velocity (flap damping); X, hub longitudinal velocity 
(inflow) ; and 6, blade pitch control. The coefficients may be grouped as 
inplane and out-of-plane forces, so the // and Q terms have similar behavior, 
and the M and T terms have similar behavior. Alternatively the coefficients 
may be grouped as inplanc and out-of-plane velocities so the coefficients with 
subscripts y and i have similar behavior, and those with subscripts 3 and X 
have similar behavior. The only difference between the coefficients within a 
particular group (say, the out-of-plane forces due to out-of-plane velocities: 
Mft, M\ y Jg, and is a factor of r more or less in the spanwise integration 
(the difference between the force and moment, and between the translation and 
rotational velocities), hence just slightly different numerical constants. Ihe 
behavior of the coefficients with a variation in the parameters (in particulai , 
with forward velocity V ) is basically the same within a group; that is, it is 
determined primarily by whether an inplane or out-of-plane force is involved, 
and whether an inplane or out-of-plane velocity or blade pitch control is the 
input. The fundamental set of coefficients is considered to be the M and // 
terms with subscripts y, 3? and 0 - one each of inplane and out-of-plane types - 
together with T 0 and Q 0 for the trim values. Then the behavior of all other 
coefficients may be inferred from a knowledge of the behavior of this set. 


Again, the inplane force due to flapping velocity is written: 
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to show explicitly the contribution due to the thrust vector tilt. 
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The blade forces can now be summed over all N blades to find the net rotor 
forces. The aerodynamic coefficients are independent of i p m , so the summation 
operates only on the blade velocity perturbations. With the definitions of the 
rotor nonrotating degrees of freedom, the flap and lag moments required for the 
equations of motion are 
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and, similarly, for the lag moments, with the M terms replaced by Q terms 
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The aerodynamics maintain the separation into lateral/vertical and 
groups because of the purely axial flow in the trim state. 


longitudinal 


Equations of motion - The complete equations 
by combining the aerodynamic and inertia terms, 
motion group, the equations of motion are 


of motion can now be obtained 
For the lateral/vertical 
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The longitudinal equations of motion are 
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Note that if, in Pg/aa, the substitution is made for Co from the equation of 
motion : 
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where the last expression holds if n^ = v is assumed for the inertias, as it 

was for the aerodynamic terms. That assumption is consistent with the use of 
only one blade lag mode; this result is comparable to that for the hub moments 
M x and M ' in terms of rotor tip path plane tilt and 3 1P- 


The present model may be used for the case of the rotor operating in 
autorotation, as is frequently done for proprotor dynamics wind-tunnel tests, 
by use of the collective lag mode for the rotor speed perturbation degree of 
freedom. It is important to include that degree of freedom so that proper 
representation of the autorotation dynamics is obtained. In fact, the 
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representation used for the powered case, assuming that the rotor hub 
rotational speed ft is fixed by the engine and that the collective lag motion 
is elastic rotor deflection with respect to the hub, can only be considered a 
limiting case of a perfect governor. A model with true engine/ governor/drive 
train dynamics would always allow some rotor speed perturbation, and hence that 
motion would also have an important influence on the powered case. 


Autorotation operation means no restraint of the rotor rotation about the 
hub. Therefore, no rotor torque is transmitted to the shaft, and no pylon roll 
motion is transmitted to the rotor. The collective lag degree of freedom is 
then rigid-body rotation of the entire rotor about the shaft. The collective 

lag mode shape is then n = v; with no restraint of this degree of freedom, 

£ 0 

v = 0. The lag motion of the blade is defined with respect to the pylon, 

which is rolled by a , rather than with respect to an inertia frame; there will 
then be a response of £q to pylon roll ql z , but the zero spring rate for col- 


lective lag motion (v 




0) assures that this motion will be only that 


required to hold the rotor fixed with respect to space. Since no torque is 
transmitted through the rotor shaft during autorotation, any Cq forcing terms 
in the pylon/wing equations of motion should be dropped. This would auto- 
matically be accomplished by the zero spring rate (v = 0, as in eq. (49)), 

C 0 

but actually dropping Cq would result in a better conditioned numerical prob- 
lem. With n = r , it follows that I* = I* = I * (and all are equal to 1 if 
So So So a ° 

= I 0 is used); and the zero collective lag spring rate means that v =0. 

So 

These changes are all that are required in the equations of motion to accommo- 
date the autorotation case. 


The Rotor Aerodynamic Coefficients 


In this section, the aerodynamic coefficients of the rotor forces 
defined in equations (40) are derived and discussed. The derivation follows 
that used for the four-degree-of-freedom model; all the coefficients required 
will now be derived, and the effects of drag and compressibility retained. 
First, the basic set of coefficients is considered; the M and H terms (out-of- 
plane and inplane forces) with subscripts y, 3, and 6 (inplane and out-of- 
plane velocities and blade pitch control), with T Q and Q 0 for the trim forces. 
The trim thrust and torque follow directly from integration of the section 
forces (no perturbation process required) as 
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From the perturbations about the trim state, the flap moment is 

J T fo dT '-j 2 £ - U p 2 ?)* 



6a 


'd 


Uu 


P 2a 


6a - Uu 


M 


P 2 a 


■] 


rt '(“r d - U P 25 ) 

“Y 


6M dr 


dr 



Me „ o 
l a U 1 M P U T 


la U 2a U 


Me 


2a U 




d M 

M 


2 a 

U 

2 

Me . 

u T 

i 


] 


dr 6tiy^ 


‘M P T 


U 2a U 2a U 

°_d( U p 2 \ U P U T %"/l, „ 

2a \ U + U ) + 2a U ~ 2a U J dv &Up B 



+ # r[Uu 

’o 


T 2 l-^P 


2 a ) 


dr 60 


= M o + M u Su T B + M ^ u P B + M q 


66 


and, for the blade drag force, 
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The radial force is simply 


42 





° d 

U ~ dr 

2a 


) 


Su 


7? 


8 — 

oa 


= R &u 
P R 


8 


r\ 

'll 

oa 


The coefficients may be identified now; the trim values of the blade 
section force coefficients and their derivatives, evaluated at the 

angle of attack and Mach number of the radical station) must be used in the 
integrands, and the trim velocities are 


u T = r> 

Up = V + v 
U = /r 2 + (V + y) 2 


For convenience, V + y will be written as V in the following expressions. If 
one substitutes for the trim velocities, the rotor aerodynamic coefficients are 
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(Eqs. (50) continued on next page.) 
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where U - (r 2 4 - y 2 ) 1 / 2 . The expression for //g* was obtained from 

tfg = Crp/oa + //£* using the result for Cf/oa. The effect of F^ in the sum 

+ R n is just a term (\?j/2 a)U in the integrand. So R. i: only adds to the drag 
terms of H^, which are usually negligible for high inflow' aerodynamics anyway. 
The following expressions for // and II q are also required, at least to find 
the other coefficients: 



The behavior of (and the expressions for) all the other coefficients can 
be inferred from a knowledge of the above set. The other coefficients follow 7 
immediately from the fact that, in the spanwise integration, the T terms have 
one less 2 1 than the M terms, the Q terms have one more r than the H terms, 
subscripts l have one more r than subscripts y, and subscripts X have one 
less /' than subscripts 0. The basic behavior of the coefficient is determined 
by whether it is an out-of-plane force or an inplane force, produced by an 
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out-of-plane or inplane velocity or blade pitch control. The combinations 
H \i + /? p and il$* = //g - C^/aa are special effects that do not have comparable 
forms in the other coefficients. There is no radial drag due to blade lagging 
velocity (no and no torque moment due to the radial drag force, so there 

are no R contributions to or any of the Q terms. The inplane force due to 
flapping velocity is written Hfe = Cy/oa + ffg* only to show explicitly the term 
due to tilt of the thrust vector, so it is not extended to H-y or any of the Q 
terms. Hence the coefficients and H g are used to derive the other 
coefficients, not and Hq*. 

If all terms but a £ and Co are dropped from these expressions (and the 
degrees of freedom reduced to 6 $ a and a x ) , the results of the 

previous section are recovered. 

Evaluation of the coefficients- Approximate expressions for the 
aerodynamic coefficients are obtained by evaluating the integrands at an effec- 
tive radius, with the techniques of the previous section. The results for the 
complete set of coefficients required are (dropping terms that are negligible 
for both high and low inflow) : 

M q = (3/4) (C T /aa) 

= (sin <t>)/6 + 2Cf/aa 

= (sin <j>)/8 + 3/2 {C T /aa) 

M g = -(cos <J>)/8 

M = -(cos <t>)/6 

= 1/ (8 cos <j>) 

- Crpf oa = a/ (6 cos cj> ) 

T = (sin <f>)/4 + 3C T /aa 

T- = (sin ) / 6 + 2 C^/oa ^2) 

Tg = -(cos <f>)/6 
= -(cos <j>)/4 
= 1/(6 cos <j>) 

H q = (4/3 )(C Q /aa) 

H v + R v = {V Sin ^ /2 ~ 47 ^T/^) + 6 (Cq/ aa) 

H;. = ( V sin <j>)/4 - (4/ 3) V (Cy/ aa) + (8/3) (C Q /aa) t 
H g = -(sin <t>)/6 + C T /oa 

(Eqs. (52) continued on next page.) 
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H. = -(sin <t>)/4 + (3/2 )(C T /aa) 

A 

H q = V/ (4 cos <(>) 

= C^/aa = 7(Cy/cra) + <?y8a 

= (7 sin <j>)/4 - (4/3)7(Cy/aa) + (8/3) (C^/aa) 

Q = (7 sin <t >)/6 - F(Cy/aa) + 2 Cq/oa 
Q = -(sin (j >)/8 + (3/4) (Cy/aa) 

P 

Q = -(sin 40/6 + Cj^/oa 

£? 0 = 7/(6 cos 4>) ; 

where the inflow angle is evaluated at an effective radius Vq {Vg = 0*75 
usually works well) so that 

sin <p = 7/Ag 2 + 72 
cos 4 > = ^g/ / jp , 2 + F 2 
o(r e ) = 0(r e ) - tan" 1 7/t\_, 


( 52 ) 


Wherever V occurs in these expressions, V + v is really meant, but the effect 
of the induced velocity V is important only for low inflow. The similar behav- 
ior of the coefficients within a group (as discussed above) appears in these 
approximations . 

The c iC [(Cq) terms are important only for the inplane forces due to inplane 
velocities at low inflow. The c^(Cy) terms are important for low V, but never 
for the out-of-plane forces due to out-of-plane velocities. All Cy and Cq 
terms may be dropped for high inflow, where V is of order 1 ; in that case, 

V + V may be replaced by V also. The compressibility influence, especially C£ M 
and may result in important contributions from the lift and drag terms 

even in high inflow; stalled flow may affect and . In general, however, 

for high inflow the behavior of the coefficients is given primarily by the c z a 
terms alone; this is expected to occur so long as significant stall or com- 
pressibility effects in the trim operating state are avoided. 

The basic features of the rotor aerodynamic coefficients in high inflow 
are obtained then the oi a terms. Of particular interest is the variation with 
forward speed V. If it is assumed that V/ttR is of order 1, that all the blade 
section coefficients are small compared with a £ , and that cijla - 1 / 2 , then 
the coefficients reduce to 
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Jv 2 + V 2 


dr 


r 2 V 

Jr 2 + V 2 


dr 


Vy/r 2 


+ V 2 dr 
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These coefficients were obtained in the discussion of the aerodynamics for the 
four-degree-of-freedom model, and the integrals were evaluated in that section 
(eqs. (31)). The variation of these coefficients with forward velocity is 
shown in figure 4. The behavior of any of the remaining coefficients is simi- 
lar to one in this set, but with slightly different numerical constants. 


If only the oi a terms are retained and with o £ /2a = 1/2, the entire set 
of aerodynamic coefficients is 


« - Vf 2 

«• . Vf 3 

M‘ = -f 
6 J 4 

M = -f 
X J 3 

M e = 32 

T = Vf. 
y l 

t- = vf 

Q J 2 

= -^3 

r x = ~fz 
t q = g \ 


H - V 2 f 
y J 0 


H- 


v 2 f 1 


H z = V - 

H x = 

H e = Vg o 
«y = v % 

% = v% 

% = - v ?3 

Q, = -Vf, 


-Vf, 


Qc 




(54) 
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where 


r n dr 


f = 


9. 


n 


/ 


1 


0 /p2 + y2 

1 


1. J r n v/r 2 + 7 2 dr 


The integrals required are 


A 

f l 

f 2 

f'3 

A 

A 

^2 


= ^ (/l + 7 2 - 7) 

- T '^ T1S - I v % 

= [/l + 7 2 (1 - 27 2 ) + 27 3 ] 

= ^ A + 7 2 (2 - 37 2 ) + | T 4 /^ 

= y A + l/ 2 + x 7 2 / 

4 2 J 0 

= i [(/l + 7 2 ) 3 - 7 3 ] 

= ± A * 7 2 (2 + 7 2 ) - | 7 4 / Q 


The aerodynamic coefficients can also be evaluated using the exact 
expressions (eqs. (50)), including the effects of the blade lift and drag, and 
of stall and compressibility - the only approximation then being in the knowl- 
edge of the section aerod>oiamics . To evaluate the coefficients requires the 
blade section force coefficients and their derivatives with respect to angle 
of attack and Mach number: o^Jla, o % /2a, a^/2a, a^j2a, e £ /2 a + Ma^Jla, 

and o^!2a + The coefficients are to be evaluated at the blade trim 

angle of attack and Mach number: 
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a = 0 - tan 1 (7 + v)/r 


M = 


M, . Jv 2 
tip 


(7 + v) 2 


The blade twist and the collective pitch required for a given flight condition 
are then required to evaluate the section aerodynamic, in contrast to the 
results with only <?£ terms where only 7 is required. The section coefficients 
may be evaluated from tables of airfoil section data, appropriate to the rotor 
blade being considered. 


Alternatively, representative analytic expressions may be used for the 
blade section coefficients. While a particular section may not be considered 
then, a representation that includes typical stall and compressibility effects 
will allow a general study of the influence of these effects on the rotor 
aeroelastic behavior. The following expressions for the section coefficients 
are used here. Below a stall angle of attack a s (typically 12°), a constant 
lift curve slope is assumed, with a Prandtl-Glauert correction for compressi- 
bility, so that eg = aa(l - M 2 ) -1 / 2 , and the lift coefficients required are 




M 2 )" 1 / 2 | 



M 2 y 1 /2 


(55) 



2a + 2 a 



M‘ 2) -3/2 j 


In practice, the Prandtl-Glauert factor (1 - M 2 ) -1 / 2 is truncated at its value 
at M = 0.95, say, to avoid numerical problems near M = 1. The drag for 
unstalled flow is 


ad = 0.0065 - 0.0216a + 0.4a 2 + A 

( 0.43 (M + | a | /0. 26 - 0.9) if |a| > a^ y 

0 otherwise 


0.26(0.9 - M) > 


(56a) 


This is the classical result of Bailey (ref. 31) with a lower as 

appropriate for current proprotor sections, and with a compressibility term 
obtained from section tests on rotating blades. The compressibility drag 
increment has a critical Mach number of 0.9 at zero angle of attack; above the 
critical Mach number, there is a large increase in ej. The drag coefficient 
and its derivatives must be divided by 2a; a = 5.7 is used for the two- 
dimensional lift curve slope. For stalled flow, |a| > a g , the following 
approximation is used: 
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°l = % s 3 n ^ 
a d = ^ s (sin a)2 

where, typically, c £ = 1.0 and ej = 2.0. Combined compressibility and stall 

effects are not included. A very important influence of stalled flow is that 
possibly c? £a = 0 there, so the c £ , cfa c^, and ty terms in the coeffi- 

cients dominate the behavior even in high inflow (if a large enough portion of 
the blade is stalled). 

With these analytical, semi empirical expressions for the section 
aerodynamics, the influence of the drag and lift terms and of stall and com- 
pressibility on the rotor aerodynamic coefficients is examined. For the 
design of a specific rotor and the prediction of its behavior, the section 
characteristics appropriate for the actual blade sections should be used. For 
the present work, it is desired only to check the relative importance of these 
effects so approximate lift and drag coefficients are satisfactory. The 
influence of these effects on the rotor coefficients is shown in figure 5, for 

and M$. The coefficients were calculated using the exact 
expressions (ecjs. (50)), with the above approximations for the section aero- 
dynamics, for two rotors. • The collective required to give the rotor thrust 
for equilibrium cruise at a given V is used. In figure 5, the results for 
these two rotors are compared with the coefficients found using only the <? £a 
terms. The coefficients with only' the c £a terms are given by equations (31) 
Cfig. 4); the approximation is the same for the two rotors since it is inde- 
pendent of the section characteristics. The most important difference between 
the two rotors so far as the behavior of the coefficients is concerned is that 
they have different tip speeds. Therefore for a given forward speed V/SIR, the 
blades have a different resultant Mach number M = M._- t (r 2 + y 2 ) 1 / 2 at a sec- 

tion. The tip resultant Mach number M = M+.p v (] + y 2 ) 1 / 2 is shown in figure 
5(b) for the two rotors. 

The exact coefficients in figure 5 show a significant difference from the 
coefficients based only on the terms; the difference is particularly large 
when the tip critical Mach number a (0.9 for a = 0 with the section characteris- 
tics used) is exceeded. The following conclusions are reached then: the o £ 

terms in the rotor aerodynamic coefficients give the basic behavior, at 01 
least so long as the section critical Mach number is not exceeded; the other 
terms in the coefficients are not negligible, however, and should be included 
to properly evaluate the behavior of a real rotor, especially when operating 
at high section a or M. When the section aerodynamics other than o £ are 
required, actual section characteristics should be used rather than a 
representative expressions . 

Three methods for evaluating the rotor aerodynamic coefficients have been 
described: 


(56b) 


Specifically, this is for the two full-scale rotors examined in later 
chapters, hence the labels "Bell" and "Boeing" in figure 5. 


50 



(a) Approximations based on evaluating the integrands at an 
equivalent radius . 

(b) Approximations based on the retention of only the ci a terms, 
with Qijla = 1 / 2 ; the integrals may be evaluated, giving 
the coefficients as functions of V alone. 

(c) Approximations based on analytical expressions for the blade 
section aerodynamics in the exact coefficients; representa- 
tive stall and compressibility effects are included, but no 
specific section is modeled. 

The coefficients based on the equivalent radius approximation (method (a)) are 
used only for the expository development, never in the calculations. Method 

(b) treats the o^ a terms correctly and exactly, and if the other terms are 
required method (c) should be used rather than method (a). The coefficients 
based on just the o £ terms (method (b)) will normally be used in the calcula- 
tions here. These coefficients include the basic behavior with inflow ratio, 
which is of primary interest here. In fact, this level is usually a good 
approximation for calculating the dynamic behavior (as shown later) . Method 

(c) , the coefficients based on the exact expressions (eqs. (50)), is used here 
only to check the influence of the terms other than . 

The method used here allows the derivation of rotor aerodynamic 
coefficients, including the influence of lift and drag and of stall and com- 
pressibility, with no more difficulty than a derivation that includes only 
terms. Therefore, a good representation of the rotor aerodynamics is available 
if one chooses to use it (and if enough information on the section aerodynamics 
is available). Evaluation of the coefficients in method (c) , or even including 
tabular data for the actual blade . sections used, requires numerical integration 
over the span, but that is no problem for numerical work. There is only one 
real complication in evaluating the coefficients by the exact expressions: the 

trim angle-of-attack distribution is required, which means that the blade col- 
lective pitch at the given operating state must be known. Hence a preliminary 
solution of the rotor performance to find the collective pitch is required 
before the coefficients can be evaluated for the dynamics. In contrast, with 
only the c ^ terms (method (b)), only VftiR is required to evaluate the 
coefficients . 

Discussion of the coefficients- Some properties of the coefficients are 
discussed here; in particular, certain useful equivalence among the coeffi- 
cients are derived. A helicopter rotor in hover (low inflow axial flight) 
exhibits equivalence of control plane, hub plane, and tip path plane tilt; 
that is, these inputs produce the same forces on the helicopter (with certain 
exceptions). This behavior translates into certain equalities among the rotor 
aerodynamic coefficients. The influence of high inflow operation on the 
behavior of the rotor, and on the basic set of coefficients in general is 
examined now. 

Consider the longitudinal moment on the rotor disk due to tip path plane 
tilt (Bic, 3x5) > hub plane tilt (oty, a x ) , and control plane tilt ( 0 j q, © 15 ); 
from equation (42), this moment: is 
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Hub plane tilt oq, gives an inplane component of V + y, hence a flap moment 
through the speed stability coefficient M p . Tip path plane tilt Bic (with 
respect to the hub plane) gives a flapping velocity in the rotating frame, 
hence a flap moment through the flap damping coefficient Mq. Control plane 
tilt Oj 7 produces a flap moment through the pitch control power coefficient 
Mq. For low inflow, 


"e - -«e = j 

so the flapping produced by blade pitch is 6 lc -/6 li? = M q /-M'q = 1 (if v„ = 1) . 
That is the familiar result of helicopter hover control (with a low inflow 
rotor): the tip path plane remains parallel to the control plane. In high 

inflow operation, however, Mq and -Mq are not equal. Based on the equivalent 
radius approximation, 


so that 




1 

8 cos <f> 


Mi 


COS (f> 
8 




-Mi 


1 

COS 2 <J> 


Pitch control power Mq increases with V, while the flap damping Afg decreases; 
the ratio then increases with V . For a V/QR up to 1 or so, there are no 
drastic changes in the magnitudes of the coefficients, but the effect is 
important. Based on just the terms, the coefficients are 



The integrand of is a factor u^/U 2 = r 2 /(r 2 + V 2 ) smaller than that of 

therefore, the ratio of the coefficients is of the order r e 2 /(r e 2 + V 2 ) = 

cos 2 <J> (as above) . This high inflow effect results from the fact that the 6a 

due to B is not the same as that due to 0 when the inflow angle <J> is large. 

The flapping term comes from 6up - r B so that 

u6u uv m 

5a - 60 — — = 69 - -=— B 

U 2 U 2 

which gives the additional factor Upp/U 2 in Afg. 
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Consider now the hub moment in terms of flapping with respect to space 
(in an inertia frame), that is, 


e lCj = 6 1 C ' a y 

01 5j = 015 + a x 


so that the longitudinal moment on the rotor disk due to flapping (with respect 
to space), hub plane tilt, and control plane tilt is 


M r 


15 


aa 


[-«• * ( V * v)M v ]<x y * «jC6 1Sl ♦ VlS 


Flapping with respect to space acts through Mg to give a hub moment as usual. 
The moment due to hub plane tilt is [-Mg + (7 + v)M y ], while that due to 
control plane tilt is Mq. In high inflow, these moments are not equal, for 
while both a and 0 tilt the control plane, only a y tilts the hub plane, 
reference frhme also. The difference is, in fact, small. From the definitions 
of the coefficients (eqs. (50)), it follows that 



Hence q 

* (V <■ V)M V ■ « e * ^ ( s7 > 

Now Mq = 1/(8 cos 4>) , so that Mq » Cq/oo. (which is of the order of VC T /aa ) for 
all V , both high and low inflow. So for all V y 

-M* + (V +v)M = 

6 n 0 

which means that hub plane tilt and control plane tilt are equivalent. For 
low inflow, ( V + v)M v is of order V 2 small, while Mg and Mq are of order 1 and 
both equal 1/8. 

Similar results may be obtained for the inplane forces. Consider the 
vertical hub force due to flapping, hub plane tilt, and control plane tilt: 
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oa 


= (V + y) (// 


/? )a 

C y 


V e !5 
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x 




+ Vi 


o 
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Tip path plane tilt 6-^ gives an inplane hub force by tilt of the thrust vec 
tor, while blade pitch control has no such effect. The corresponding hub 
forces due to flapping and pitch, II ^ and //q , are equal for low inflow where 


In high inflow, these forces, like the flap moments, are no longer equal. On 
the basis of the equivalent radius approximation. 


so that 


H = 

0 4 cos f 

//.* = - sin 
(3 6 


// 

“0 . 1 



COS 2 f 


Both forces increase with 7, but //q increases faster bv a factor 

-cos 40 ‘ - O'V + 7 2 )/?^ 2 . On the basis of just the an terms, the 
coefficients are a 




UUy 


:i7* 


Vp 


la 


The integrand of 


tfg* is a factor wpp/U 


p 2 / O 


) smaller than that of 


cq, which agrees with the result from the equivalent radius expressions. 

Also, it is the same factor found between the integrands of and Mr , as 
expected since the difference between the coefficients has the same source in 
both cases, namely, the difference between 6a due to 0 and due to (3 when the 
inflow <p is large. In terms of flapping with respect to space, the hub 
force is 
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The equivalence of hub plane tilt and control plane tilt may again be 
demonstrated; from the expressions for the coefficients, there follows exactly 


-#•* + (7 + v )(// y + V = H q + W + 
which is, for all 7, approximately 


-y + (V + + V = H e 

and if only high inflow is considered, so II' = H , it follows that 


+ (7 + vH H v + RJ = H q 

The speed stability coefficients and + /? y produce moments and 
forces due to hub inplane velocity - that is, due to shaft linear velocity, 
aerodynamic gusts, or angular velocity of the pylon about a pivot aft of the 
hub. For low inflow, these coefficients are an order 7 smaller than the other 
moments and forces. For high inflow, however, the speed stability coefficients 
are of the same order as the others, and this difference between the low inflow 
and high inflow aerodynamics may be expected to have an important influence on 
the aeroelast ic behavior. The speed stability coefficients usually pioduce 
positive damping of the motion (see, e.g., the four-degree-of-froedom model, 
eqs . (26)). With high inflow, hub plane tilt results in an inplane component 
of 7 that produces forces and moments through the speed stability coefficients. 
These forces are in the same direction as the displacement, so they act as 
negative springs (eqs. (26) again). The speed stability coefficients often 
appear in the equations of motion as coupling terms rather than as direct 
damping or spring terms. In such cases, the effect of the coefficients may be 
significant since their magnitude is of order 1 in high inflow, the overall 
effect may be stabilizing or destabilizing, depending on the coupled motion 
involved . 

The inplane forces on the hub due to tip path plane tilt have 
contributions from the direct tilt of the thrust vector (2 C'j’/occ term) and from 
the inflow (ffg term). From the equivalent radius approximation, the relative 
size of these two terms is 


2 C/oa 



a/ (3 cos j>) 
(sin <|>)/6 


3_ a 
2 7 cos 2 <J> 


Thus the ratio is of order a/7. In high inflow, C T /oa may be neglected since 
7 is of order 1, and also because the working angle of attack a of the prop- 
rotor is low. These terms are very important for flapping rotor dynamics. 

The response of the rotor to low-frequency excitation, say, control plane tilt 
or shaft motion, is tilt of the tip path plane. Through 2Cf/oa + H$* , the tip 
path plane tilt gives a hub inplane force that couples the rotor and vehicle 
motion or is used to control the vehicle. For high inflow, dominates the 
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thrust term, so the inplane force on the hub is not only an order larger but 
also of opposite sign. This effect is expected to produce dynamic behavior 
significantly different from that of low inflow rotors. Hub plane angular 
velocity and tip path plane angular velocity act through Ho to produce inplane 
iorces also. In high inflow, the thrust term may be neglected, so Ho = Ho*- 
the coefficient is then the same as for tip path plane tilt. ** P ’ 

Consider now the cross derivatives M and H», which are out-of-plane force 
due to inplane velocity and inplane force due to out-of-plane velocity respec- 
tively. On the basis of only the terms, these coefficients are equal: 


M U = -H, = 


f 

J 0 


r 2 V 


/r 2 + V 2 


2 a 


ar 


This equality is useful in simplifying the whirl flutter equations of motion 
m a ater chapter. Examine, therefore, whether it remains valid with the 
more complete expressions for the coefficients. From equations (50), 


M + H- 

a 6 




(59) 


The right-hand side is negligible for large V, compared with M u or Ho- 
nigh inflow, one has always then y 


for 


M 

P 



High inflow, V/SiR of order 1, influences the rotor aerodynamic 
coe lcients substantially. It follows then that the features of high inflow 
aerodynamics are an important factor in the aeroelastic behavior of the rotor 
and wing system. In summary, the combinations of the coefficients derived are 


-M- + {V + v )M = M \ 
6 a 0 1 


-//•* + ( V + y) (H + R ) = h 
p a a 0 


M = -//• 
a B 


I 


(60) 


The first two approximations are valid for all V, while the last 
high inflow; to these may also be added ffg = Hg* for high inflow 


is only for 


Performance considerations- The 
coefficients for the analysis of the 
tion of the proprotor performance. 


evaluation of the rotor aerodynamic 
proprotor dynamics requires a considera- 
First, to obtain the rotor collective 
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pitch value requires a solution for the performance in a specified operating 
state, for example, cruise flight (C T given by the airplane drag) or autorota- 
tion (Cq = 0). The rotor collective pitch is needed only to evaluate the com 
plcte expressions for the coefficients, however (eqs. (50), described as 
method (c) previously). Second, the total axial velocity V _ + v, a major param- 
eter in the coefficients, includes the induced inflow V, which is related to 
the rotor thrust and operating state. Two topics are now considered: the 

expressions required for an elementary analysis of the proprotor performance 
(Crp and Cq) to find the collective pitch; and an evaluation of the rotor- 
induced inflow. 


In the previous analysis, expressions were obtained (eqs. (50)) for the 
rotor thrust and torque coefficients in terms of the blade section aerody- 
namic forces. Using the identity of the power and torque coefficients for the 
rotor (C p = Cq since P = QSl) , these expressions are 


°T 

aa 

C _P 

aa 


/: ■[ 


>• Ta - ^ ^ Ta 




H 

J dr 
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where l/ 2 = r 2 - + (V + t?) . If uniform induced inflow v , is assumed, the power 
expression is 

• 1 


C p = (V + v)C T + . 2 

/ 0 


/ 


ac 


d U 3 dv 


(62) 


which is the usual result for the power required by a rotor operating in axial 
flow. The first term is the sum of the induced power loss and the useful work 

done : 


C p = (V + V)C T = (V + C t /2V)C t (63 

i 

For high inflow, the induced loss term vCj, = Cf 2 /2V is negligible (see 
discussion below); then Cp^ is proprotional to Cf, in contrast to hover 
operation where the induced inflow is important and for which Cp ^ is 
proportional to Cy 3 / 2 . 

The second term in equation (62) is the rotor profile power loss; if the 
section drag coefficient is constant over the span (or if an effective mean 
value of Oj is used), then 
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(6.27) 

V + v = 1 


oo 

ct 

8 
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high V 

operating conditions of 

primary interest here are: 
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operation 

‘•" c F^vruing trie propulsive force required for equilibrium cruise 

light, and autorotation operation (no net power supplied to the rotor through 
the shaft). The latter corresponds to the condition in which dynamic tests of 
the proprotor and wing arc often performed. The performance problem involved 
is then to find the rotor collective pitch required (and the other coefficient) 
tor a given Cf (thrust required for equilibrium cruise flight) or Cn ( 7 ero for 
autorotation) . H 

The rotor aerodynamic coefficients require an estimate of the rotor 
induced velocity V . For high inflow, however, the induced velocity is much 
less than the forward speed of the rotor (as shown below), so great attention 
to the inflow calculation is not required to satisfactorily evaluate the coef- 
ficients. The assumption of uniform inflow is adequate then, and it may, in 
fact, be possible to neglect the induced inflow entirely. The rotor Ct 
required in cruise flight is obtained by equating the rotor thrust (for two 
rotors) with the aircraft drag, and expressing the drag in terms of an 
equivalent flat-plate area f for the aircraft (D - (l/2)p V 2 f). Therefore, for 
the thrust required of the rotor. 


C„ 


- L ~L r/2 
4 A ■' 


where A is the disk area of the rotor. Momentum theory can be used to estimate 
the induced inflow; the usual result (ref. 32) for axial flc 
thrust Cm is: 


low operation at 


V + v = 7/2 + J~[V / 2 ) 2 + C T /2 


(65) 
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Substituting for Cf yields 


y _ A 4 C//2.4) - l 
V ~ 2 

or, since f/A and hence Cf/V 2 is small, 

V ~ S’ _ I £ 


Typi cal values of proprotor aircraft drag and radius give j/A O.OO.vi, it 
follows then that v/V = 0.0017 typically, which is indeed much less than 1. 
for the proprotor in equilibrium cruise flight, neglecting i /. is a rtusona 
approximation then. 


When the rotor is operated in autorotation, the performance 
is that Cq = 0, which requires that Cf = -Cn / ( v + A = ~ c QoA' ■ 
inflow result for the profile power, Cq q 2 (crcj/2) 7 > the thrust 

autorotation is 


requirement 
With the high 
required in 


For the proprotor in axial flight at high in flow autorotation occurs in the 
windmill brake state fi.e., at V > 2/\C T \/2, or \C T /2V 2 \ 2 \v/V\ < 1/4); hence 
momentum theory may again be used to estimate the induced inflow (ret. 

The same conventions are used for the directions of V and Cf (so the Cf 
required in autorotation is negative, as given previously); the momentum theory 

result is 


V + v 




c 


( 66 ) 


Substituting for the required Cf yields 


y 

V 




Typical values of proprotor drag coefficient and solidity yield 

oe j /2 = 0.0004. Hence the typical induced inflow v/V = -0.0002, which is 

smaller than in cruise flight. Again, the induced inflow may be neglected. 

The effect of the induced inflow on the blade load distribution (that 
required to evaluate the aerodynamic coefficients for the dynamics analysis), 
may be investigated by considering the change in angle of attack due to v/V: 


6a 


6 tan 1 


V + V 


v_ T/Elr 
V 1 + (V/ttr~) 2 


which has a maximum value of 6a = (l/2)y/7. Use of the momentum theory 
result v/V 2 Cf/ 2k 2 (valid for both powered and autorotation operation in 
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high inflow) and the mean angle of attack a 
in angle of attack: 


6 Cjt/aa yield the maximum change 


6a 

a 


For typical values of rotor solidity, 
for high inflow ( V of order 1). 


. gq 
24 V 2 

Sa/a = 0.02V' 2 , 


which is a small fraction 


Numerical calculations were performed to verify that V + v = V is a good 
approximation in the calculation of the aerodynamic’ coefficients . The impor- 
tant consideration in the dynamics analysis is that the performance calcula- 
tion and the calculation of the aerodynamic coefficients be consistent, either 
neglecting the induced inflow v/V or using the same estimate of v/V for both 
calcuiations Any error in estimating the rotor performance or the collective 
pitch required is not relevant to the dynamics analysis. The aerodynamic coef- 
lcients that correspond to the operation of the rotor at a given value of Cw 

or C Q are obtained, the only error being a small change in the angle-of-attack 
distribution over the blade. 


SECTION 3: BEHAVIOR OF ROTORS IN HIGH INFLOW 


„ " the next four chapters, several topics on the behavior of high inflow 

1 p otors are investigated, based on the equations of motion derived pre- 
viously. The development of the proprotor and cantilever wing model is 
resumed in section 4. The reader interested in that topic may skip the four 
chapters m this section. P 


Elementary Dynamic Behavior 


Some aspects of the dynamic behavior typical of proprotor aircraft are 
examined First, the fundamental stability of the blade motion is examined 
through the eigenvalues of the uncoupled blade motion. Then the influence of 
the transformation to nonrotating degrees of freedom and equations of motion 
on the eigenvalues of the rotor is examined. The actual coupled motion of the 
proprotor and wing system is considerably more complex, but these considera- 
tions are useful m the interpretation of the results for the complete model, 
he equations of motion and the hub forces for the rotor were found to separate 
into longitudinal and lateral/vertical groups Ceqs. f44) to (48)). The wing 
motion in f a c tj couples these groups, but it is useful for a preliminary 
s udy of the dynamics to neglect that coupling, and examine the rotor response 
to shaft motion, gust, or blade pitch control in the longitudinal or in the 
lateral /vertical systems. Attention is directed to the low-frequency response 

t^ c :i^ft e :Sin;:. useful in evaiuating the infi — ° f the - ’ 
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Blade stability- Consider the uncoupled, shaft fixed flapping motion of a 
single blade" The equation of motion ( in the rotating frame) is 

( 67 ) 


Ip* 6 - Y^B + (I 


8 0 


+ £ p YM Q )0 = yM 0 6 


The roots (eigenvalues) of this equation are 


yM- 


A = 


21 * 


± % 



+ 


YW„ 


■6 


(_ 3 ) 

WbV 


( 68 ) 


For low inflow -Mo = Mp, = 1/8, and equation (68) reduces then to the usual 
result for the flapping motion of a hovering rotor. The flap damping 1S P°si 
five. Mo < 0, so the real part of A is negative and the flap motion is stable 
As increases, the flap damping -M 0 decreases and the pitch control power 

Mq increases . 


ow 


and* the stability of the flapping motion decreases. The change is not great 
for V/9.R of order 1, however; and Mg is always negative (at leas e i 
contribution is) so the motion remains stable for even extremely high in 

Pitch/flap coupling K P introduces a flap spring term through the 
aerodynamic force M e , which changes the effective flap natural frequency. 


V + K P 


yM 0 


(69) 




Negative pitch/flap coupling, K P > 0, increases the effective flap frequency 
v fi . Increasing V/UR increases M 0 , and so increases the effectiveness of Kp. 
Again, the influence is not great for V/SIR of order 1 or less. 

Consider the uncoupled lag motion with the rotating equation of motion 
(homogeneous form) : 


i *i + y q-jA + = 0 


(70) 


The roots for this motion are 


A = 



(71) 


The lag damping Qj- is positive, Q £ > 0, so the motion is stable. e c e 
cient Qx increases with V/QR, hence the stability of the lag motion _ increases . 
For low inflow, the lag aerodynamic damping is very low. For high inf 1 , 

however, there is a significant increase in lag damping, which is important 
proprotor dynamics. In high inflow, the source of lag damping is the same as 
for the flap damping, namely, the lift change due to angle-of-attuck 
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perturbations (i.e. , the terms). Therefore, both flap damping and lag 
damping are of the same order in high inflow. 8 

The flap and lag motions of the blade are strongly couplc-d by the 
aerodynamic forces in high inflow; specifically, they are coupled' bv the cross 

rrssir * and * ,<hich aro of ° rder 1 

Nonrotatjn^_s^s t cm eigenvalues- Consider the blade uncoupled flap motion 

n °r° r tat \" s fr r- Th ° coni " E ^ s » *** 

motion Ceqs (47)) identical to that of the blade in the rotating frame so 
ttie eioenvalues of its motion are the same as those given previously. The 

mi! l i' 1 ° ^ i ° j 71011011 for t,ie t;i P path plane tilt coordinates, 6i r and 8i ~ are 

modified by the transformation to the nonrotating frame; centrifugal and’ 

f^moJiois™ 5 are lntr ° duCed > Which have the ef ^ct of coupling the 8 1C and 


The homogeneous equations of 
( from eqs . (44) ) : 


motion in Laplace form for 8 and 8 1C art 

4- I O 


*s 2 -yM^s * C 2 - ! ) +K p yM Q 


21 o *S~ yM- 


' (2-fg*s-yA/g) 


I B *s 2 -YM-s + I^Uv^-l) + K pr /.? f 


0 (72) 


The characteristic equation for this system is the fourth-order polynomial: 

A = [ V a2 ■ Y V + V (v 8 2 ' 1] + V A/ 0 ]2+ (2 V A " yM b )2 = 0 f73) 

of 6 tbo g n ValUeS are . then the solutions of this polynomial; writing the roots 
of the flapping motion m the rotating frame (eq. (68)) as A = \ p Ap the 

nomotatmg system eigenvalues (roots of eq. (73)) are then R R> 

X NH X R ± Z (74) 

and the conjugates. The nonrotating eigenvalues have the same real parts as 
the rotating roots, but the frequency is greater or less than that of the 
otatmg roots by 1/rev. On the root locus plane (TmA vs. Re A), the nonrotat- 
ing roots are then shifted vertically by ±l/rev from the rotating roots. 

trinle°ooles h at ee A blad f ^ rotor, the rotating roots on the A-plane are a pair of 
at l a P n -f and 1 ^ S con J u S a te. The nonrotating roots are a single pair 

at x R and its conjugate for the 8 0 mode, plus a pair at A* + f and A R -i fnd 

their conjugates for the 8^ and 6 1£ . modes. In both the rotating 2nd non- 
mg rame, there are then six roots, corresponding to the three blades 
"on Tvn differential equation. "ihe eigenvectors of ^ 

raises mirr ! 6 " UP j : ’ l0 ' 10: ; of and 6 1.1 corresponding to the eigen- 
Khi- e are two modes (and their conjugates); each eigenvector defines 
a wobbling motion of the tip path plane, at frequency Im\ R - 1/rev for the 
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1 ow-freauency mode and at frequency Im\ R + 1/rev for the high-frequency mode 
The high-frequency mode wobbles in the same direction as the rotor rotatl ™j\ 
direction, and so is called a progressive mode. The low- frequency mode wot - 
opposite the rotor direction if Im\ R > 1/rev, so it is a regressive mode. If 
Im\ R < l/rev, however, the low-frequency mode wobbles in the same direction 
therotor rotation and so is another progressive mode. The two modes are 
denoted here by B ± 1 for the flap modes at frequencies ImX H ± 1/rev, 
respectively. 

The behavior of the roots of the lag motion in the nonrotating frame is 
similar to that of the flap motion (the influence of the transformation from 
the rotating to nonrotating frame is a general result for all rotor degrees of 
freedom). The coupled motion of and 5^5 is a low frequency an h 
frequency mode, denoted by 5 ± 1. each a whirling of the net rotor center of 
gravity about the shaft. This motion has an important role i n proprot o 
dynamics. The blade flap and lag motions are, of course, highly coupled, t 
each other and with the shaft degrees of freedom The basic flap and lag roots 
remain identifiable even for the coupled motion however, both by the _ frequ , 
of the eigenvalue and by the participation of the degrees o , 

eigenvectors. The characteristic location of the nonrotating roots, at t 
low-frequency and high-frequency rotor modes, remains an important featu 

the dynamics. 

Longitudinal system - Consider the longitudinal dynamics of the rotor,^ 
descriFed by equations 747) and (48). It consists of two degrees of freedom 
(coning Bn and collective lag S 0 ) excited by longitudinal shaft motion (hub 
velocity z P and roll angle o 3 ), longitudinal gust (Vug) , and rotor j-° | 

pitch control (6 0 ). The influence of the rotor is transmitted to the shaft as 
a thrust force and a torque moment (Cy and Cq) acting on the hub. 

The equations of motion for the longitudinal system are, in transfer 
function form (and matrix notation). 


X 


* s 2 -yM'S+I* v 2 nV' f 

Bo B Bo Bo * e 

-yQ^s+yK p Q Q 


yW-. 

I* s 2 +yQ's+Il v 


Co 


* v 2 

C~ Co Co 



Y 

Y Q, 


0n + 


yVM 

/ 

Y VQ, 


U G + 


y M'S 

I* 3 2 +yQl* 
Co“ C_ 


a + 


■S* s 2 +y M.a 


Bo 
Y Q 


r 


'p 
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The thrust and torque are 
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c. 


aa ~ Vo + VT y u n + Tis<x_ + 


M, * 


rff 


C 3 


'b 2 


s 2 + r,s 3. 


y (s e 1,2 • Y V * W 6 o - s >S 


c 

J? _ /* 


era 


5 0 * 


( 76 ) 


To simplify the notation, subscript 0 on the inertias and frequencies is 

lonHtilr",!^ reSt .° f thc . discussion of the longitudinal system; only the 
g ma ynamics are involved, so the values appropriate to the conine 
and collective lag modes are implied throughout. 8 

in a reLce e ti° W ' freqUCnCy reSP ° nSe ’ the e ^tions of motion to lowest order 


X *V 2 
8 6 


Y V ? 0 


I *v 2 
5 ? _ 



0 A + 


The rotor thrust and torque are 



( y U G + 3 p ) + 



(77) 


C T 

oa 


°R 

oa 


= T 0 


e 9 o + t x (Vu g + V + V 3 - Vi 


0 P O 


J *v 2 




/ 

The solution for the low-frequency response of the rotor coning mode to 
collective pitch, longitudinal gusts, and shaft motion is then 


(78) 


8 


0 



fyw 0 0 o + y M x (Vu c + z p ) 


+ yM • a 1 

t 3 J 


(79) 


The response of the 
the rotor dynamics, 
Cf/aa (eqs. (78)): 


rotor thrust to these inputs, including the influence of 
is obtained by substituting for 6 0 in the expression for 
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c„ 


= I T n - M„ 


K y T 
P Y 0 


aa~\ Q 9 j * v 2 + KpyM& 


+ lF X~ M X 


+ I T- - M 


K y T 
P Y 0 


I *v 2 + K-xM. 

0 6 P 6; 


\(Vu G + ^p) 


V r 0 


c 5 VV + V^' 2 

This result can be simplified by use of the relations. 


( 80 ) 


MT* - APi 


1 sin <t> sin 4> 1_ 


0 5 ' 'T0 8 cos 


8 6 cos <t> 


= 0 


and, similarly, M Q T X - M X T Q = 0. With these relations, the expression for 
thrust reduces to 


C V o* 

3§ * IVo * t x (Vu g * V * J 

0 . 


(81) 


The factor in brackets is the direct response of the rotor thrust to the 
inputs. The only influence of the rotor flap dynamics (Bg term in cqs. (. )J 
is to reduce (if Kp > 0) the direct low-frequency response by the ratio 
v 2/vo . The effective flap frequency v« includes the aerodynamic spring due 
to pitch/flap coupling Kp. The ratio vg/v^ = 1 if K P = 0 and, for values o 
Vg 0 and Kp typical of proprotor coning mode, it will always be nearly 1. 

With a procedure similar to that used for the rotor thrust, the flap 
dynamics are eliminated from the low-frequency lag response to give 


1 


J *v 2 

1 


[y<2 e 0 o + vQ\( Vu n + z p^ + Y ^c a s Y ^p^0^o^ 


C 3 


— [v« e 0„ + y^ Vu a * V * 1 


(82) 


V V C 


B 


Again, the only influence of the flap dynamics on the direct response is the 
factor v 6 2 /vg . The rotor torque follows directly from the solution for 

The low-frequency response of the rotor coning and collective lag motion 
to collective pitch, longitudinal gust, and shaft motion is then 
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i mP °^ tan !- term in the res P° nse is shaft torque due to longitudinal 
g st. The direct response has the magnitude VQ X = -(v sin cj>)/ 6 . The aerodv- 

r^ha °o r n- % h-T reaSe f greatl >;. with inflow behavior of Q x is similar 

to that of h 3 which is shown m fig. 4; in fact, with only the o\ terms, 

Q x ~ exactly) - the source of a fundamental problem of the higliinflow 

drive tral^ln Y- l0 ^ gl ^. udl na ^ gusts, namely, a great increase in shaft and 
drive t id In loads with t orward speed. 

tny . rotor aerodynamic coefficients ? x and Q‘ r result in hub thrust and 

m d ue crms that damp the Sp and a„ motions, respectively. Since 

• it . cos * d a J ld sdn 4>)/6, the damping of r,p decreases slightly as V 

increases (its behavior is like that of «y, while the damping of c increases 
fnnoi V " Uts behavior is like that of (/„). The thrSst force"duc to 

low is a usual feature of rotor dynamics; the influence of the torque due 
to rotational speed <?£ is much different than for low inflow rotors, however 
because the coefficient is orders of magnitude larger in high inflow. 

!n the derivation ol the low-frequency results, it was assumed that 
C 0 * °’ for the autorotation case, however, = 0. Then C Q = o always; and 

the characteristic equation has a factor a, that is, one root at the origin 
which indicates that the proper degree of freedom is ,; 0 , the rotor speed per- 
turbation. The other root of the collective lag mode is at A = -yQ t /2I* 

(uncoupled), that is, a real root with time constant t = 21* /yQ- = 12/yV sin <+> 
The low-frequency response of ? is then C 0 ? -in <p. 



There is a response of to shaft roll even though 

been defined with respect to the rolled shaft position; 
that c Q - a 3 , the rotor speed perturbation with respect 


~ because has 
j 0 0 

the response is such 

to space, is exactly 


zero. 


From the expressions for the coefficients with only the a z terms, it 
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follows that VQ X /Qi = -V 2 f 2 /V z f 2 = -1 exactly (eqs. (54)); 
of the rotor speed perturbation to a longitudinal gust (.a 


hence the response 
low freciuencv) is 


JL = _l (86) 

U G 

This is the response required to keep the inflow ratio constant including the 
velocity and rotational speed perturbations, as shown in the following 
analysis : 

/ _v_ \ = &v vnm 

\QR J SIR 2 

_V_ T6F _ 6(f Iff) 1 

“ OR L V SIR J 


SR ( U G + C 0- 


Hence &(V/SIR) = 0 implies i 0 /uq = -1 as given above. Ibis result holds e\cn 
for the complete model of the rotor and cantilever wing. 

Lateral/vertical system: Flap response - Consider the iaterai/vertical^ ^ 

group of degrees of freedom, inputs, and hub forces (eqs. (44) to (46)). Ihc 
cyclic lag degrees of freedom, S lC and are dropped to obtain a manageable 

set of equations. The rotor motion is then described by two degrees of free- 
dom (tip path plane tilt & lC and B 15 ) excited by lateral and vertical shaft 
motions (shaft tilt, a* and a y and hub inplane velocity, xp and y P ) , aerody- 
namic gust (angle of attack a G and sideslip B&) . and cyclic pitch contro ( lC 
and 0i e). The influence of the rotor is transmitted to the sha.t as 10 
drag force and side force {C H and Cy) and pitch and yaw moments {C My and C Mx ) 

acting on the hub. 

The low-frequency response is of interest here so that the basic response 
of the rotor in high inflow can be examined, and also as a basis tor investi- 
gating (in a later section) the influence of the rotor forces on the aircraft 
stability and control characteristics. For stability and control problems, it 
is convenient to obtain the equations of motion m a body-axis system rather 
than in the inertia axis system used so far (the inertia axes are most conveni 
ent for the cantilever wing problem, which is the mam subject of is p 
In body axes, the trim velocity of the rotor (7 + v) remains axial during 
shaft tilt by and a v ; then converting to body axes requires that the terms 
in the equations of mofion and hub forces due to the inplane component of 
(V + v ) produced by the shaft tilt (i.e., the (7 + v)a x and (V + v)a y terms) 
be dropped. The shaft angular velocity a x and d must then involve a centrif- 
ugal acceleration in order to turn the vector (V s + v) to keep it aligned with 
the shaft. This may be accounted for by simply adding the rotor mass to e 
aircraft mass; hence the equations presented here do not incorporate these 
centrifugal acceleration terms. The perturbations to the aircraft velocity 
( v + v) are given by the hub velocities x p , y P , and a P . These motions and 
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ane velocity of the hub in body 


aerodynamic gusts are the only sources of inpl 
axes. 


The equations of motion for B 1C and Q 1S (in transfer function form with 
the a x and ay terms dropped as discussed previously) are 


iys 2 -Y«-s.i *(v 2 -i )+ x pY m 


b b 

-( 2I £*s-yM-) 


+ y VM 


2 V°-y“k 
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velocity °V he r ? t0r t0 CyCliC C ° ntrol > inplane hub 

iclocity, and shaft angular velocity. For the low-frequency response only 

the lowest order terms in s are required and the equations become ’ > 
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V f V- 1)+ V w e 


yM- 


~ yM k 

yc v e 2 -D + Y^ P . 
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IS/ 



+ Y M 


y p *n G -hi x 


- y M- 


- 21 * 
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-21* 
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Y M* 
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Only the steady-state terms in the flap response are retained. A steady tip 
path plane tilt gives a flapping velocity in the rotating frame so there are 
moments through HfJ; since shaft angular velocity involves a Coriolis accelera- 
ion of the Rotating blades in the hub plane (S 1C and e 1<7 are flapping with 
respect to the hub plane), inertia terms also appear. 

Inverting the matrix on the left-hand side of equation (88) gives the 
solution tor $ and J M tne 

1C lo 



where 


*■' * 


J E ^ v e 2 - ^ + V A; e 


r 9 oi 


The parameter /■/* is a measure of the lateral/vertical coupling of the rotor 
response (lateral/longitudinal for helicopter orientation); A ? * = 0 for vg = 1 
and Kp = 0, that is, for an articulated rotor with no hinge spring or offset, 
and no pitch/flap coupling. Pron equation (69), it also follows that can 
he written N* - / . ' ( v p - l)/-yA/g, so the parameter is a mcasuie of the 
(effective) structural restraint oi the blade flapping motion. 


In the flapping response (eq. (89)), the bracketed quantities give. the 
basic effectiveness of each kind oi excitation, while the factoi involving a * 
accounts for the rotor flapping dynamics. Inside the brackets, most of the 
lock number factors cancelled, indicating that the response is mainly a bal- 
ance of aerodynamic forces; the exception is the third term, which is a bal- 
ance of the inertial forces (21* ) due' to shaft angular velocity T,- and.-.qy and 
aerodynamic forces ( -yA/p) due to tip path plane tilt Spy and Bp:- Cyclic 
pitch control , hub inplane velocity, and hub plane angular velocity result in 
aerodynamic moments on the disk. These moments cause the rotor to flap until 
the moment on the disk due to the- flapping cancels the appl ied moment , thus 
achieving equilibrium deflection. The coefficient "g is the effectiveness of 
flapping in producing a moment or the rotor disk, hence its appearance in the 
denominators always. Blade cyclic pitch produces an aerodynamic moment on the 
disk through Mq, so the net control effectiveness is M Q /-Mh . .. Hub inplane 
velocity produces an aerodynamic moment on the disk through with a net 

effectiveness of With the present model, there are three sources of 

hub inplane velocity; shaft velocity x r and i/ p , vertical and lateral aerody- 
namic gusts a r and (which are angles, so Va G and VB C are the corresponding 

velocity perturbations), and angular velocity of the. shaft about a pivot a 
distance h aft of the hub producing hub velocities ha x and da,,. 


Consider the response to shaft angular velocity and op,. The inplane 
velocity due to angular velocity about a pivot aft of the hub" gives a moment 
on the disk (hence flapping) through M u as discussed previously. The third 
term in equation (89), with effectiveness 2T| a /-yA^g_, is the lag of the tip 
path plane required to precess the rotor to follow the shaft (Bpy and 8jy are 
flapping with respect to the shaft), for the tip path plane to. follow the 
angular velocity of the hub plane requires a Coriolis acceleration of the disk 
of magnitude 2I$ a ct,, because of the rotation of the blades. A moment is 
required to produce this acceleration, that is, to precess the rotor disk, and 
that moment is supplied by aerodynamic forces on the disk that result from the 
flapping velocity of the rotating blades. The tip path plane tilts back, lag- 
ging the shaft tilt, until the aerodynamic moment due to flapping (of magni- 
tude -yWaBic) is just large enough to provide the required acceleration. 

There is a 90° azimuthal phase difference in the response, characteristic of 
rotor or gyro dynamics: shaft pitch a y requires a lateral moment on the disk 
to precess it; that lateral moment is supplied by the flapping velocity in the 
rotating frame due to steady tip path plane pitch 6 ^ ^ • lhe flapping dynamics 
also introduce a phase shift if N * f (). 
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The last term in equation (89) is also a flapping response due to shaft 
angular velocity. Hub plane angular velocity due to a or d x produces a flap- 
ping velocity of the rotating blade, hence a moment on the disk through 
The rotor flaps until the moment on the disk due to tip path plane tilt, which 
also acts through -yMg, is just large enough to cancel the exciting moment. 

Both moments in the equilibrium are aerodynamic,, so the Lock numbers cancel. 
Both are due to the flapwise velocity of the blade so the coefficients Mo also 
cancel. The net effectiveness, then, is exactly 1. The last term in equa- 
tion (89) results in a coupling of the lateral and vertical response of the 
rotor to shaft angular velocity (ot x and a ), even if N* = 0. The basic flap- 
ping dynamics in the nonrotating frame (as discussed previously) consists of 
low-frequency and high-frequency flapping modes (8 + 1), each of which appears 
as a wobble in the tip path plane. The basic response of the rotor then 
involves coupling of the lateral and vertical motion, a fundamental character- 
lstuc thar results from the rotation of the blades. The low-frequency response 
should involve only excitation of the low-frequency flap mode (8 - 1) but the 
lateral/vertical coupling is still expected. 


When N* - 0, the first factor in equation (89), representing the influence 
ol the rotor flap dynamics, reduces to a unity matrix. With a hingeless rotor 
(vg > 1) or pitch/flap coupling ( Kp £ 0), N A is no longer zero and the flap 
ynamics matrix introduces a phase shift and a magnitude change in the flapping 
response. An example is the response to cyclic pitch control, Qiq and 0iy. 

The diagonal terms in the flap dynamics matrix produce a tip path plane tilt 
following the control. plane tilt, through equilibrium of the aerodynamic 
moments due to blade pitch and flapping. Hie off-diagonal terms correspond to 
lateral/vertical coupling of the tip path plane/control plane tilt, which 
results from the aerodynamic pitch moment being reacted by a moment due to 
structural flap restraint (v e > 1) or pitch/flap coupling (Kp ? 0), hence due 

1,16 two sour ces of tip path plane t ilt produced by control plane 
i t (or any other excitation) have a 90 azimuthal phase difference because 
the diagonal terms are due to flapping velocity (in the rotating frame), while 
the off-diagonal terms are due to flapping displacement. 

Consider the individual terms in the flap responses - each has the form 
of a basic control power multiplied by a factor (involving N,) that accounts 
for the flap dynamics. The effectiveness of cyclic pitch is given by (using 
the approximations based on the inflow at an effective radius) 


0 . 
■M 


1 

cos 2 <{> 


1 + tan 2 cf> = 



(91) 


The pitch control power increases with V then, due to both an increase in the 
pitch moment M Q and a decrease in the flap damping M' R (as shown in fig. 4) • 
the increase is slow for the usual range of V (up to 1 or so). With the 
influence of the flap dynamics, the response to control is 
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( 92 ) 


90 1C 

V' M 6 

"is 

1 + N 2 

90 1C 


"is 

1 + N* 2 


and similarly for 8^. The total flapping response to cyclic pitch then has a 
magnitude 


iij! = ( 93 ) 

3,01 A + n/ 


and an azimuthal phase shift of 


Aip = tan 1 N * 


(94) 


Then if N* t 0, the magnitude of the response is always decreased (by a factor 
(1 + ff* 2 ) -1 / 2 ) whether is positive or negative (vg 2 - 1 is always positive, 
but the pitch flap coupling Kp can be positive or negative); and N * produces a 
phase shift of the response, the sign of which depends on the sign of N The 
coupling parameter N * is approximately 
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so N, increases somewhat with the inflow ratio V. The effect of the lap 
dynamics (N*) on the response to all the inputs follows that on the cyclic 
pitch (a magnitude decrease and an azimuthal phase shift), so, for the 
remaining terms, only the basic control power is examined. 

The response to inplane hub velocity is given by 
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The flapping due to hub inplane velocity increases substantially with V then 
primarily because of the increase in the speed stability but also because 
of the decreased flap damping Mg (fig. 4). This important influence of the 
high inflow aerodynamics produces significant differences in the dynamics from 
that of low inflow rotors. The response to aerodynamic gust is 
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This is the response to angle-of-attack changes rather than just velocity 
perturbations, hence the extra factor of V (though both are hub inplane X 

nrid. Clt H k £ 8 ! S) 'i 7116 substantial increase in high inflow of the ? flapping 
produced by hub velocity results in problems characteristic of the proprotor 
ireraft: high flapping or blade loads in maneuvers and gusts and substantial 

coupling of the rotor dynamics with the shaft motion. 


The flapping response to shaft angular velocity a,, is (for N = 

y * 
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Sout e thS d n?v^ iS A dU t t0 t ? e i Ub inplane velocity with shaft angular motion 
viouslv'l P Th P f? ^ ance .^ a 5 t ° f the hub (its behavior was discussed pre- 

precess the^i.k to f ?T m required t0 P roduce a moment to 

L th , k . t0 P°ii ow the hub plane; it increases somewhat with V because 

of the decrease m the flap damping H • . The response to d x is 
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^H eq H ilibr r ^ tHe m °I nent due t0 fla PPing velocity due to a x and g ]r , 
which does not change with 7. For values of the Lock number, mst height and 

bution to^h^fl I' kl 7) typiCa J 0f P ro P r °tor operation, the primary contri- 
bution to the flapping response due to angular velocity of the shaft is the 
term required to precess the disk: 


9 I 6 

3d 


96 


1C 


16 


y 


9a^ y cos 


( 100 ) 


With oniy the inertia term of the response to shaft tilt, the influence of N 
is the same as for the response to cyclic pitch. * 


Lateral/vertical system: Hub 


vertical hub forces and moments: 


forces and moments- Consider the lateral/ 
rotor drag force C H , side force Cy, hub 


C ^’ a u d hUb yaw m ° ment ° M x- Retaining only the flap degrees of 

convert’ fn^nd 1118 ^ termS fgiving in P lane components of V + v) to 

convert to body axes, and then retaining only the lowest order terms in s 

yields the following hub forces: ’ 
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The hub moments are simply 



If one substitutes for the low- frequency flapping response, the hub forces are 
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Basically, the response for each input is a direct H force plus a contribution 
from the rotor flapping motion; the last term is the aerodynamic force pro- 
duced by the flapping velocity due to shaft angular velocity, combining the 
direct and flapping terms. Note that the hub force response is considerably 
simplified if N* = 0; then the tip path plane tilt produces hub forces through 
only 2Cj>/oa + and the matrix in the last term reduces to simply Cf/oa 

(which may then be neglected for high inflow operation) . The hub moment is 
simply a constant times the flap response. 


Pivot moment - The total moment about the pivot a distance h aft of the 
hub is given by the hub moment plus the hub forces acting on the arm h: 
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The hub moment is directly proportional to the tip path plane tilt, $\q and 
3i s> hub forces (eq. (101)) are also produced by the tip path plane tilt 
through ICrpfoa + tfg*, the first term being the thrust vector tilt and the 
second, the negative H force contribution of high inflow. (The off-diagonal 
terms that result from pitch/flap coupling are not considered now, that is, 
it is assumed that Kp = 0.) The total moments about the pivot due to the tip 
path plane tilt are then (after the hub moment and hub force terms are 
combined) : 



(105) 


The combination of terms appearing on the right-hand side has already 
appeared in the equations of motion for the four-degree-of-freedom model 
(eq. (26)). With this result, a simple design criterion can be derived for 
optimum proprotor/pylon dynamic stability. If the total factor of B \q and Big 
is set to zero, with a proper choice of the parameters, the moments due to the 
rotor flapping dynamics (at least the low-frequency response) will not be 
transmitted to the pylon. Then the pylon is effectively decoupled from the 
rotor and will not respond to rotor flapping motion: this is expected to 
improve the stability of the proprotor/pylon system. When the rotor lags the 
shaft in response to shaft angular velocity, that tip path plane tilt produces 
a hub force through 2Crp/oa + tfg*, which is dominated by the negative H force 
term hence it results in negative damping of the shaft tilt motion. This 
effect, in particular, is nullified by setting the total tip path plane tilt 
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factor of equation (105) to zero. The thrust vector tilt term 2 C T /aa 
contributes positive damping, but it is negligible compared to H g* for high 
inflow. Thus the criterion essentially is to use the positive damping avail- 
able from the hub moment of a hingeless rotor to counter the negative aerody- 
namic damping of H g*, and the design rule becomes a specification of the flap 
frequency vg required to satisfy that criterion. 


Setting the coefficient of 0^ and 0^ in equation (105) to zero gives the 
design choice for v^: 
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For typical values of y, h> and V for proprotor operation, equation (106) 
requires that vg - 1.1 or so such a frequency is easily obtained with current 
hingeless rotor technology or even with an offset hinge on an articulated 
rotor. An obvious limitation is that is a function of V while is not, 
so the design criterion can be met only at one design speed. 

This result was first obtained by Young and Lytwyn (ref. 18). On the 
basis of optimum flutter stability with the pylon decoupled from the rotor, 
they obtained the above result for the optimum flap frequency Vg. Numerical 
calculations of the flutter boundaries for a four-degree-of-freedom model (a 
flapping rotor on a pylon with pitch and yaw motion) show that the optimum 
flap stiffness is actually somewhat above this. They suggested that a value 
about 5 percent higher be used: 


v 2 = 1.05 
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With this choice of Vg, the remaining rotor force acting on the pylon is 
mainly a negative spring force due to the speed stability coefficient #g*; 
this force leads to a divergence criterion, which does not require very much 
pylon stiffness. They also point out that the only parameter of the pylon to 
enter this criterion on Vg is the mast height h, because it is essentially a 
static decoupling criterion, so the pylon dynamics do not enter at all. 

Unfortunately, the problem of proprotor dynamic stability is not so 
easily resolved. This criterion for vg is based on pylon pitch and yaw motion, 
while, for the rotor on the cantilever wing, the degrees of freedom produce 
much different shaft motion and respond differently to hub forces and moments. 
This criterion is based on only the low-frequency response of the rotor, and 
only the flap motion at that; the dynamics involved with high inflow proprotor 
instabilities is more complex. However, many other design considerations are 
involved in the choice of wing stiffnesses and blade frequencies; therefore, 


75 


while the proprotor dynamic stability must be verified, and perhaps the 
parameters modified somewhat to increase the margins, the actual choice of the 
rotor and wing configuration (including the placement of the fundamental fre- 
quencies) is not necessarily based on the criterion for optimum dynamic sta- 
bility. The above result indicates that, for practical applications, a 
flapping frequency greater than 1/rev should be favorable for rotor and wing 
dynamic stability. 


Whirl Flutter 


This chapter examines the classical whirl flutter problem: a truly rigid 

propeller on a pylon with pitch and yaw degrees of freedom. With the infinite 
blade stiffness of a truly rigid propeller, there is no blade flapping motion. 
The theory developed here may accommodate the rigid propeller by letting the 
flap natural frequency go to infinity; in that limit, the four-degree-of- 
freedom model reduces to two degrees of freedom, pylon pitch a and yaw a x , 
which is the usual whirl flutter formulation. Nomenclature particular to the 
whirl flutter analysis is given at the end of this section. 

The whirl flutter model may exhibit instabilities at high forward speed 
or low pylon stiffness, primarily because of the high inflow aerodynamics. In 
fact, it has much different behavior than the flapping rotor case. The whirl 
llutter problem is considered here because it is a special limit of proprotor 
dynamics, often mentioned in the literature as an ancestor of the current 
proprotor analyses; and because, as a two-degree-of-freedom problem, exact 
expressions for the stability boundaries may be obtained. However, the 
dynamics and aerodynamics of the whirl flutter model are found to bear little 
relation to those of the flapping proprotor. 

The whirl flutter equations can be derived from the four-degree-of- freedom 
equations of motion already given (eqs. ( 26 )), with the limits vg -> <*> and 

0. However, (vg 2 - l)Bi£ and (vg 2 - 1 J B 1^7 must remain finite in this 
limit because the propeller can still transmit a hub moment to the pylon; then 
the $^0 ar *d 6 ec l ua tions become simple substitution relations for the pylon 
equations. A more direct approach is to first substitute for the hub moment 
(vg 2 - 1)S 1( 7 and (vg 2 - 1)8^ in the pylon equations of motion, and then take 
the limit of Bi^uBjg **• 0 . The pylon equations of motion for the four-degree- 
of freedom model (eq. (8)) are 
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which is simply a mass/spring/damper system for pylon pitch and yaw, forced by 
the rotor moments and inplane forces applied at the hub. The equations have 
been made dimensionless (using p, ft, and R) and the inertias, damping coeffi- 
cients, and spring constants have been normalized by dividing by The 

hub moment is 


Y 


The flapping equations of motion from the four-degree-of-freedom model (eq. (8) 
for N > 3) are rearranged as 



Cv 
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+ 




( 110 ) 


With this result, the hub moment in the pylon equations of motion may be 
replaced; then the indeterminant combination (v^ 2 - 1)8 no longer appears, and 
the limit 6ic,3lS -> 0 may be taken properly. In other words, the and 
degrees of freedom may be dropped after substituting for the hub moment, thus 
the equations of motion are 



The pylon is now forced directly by the flap moments on the disk ^F^q and 

which are now transmitted directly to the hub rather than through the flapping 
response. 
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The aerodynamic forces follow directly from equations (23) for the four- 
degree-of-freedom mode, with the S 1(7 and B 15 terms dropped, of course. The 
hub moments and hub forces are then (hub moment is obtained by use of eqs 
(109) and (110) with eqs. (23)): 




(113) 



A comparison of these expressions with the results obtained for the 
flapping proprotor dynamics shows that the whirl flutter problem involves con- 
siderably different behavior. Comparing the hub forces with equations (101) 
indicates that the response to shaft tilt, cyclic control, and gusts is the 
same (except that eqs. (101) are in body axes, so the a x and a y terms are 
dropped). With the rigid propeller, however, there are no longer hub forces 
caused by flapping, acting through the thrust vector tilt and the negative H 
force due to inflow (the 2 Cy/oa + factor). MCtb forces due to the tip path 
plane tilt are naturally not part of the whirl flutter problem since they arise 
only with the flapping rotor, and that difference is very significant. 

Comparing the hub moment with equation (89) (the hub moment for the 
flapping rotor if it is multiplied by (vg 2 - 1 )/y) shows that the role of the 
flap moments is altered. Hie rigid propeller limit of vg 2 -> °° corresponds to 
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the limit N* -*• » in equation (89), and with that the two results are identical 

(except that in eq. (89) the a x and o„ terms were dropped again for body axes, 

and the inertia terms and d y were 'dropped for the low-frequency approxima- 
tion). The more appropriate limit for studying the flapping rotor dynamics is 
vg 2 -+ 1, that is, N * -* 0. Between N* -+ “ and N* -*■ 0 there is a 90 azimuthal 
phase shift in the influence of all flap moments. For the rigid propeller, 
the flap moments on the disk are transmitted directly to the hub while, for 

the flapping rotor, there is a 90° azimuth lag (for vg = 1) from the applica- 

tion of the flap moment to the achievement of maximum tip path plane tilt. 


Hence the role of the rotor in the whirl flutter behavior changes greatly 
from that for the flapping rotor. The Coriolis flap moment and the speed sta- 
bility moments on the disk (My) now are off-diagonal terms (i . e. , coupling 
terms"), while the flap damping disk moment (Mg) now contributes directly to 
the damping of the pylon motion. The hub inplane forces are given solely by 
the direct terms and H g, with no contributions from tip path plane tilt. 

The speed stability coefficients, My and have a dominant role in whirl 

flutter dynamics. 


If one substitutes for the hub forces and moments, the equations of motion 
for the pitch and yaw motions of a pylon with a rigid propeller are 
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These equations are valid for propellers with three or more blades. As usual 
the notation of the aerodynamic coefficients indicates their source: H, f° r 

hub forces; M, for flap moments; subscript y, for inplane velocities; and sub 
script 8, for flapwise velocities. 

The inertia terms are the sum of the pylon and rotor inertias: 


I I + {N/2M b 

V + 1 = + 1 = V {R/2)I b 
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where I y + (N/2)I b is the total moment of inertia of the pvlon and rotor (for 
N > 3); I already includes the contribution h 2 NM b of the rotor mass to the 
moment of inertia about the pylon. The damping terms (diagonal) have contri- 
butions from the structural or mechanical damping of the pylon (C * or C x * 
and from the aerodynamic forces on the rotor. All contributions to the 
diagonal damping terms are positive (H u and -Mg are positive) so the net damp- 
ing is always positive. The aerodynamic contributions result from the flap 
damping of the rotation of the hub plane by a x or d„, and from the inplane 
force due to the hub inplane . velocity during pylon angular velocity. Based on 
the equivalent radius approximations, the aerodynamic damping terms are 


Y [~M- + h 2 {H + /?)] = y(S°*± + h 2 V_sin »\ 
p \i V J \ 8 2 ) 

The first term decreases some with V, while the second term increases. The 

spring terms (diagonal) have contributions from the structural restraint of 
the pylon deflection (Ky * or K x *) and from the aerodynamic forces on the rotor 
rhe aerodynamic term, a negative spring (// p > 0) from the inplane force due to' 
the inplane component of ( V + y) produced by pylon pitch or yaw, is 
approximately 


-hy(V + v) (H + R ) 

a p' 


= -hy 


V' 1 sin $ 
2 


which increases with V. The coupling spring term (off-diagonal) is entirely 
aerodynamic, from the flap moment on the disk due to the inplane component of 
V + v produced by pylon pitqh or yaw; it is approximately 


which increases with V. The coupling damping terms (off-diagonal) have inertia 
and aerodynamic contributions. The inertia term (-2A X or 2<4 y ) is the Coriolis 
acceleration produced by the pitch or yaw angular velocity of a rotating body. 
The aerodynamic terms are from the flap moment due to the inplane velocity 
produced by the pylon angular velocity, and from the inplane force due to the 
angular velocity of the hub plane. From equation (59), it follows that 

3 C 

hy(M + #•) = hy — - « 2 
M p oa 

hence the aerodynamic contributions may be neglected compared with the 
gyroscopic coupling. 

The whirl flutter model then consists of a mass/spring/damper system for 
pitch and yaw, with positive aerodynamic damping, a negative aerodynamic 

spring, gyroscopic coupling in the damping, and aerodynamic coupling in the 
springs. ' r 6 
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Consider the single-degree-of-freedom problem: K x * -> « so a x -*■ 0; the 

equation reduces to 

(. Iy * * 1 )Sy * [Cy* + h 2 Y + i? y ) - Y^ ] ky * [*/ - h? ( 7 + tf ) (fly + 

= hyH B Q ls + + V 0 ^ (115) 


The damping is always positive, so the system is stable in this limit so long 
as the net spring constant (structural plus aerodynamic) is positive: 


K y * > hy(V + v)(H v + fl y ) 


(116) 


This is simply a static stability boundary, that is, the divergence criterion 
(for K x * -> “) . With one spring stiff enough, the system always has dynamic 
stability; the only stability criterion remaining is the divergence boundary 
for the other spring, due to the negative aerodynamic spring term. 

Consider the case of no aerodynamics, but with gyroscopic (Coriolis) 
coupling included. For the isotropic case (same mass/spring/damping constants 
for both pitch and yaw) , the characteristic equation is 

[ (J* + 1)X 2 + C*A + K*] 2 + (2A) 2 = 0 (117) 

The four roots will have negative real parts so long as y, K* s and I* are all 
positive. Therefore, the system is always stable, even with gyroscopic 
coupling, if there are no aerodynamic forces on the rotor. 

It follows then that any dynamic whirl flutter instability can only be 
due to the aerodynamic cross-coupling of the pitch and yaw motion; there is 
only one such term, namely, the off-diagonal spring terms yVM^. Since this 
aerodynamic force increases with forward speed V, an instability is expected 
to occur eventually as the inflow is increased for a given pylon and rotor. 

Of course, there is the possibility of a static instability (divergence) due 
to the negative aerodynamic spring term. Whirl flutter instability is then 
the result of the high inflow rotor aerodynamics. 


In Laplace form, the homogeneous equations of motion are 
(Iy *+ 1 ) s 2 + [Cy *+h 2 ' y (H „ +R ]1 ) - yWg ] s 
+ [Ky*-hy(V+v) (W ] 


-2s+y(V+v)M Vi 


2s-y (V+v)M u 


(I x *+l)s 2 + [C x *+h 2 y (H v +f? y ) -y Mg ] s 
+ [K x *-hy(V+v) (fl u +/? y ) ] 



= 0 


(118) 
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To simplify the notation, write V for V + y, // p for l! v + R v , and I* for T* + 1 


I*s 2 +(C*+C u )s+K*-K u 


s ^\i 


-s£> p +L p 





= 0 


where 


C u = - Y^e 

D v = 2 
L M = YWtfp 


(119) 


These are, respectively, the aerodynamic damping and negative spring, the 
gyroscopic coupling (which is given a symbol so that its influence can be 
identified in the results), and the aerodynamic cross spring. Combining the 
structural and aerodynamic terms in the spring and damping coefficients yields 

IyS 2 + Cy8 + K -sD + L 

sD - L I x s 2 + CpS + K, 



This last and simplest expression is used for the derivations. 


The characteristic equation for the whirl flutter eigenvalues X is then: 


(TyX 2 + Cy\ + Ky)(I x \ 2 + (T^ A + Kg) + (~\D + L ) 2 = 0 (121) 


The stability of the system may be investigated by examination of the root 
loci, that is, the behavior of the eigenvalues with variations in V or some 
other parameters. For a two-degree-of-freedom problem, however, it is possi- 
ble to solve explicitly for the stability boundaries. The stability boundaries 
on the K x * , K y * plane (i.e., as a function of the pylon stiffnesses) is the 
usual form for the presentation of the whirl flutter solution. The other major 
parameter is the inflow ratio V. Two kinds of instabilities are possible with 
a system that can be described by linear constant coefficient differential 
equations. The first is divergence, a static instability that, on the X plane 
(root locus), takes the form of one root on the real axis that passes through 
the origin into the right half plane. The second type of instability is 
flutter, a dynamic instability that takes the form of a complex conjugate pair 
of roots that cross the imaginary axis (at finite frequency) into the right 
half plane. The stability boundaries for these two types of motion are 
examined. 
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Divergence- The divergence boundary is defined by the requirement that 
one of the roots of the character Lstic equation be at the origin. For X = 0, 
the characteristic equation gives 


Vx * Ll - 0 

or 

<V - V cv - V + V = 0 

therefore, 

(K y * - hyVH v ){K x * - hyVH v ) + (y® y ) 2 = 0 

This sum must be greater than zero for divergence stability. An instability 
may be encountered if one of the spring constants, Ky* or K x * , is too small, 
so the negative aerodynamic spring r,aake:S the first term and hence, perhaps, 
the sum negative. Sufficiently large structural spring restraint of pylon 
pitch and yaw motion guarantees divergence stability. 

The equation for the divergence boundary is a hyperbola on the Ky* , K x * 
plane; the divergence hyperbola is shown vn figure 6. The asymptotes are the 
K x , Ky axes, that is, Kx* or K y * = K v = hyVH . The asymptotes give the diver- 
gence criterion when one of the springs is very stiff; for K x * 00 , equation 
(124) reduces to (for stability) 

Ky > 0 (125) 


( 122 ) 

(123) 

(124) 


or 


Ky* > K v = hyVHy (126) 

which, indeed, is the criterion obtained previously for the single-degree-of- 
freedom limit. This criterion is that the total spring constant - aerodynamic 
plus structural - be positive. For finite K x * , the cross spring L y = yWy 
introduces coupling of the and oq, motion, which eases the criterion on Ky*. 
The hyperbola intersects the line K^. = -Ky (the? -45° line on the K Xi Ky plane) 
at points Ky = = fyFiWy, where the closest stpproach of the two hyperbola 

branches occurs. The minimum width of the divergence-stable corridor is thus 
2/2 = 2j2y VM^. These points occur in the fir:;t quadrant of the K x *,Ky* 
plane (fig. 6), only if Ly/Afy < 1 (which is., in fact, true only for rather 
large 7). In any case, the cross spring L y is a good measure of the corridor 
width at low K x * and Ky * . 

In terms of the approximation based on the equivalent radius, the 
negative aerodynamic spring - which determines the asymptotes of the divergence 
boundaries - is 


K u = hyVH u = hy 


V* ' - sin 
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which increases with both V and h, so both decrease divergence stability. The 
minimum width of the divergence-free corridor at low K* is 


2/2L u = 2/2yVM u = 2 /2y 


which increases with V. 

liie divergence boundary is then determined by two aerodynamic influences: 

the an 'ct negative spring K y[ = hyVH^, which gives the location of the asymp- 

totes, and the cross spring coupling - yVM which gives the width of the 
corridor at low K* independent of A' p . The first is the direct hub force due 
to the rotor and the second is the flap moment on the rotor - both react to 
the inplane hub velocity component of V •+ y due to the shaft tilt by pylon 
pitch or yaw. Both spring terms increase with inflow ratio V; the effect of 

the direct negative spring also increases, with h since it is a hub force. 

The pylon and rotor will always be statically stable if both K* and K * 
ai e greater than the negative aerodynami e spring AT, which is usually not a^ 
very stringent criterion on the pylon stiffness. Static stability is also 
achieved always for the nearly isotropic case, K x * = K u * , even if the stiff- 
nesses are smaller than A' y . The cross spring L u = y7< stabilizes the' diver- 
gence motion if the springs are not. too stiff, so that both degrees of freedom 
are active and the coupling of the aerodynamic spring may be effective. This 
effect produces the divergence- free corridor at low A'*. 'however, the low K u * 
and K x region will be within the flutter instability region, so the presence 
of the divergence-free corridor has little practi cal' application. 

Fl utter: Some preliminary c o nsiderations- Some general results for whirl 

flutter dynamic instability are presented here; they were originally obtained 
by Young and Lytwyn (ref. 18). The equations of motion are 



Forming the vector product wi'th [eg d j: ] (i.e., the sum of a u times the first 
equation and a x times the second)," one obtains 


ip [2 W + I w + 1 K i' a y 2 4 I W] 


[Cydy 2 + C x d x z + A (1 (a x a.y - a^d^) ] (128) 


Fhis equation is an energy balance for the whirl flutter motion. The left- 
hand side is the time rate of change of the sum of the kinetic and potential 
energies of the system, which must always be positive (assuming that the diver- 
gence criterion is satisfied). If the right-hand side is negative, the total 
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energy is decreasing with time, hence the system is stable; if it is positive, 
the total energy is increasing with time, so the system is unstable. The damp- 
ing contributions are always stabilizing and extracting energy from the system. 
An instability is possible only if L y (a x d y - a y d x ) is sufficiently large and 
negative. Again, whirl flutter is a direct result of the aerodynamic spring 
coupling L v = yVW p = yF(sin <t>)/6. Since this coupling increases with F, an 
increase in the forward speed eventually produces an instability (if the term 
is negative) . The gyro coupling due to the rotor (C*y) does not appear in this 
result, but it does influence the mode of motion of a x and o.y and hence may 
influence the stability. Without L y , no instability is possible no matter 
what influence has on the motion. 

Since the aerodynamic coefficient L y is positive, a necessary condition 
for an instability to occur is that a^d,. - a l ,a x < 9 . Writing a.y = q sin 6 
and cep = q cos 0, it follows that a. x a.y - &yO- x = < f®* Hence the requirement is 
that 6 < 0, which means a mode of motion in which the shaft whirls in the 
direction opposite the rotor rotation. This mode is called a backward whirl 
mode; the motion with the shaft whirling the same direction as the rotor is 
called forward whirl. 

These considerations show that a whirl flutter dynamic instability for a 
rigid propeller on a pylon occurs only in a backward whirl mode, and that it 
is a high inflow instability due to the aerodynamic spring coupling 

Flutter- The flutter boundary is defined by the requirement that a complex 
conjugate pair of the roots of the characteristic equation be on the imaginary 
axis, that is, that roots with zero real part be a solution. For X = iu, 
where w is real (u> 2 > 0), the characteristic equation becomes 

(iuCy + K y - u> 2 I y )(iuiC x + K x - u 2 I x ) + i-iuD + L) 2 = 0 (129) 


The real and imaginary parts of this equation are 

(. K y - m 2 Jy) {K x - a > 2 I*) - u 2 C x C y + L 2 - D 2 m 2 = 0 ’ 

- b^Iy^Cx + ( K x - t^ 2 1 x) Cy - 2DL] - 0 

If it is assumed that w ^ 0, the imaginary part may be solved for uT 

2 tyx + _ K x C y ~ 2DL 


Substituting the solution for co 2 from the imaginary part (eqs. (130)) 
into the real part produces 
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K 2 I 2 C C + K 2 I 2 C C - 21 I C C K K + K [C (I C + I C ) (C C + D 2 ) 
y x x y x y x y x y y x x y y a y x x y y x 


+ 2DLI (Cl - C I )] + K [C (I C + I C )(C C + .D 2 ) - 2 DLI (Cl - C I )1 

X X Ij It X X If u X XU U X U X U If x J J 


~[L 2 (I C + I C ) 2 + 2DL(I C + I C )(C C + D 2 ) + T I 4D 2 L 2 ] = 0 (132) 

y x x y J y x x y y x x y J v } 


Tliis general second-order equation for K x and Ky defines the flutter stability 
boundary; it is, in fact, the equation for a parabola. The equation is 
simplified by a transformation to a new axis system .Y and Y defined by 


T y 

x = t K y + — K z 


z x 


I, 


7 = T h - T X* 


(133) 


where I 2 = J r 2 + J,, 2 . This transformation is a rotation by the angle 
tan~* I U J I x from the Ky and K x axes; for isotropic inertias, ly = the 

rotation angle is 45°. r rhe X and Y axes remain orthonormal uncler this rota- 
tion. The relation of the X and Y axes to the original K x and Ky axes is 
shown in figure 7; the K x and Ky axes arc, of course, already shifted by K y 
(the negative aerodynamic spring) from the structural spring axes, K x * and K u * , 
With this rotation of the axes, the equation for the flutter boundary becomes 


+ 7 


T I (C 2 - C 2 ) + (I 2 - I 2 ')C C 
x y x // v x y_ x y 


r-c c 


C + K/ 21 'L (i c - i c ) 

C u u x X u 


OCT 

X u 


+ IC) : 


+ u 


I Z C C 


I 


I C + I C 


T o U X 


a if 

+ 


rl c c 


)L(r c + i c ) c c 

J T* T' i y r Y' * I 


+ D 


n 2 


y x x y x y 

~fr~c k 

x y 


I 


T x I u 4D 2 L 2 


T±C C 
1 l 'x L/ y u 


= 0 (134) 


which is, indeed, the equation for a parabola. 


Tli is equation may be simplified further by expressing it in a form based 
on the isotropic case. Let zj and zq be measures of the anisotropy of the 
moments of inertia and damping coefficients that 


tan Zj. 
tan e 






i,, + i. 


■'x 


C, 


r + c 

■ x y 


( 135 ) 
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and, as representative inertia and damping, use 

J = Jl x 2 + 1/ ' 


(136) 


C - + 

The angle er is the rotation of the X axis from the 45° line on the K x ,Xy 
plane (fig. 7). With these definitions, the equation for the boundary becomes 


Y 2 + Y 


+ X 


I cos 2e 
.2 


sin 2(e + £„) ( 

1 C (t 2£ c * 02 


C 

Z C ] (c 2 


) 


ADL 

C 


sin (e^, - e r ) 


J 


cos 2e 


C 


2 cos 2 (ej + e^,) / 
I cos 2e p \ 


Y COS 2e p + D 2 


) 


2 L' 


COS 2 (£j + Zq) 


cos 2e 


C 


ADL 

IC 


COS (£j * E ) / C 2 

cos 2 e„ 

cos 2e^, \ 2 c 


4D 2 L 2 
+ D 2 j + — 

/ C 2 


cos 2e 


I 


COS ZE 


C 

(137) 


The equation is now in the standard form for a parabola, namely, 

(7 + ey) 2 + 4dZ = Adh v (138) 

The geometry of the parabola is described by the shape factor d and the vertex 
h (fig. 7). The axis of the parabola is parallel to the X axis, but shifted 
a distance Ey below it. The maximum value of X on the boundary occurs at the 
vertex, a distance h y from the Y axis. The parameter d describes the shape of 
the parabola: at a distance d from the vertex, the parabola has a width of Ad 
Moreover, the isotropic inertias (so Ey = 0 and the parabola axis is parallel 
to the 45° line on the K y * } K x * plane), the shape factor points (points A m 
fig. 7) give the maximum values of K y * and K x * on the parabola. 

To first order in ej and e c , that is, for small anisotropy, the parabola 
axis shift, vertex, and shape factor are given by. 


ey = 


C 2 + 2D 2 

21 


(ej + e c ) 


2DL 


( G , 


Ej ) 


(139) 


, „ C 2 + 2D 2 
d = 4J 

For small anisotropy, there is no first-order influence on h v and d, which 
describe the parabola shape and position; there is only a small effect on the 
orientation of the parabola axis, which is rotated by Zj and translated by £y 


(140) 

(141) 
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from the 45° line on the K y ,K x plane, 
the equivalent radius, the vertex is 


In terms of the approximation based on 


where 


( i k + 

eye 


■) 


/ 

(7* 

/2 (C* + C u ) \ 


C* + C + 


) 


= /?- 


Y VM 


(n 


Y VM 


C + h yHy - y M q \ (17*+ hp-yft-y - yM'^ 


* 2 ) 


(142) 


Y VM 


C* + C, 


Y^fsin <t>)/6 


C* + h 2 yH, 




C* 


+ Y [ (cos <|>)/8 + h 2 (V sin <p)/2] 


(143) 


The farthest penetration of the parabola that defines the dynamic stability 
boundary , s governed by the vertex V The size of ^ is determine? pj^arily 
i' / ; hp h ! la ° g t J e coupl mg aerodynamic spring to the total pylon damp-' 

decreases b^^The 0 pylon ’ Wlth structural and aerodynamic contributions, 
decreases h v The aerodynamic damping probably does not change much with V, 

^ mast hei^ht^ f? h while the ° ther “creases. An increase 

damoinv H m • he i pful for stability since it increases the aerodynamic 
.im P ing The major influence of V is on the coupling spring 7 U = yVM u = 

Y^sin J)/6, which increases with V , and so h v does also. The gyroscopic cou- 

sS/l t'S K he inStability ’ but the Primary source is the cross 

dn^H g K P +w H ? ir hub m ° ment due t0 the in P la ue component of V + v pro- 
duced by the pylon tilt, 7, = Y ^). The shape factor is approximated 


d = 


C 2 + 2 7 2 
4T 


(C* * ♦ D„2 

2/2 I * 


(c* + h z yH v - Y A/g) 2 + 4 
2/2 7 * 


= K* * yT (cos »)/8 + 7 2 (T sin j>)/2] } 2 + 4 


2/2 7 * 


(144) 


Smail d means a flatter parabola; therefore, in this respect, damping 

77 E the inertia J* is good, in contrast to their 


hon 1Ue ?h 6 ° nhv \ - The shape fac tor is relatively independent of V, however- 
J™ , the mam influence of an increase in the inflow ratio is to increased 
hlch shlfts the Parabola outward without changing its shape much 


v> 
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It may be concluded then that the whirl flutter dynamic instability is 

due to the aerodynamic spring coupling - yVM v , as suggeste y ea 

aue to tne ue y f * parabola is relatively independent of V, 

while°the v“l increases greyly »lth P; forward speed therefore decreases 
the flutter stability by shifting the parabola outward. Since an increase 
damping either structural or aerodynamic, makes the parabola broader, u 
decreases 6 the vertex size, the overall effect on stability is favorable Con- 
versely, an increase in the inertia is unfavorable is favor- 

enters into the aerodynamic damping of the pylon, where i cannot 

able. The gyroscopic coupling participates m the mstcb 
cause it if is small enough. 

Construction of boundaries- From the previous solutions, the whirl flutter 
stabilSfra^ ies can faT Sftly sketched on the K x *K y * plane for a given 7. 
The relevant axes for both flutter and divergence are the K x ,K y axes on tne 
* Plane, which are offset from the V.V (structural spring) axe by 
the negative aerodynamic spring K v - hyVH v . felines K x - K v and K y 
are constructed first . 


The divergence hyperbola has as asymptotes the ^ ant 
minimum width for small Kr* and K u * of 2/2 L u - 2/2 yVM y . 


minimum widtn tor smaii K x * and Ky* ... y - ■ m f nrqr tical 

of the divergence boundaries within the flutter boundary are not of practica 

Interest? The divergence boundaries are defined primarily by the asymptotes, 

° line from the K u and 


and a 
Since the portions 


and Ky axes, 


The flutter boundary construction begins at the 45 

a 1 .4-1 O "t* 


axis? ThTFald r axis are constructed, then rotated an angle e T from the 
45" line; the parabola axis is then constructed a dl f c , ic:low th ^ -• 

The vertex of the parabola is given by the distance y , 

the shape factor d, four more points on the parabola are easi y oun 
“axis and points A as in fig. 7). The shape factor points (Am fig. 7) 
give the points of maximum X y > and K x \ at least for isotropic inertias 

(EJ =0). 

It is almost as easy to construct the boundaries exactly as it is to 
sketch them. From the equations for the divergence hyperbola and fl u t 

parabola (eqs. (124) and (138)), a number of points on each boundary can 

obtained quickly. Typical results are shown in figure 8 for V/W - 

h - 03 Y = 4 It* = Iv* = 2, and c x = Cy* = 0 t hl 8 h inflow typical prop 

rotor mast height and lick number, isotropic inertia of the pylon, and no 

nvlon structural damping). Also shown for the K x * axis is the equivalent 
pylon structural a<*mpxug, * f n to aid in the interpretation 

pylon natural frequency, u x - K x /I x (.per revj , 

of the results. 


W hirl Flutter Nomenclature 

q _ ^2 direct aerodynamic damping 
v = hyVH^ direct aerodynamic spring 

D = D = 2 Coriolis coupling 

L = L y = yWp cross aerodynamic spring 
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Whirl Flutter Nomenclature (Concl.) 


K y ~ K y* ~ K \i 


%X ~ ~ & 


c y ~ c y* + ^ 


~ + C 1 , 


x 


u 

I„ - I 


x 


e T = tan " 1 


y 


y 


£ q - tan" 


^ ~ c y 
^x + ^ 


total pitch spring 
total yaw spring 
total pitch damping 
total yaw damping 

inertia anisotropy 
damping anisotropy 


4 2 * V 

" - -C 2 ” C 

X,Y 

e Y 

d 


representative inertia 
representat ive damping 
whirl flutter parabola axes 
whirl flutter parabola axis offset 
whirl flutter parabola shape factor 
whirl flutter parabola vertex 


Two-Bladed Rotor 


experiment s^fthe^XV-S^and a^sIr^Thlve^S ^th' ' 

differential equations. 7 described by periodic coefficient 

of freeSom^^e'rot'o^Ha flappin S. rotor on a pylon with pitch and yaw degrees 

coning mode. ^ TolVt iLFT** ° f * teeteriag m ° de and a 

moment to the pylon to excite its n^cb ° 0t tranSmit net force or 

case of a two-Maded rotorls the Veit m ° tl0n : Moreover > the usual 

levered to the hub which is then a rt W rot or, which has both blades canti- 

For the teetering mode JheroJor acts" t0 the Shaft by 3 sin S le fla P Mnge. 
coning mode it acts a^ l hi i 30 artlcula ted rotor, while, for the 

Thus the coning mo ™ y \e "fg S ^TtS^T' h ‘*V““» I fluency, 
to three degrees of freedom: lotorleMer^g and '7 reduCes 

differs from that for three or more blades, where L ^W^e" mo^i ha?” 
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four degrees of freedom.) For N > 3, the rotor motion in the fixed frame is 
described b) tip path plane tilt degrees of freedom (plus N-2 other modes that 
do not couple with ax and «„) . For ill = 2, however, such a description over- 
describes the transient motion since only one degree of freedom is mvo \cc 
the teetering mode. For N >- 3, the use of tip path plane coordinates results 
in constant coefficient differential equations, even though the equations o 
motion for the individual blades in the rotating frame involve a Periodic 
variation of the inertia and aerodynamic forces. For il 2, the 1 _ . 

HcJenJs remain in the equations of motion since the teetering mode is really 
a rotating degree of freedom. The fundamental difference between the dynamics 
of the N > 3 and N = 2 cases, which makes the two-bladed rotor much more c 
ficult to'analyze, is that fact: for N = 2, the equations of mot:on for the 

coupled rotor and pylon system involve periodic coefficients even m pure > 
axial equilibrium flow. 

The equations of motion for the pylon and rotating blades are the sameas 
for the four-degree-of-freedom development (eqs. (4)). If the coning mode is 
neglected since it does not influence the coupled dynamics anyway, the flap 
ping angle of the mth blade (m = 1,2 here) is given by the teetering degree 
of freedom 6: g( 2 ) = B for the blade at azimuth position i|>, and g( J " /J 

for the blade at azimuth position + it . The equation of motion 
tering mode is obtained by operating on the rotating equation with 
1/2 y (-l) m (. . .); this operation yields the total teetering moment on the 

rotor. With sin \p m = (-1) sin i p and cos ip m 
is 


(-l) m cos ij), the equation for 


M r 


3 + v 2 6 - (cty - 201^)005 + Ca x + 2ct y )sin ip - Y 


a.o 


(145) 


where Mp = 1/2 ^ (-1)^? is the aerodynamic teetering moment. Already, 

neriodic 1 coefficients appear with the pylon inertia terms as seen by the rotat- 
ing blade. The evaluation of the hub moment and force for the pylon equations 
and the flap moment for the rotor equation follows exactly that for the ou - 
degree-of- freedom model, up to the point where the summation over the N K lade 
is performed. At this point, the influence of N = 2 appears, usually the 
introduction of periodic coefficients. 


The following 
written for V + v , 
for H v + R u ) : 




equations of motion are obtained (for convenience, V is 
except in the gust terms, where it really is just V); and 
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’ YA?e (^V^)si„ HWb cos * (2 + y^ p )cos (,-YAf- Sin 0 ~| /6 

-;^ 2 Sin * Sin 2, , V(J sin 2 ,.^- (1 . cos w / 

-nyli 6 l cos 4 - h 2 yH v sin 20 + ; IY ffg (l + cos 20) <V+*V/p (l*cos 20)-^ sin 20 |\/ 


v +y Vp 


-YW p sin 


( V 3 * 1+y/? aa) 2 cos ^ +;z y // 0 A >- si ^ * Ky*-hyVH v (1 -cos 2*) 

-(v B 2 -l + y« -~J 2 sin ^h y H e K p 2 cos 0 -AyW sin 20 


?p sin 2\l> 


-yVM cos ^ 


"^Y^p sin 


Kx*~hyVH^ (1+cos 2^) 


^y^02 cos ^ 


Y*W sin ip 


yVM v cos 2^ 


= hyH 6 2 sin 0 01 + ftyW^l-cos 20) * Y W si 


p sin 2^ 


sin 2 * hyVH (1+cos 2^)1 X Gi 


. (146) 

Inis three-degree-of-freedom set of differential pnnatinnc u . , 

coefficients, with a period of 2ir tv.cc d e 1 uatlons has periodic 

techniques of Floquet theory for the e' 10ns can be solved, by the 

cate the stabilit? (as dlSe" T g . "Lf ^ 

involve, the XV-3, whirt L^letek "^r T 'TntTV ( T ' ”* 

^"•frss th q r?^:ieirr r; he h sol r d 

a more general form and even to delude Une1r h aerod°yJa e mics!' e eqUati ° nS in 

corresponding^ t^the t^ej^g [e ScpJ/g* "g™ 

WC ST svrrr-;!-: 1 jsrjatss 


— ~ cuu v^ii L±uua i swasr 

lly 0 1 C and e i 5 > where 01 = e cos </, + 


sin ip, 


.ii.i2; v re ob,ained " uch *? *»• * i 3. 

the flapping equation then dron the* a ° n e ^ U ^ tlons °f m °tion by use of 

a truly rigid propeller The result for^ 66 ° freedom ln taking the limit to 
v 1 me result, for homogeneous equations, is 
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I. 

_ 


i"^*+l+cos 2\l> 


-sin 2\ If 


I x * + 1-cos 2^ 


E v *^ 2 Y// u (l-cos 2^) -yAfg (1+cos 2*)- (2+2^) sin 2 * 
^ 2 Y J/ p sin 20+Y^g sin 2^+2(l-cos 2*)-2feyAf u cos 2* 


/-j 2 y gjj-j 2^+Y-^g sin 2i^-2(l+cos -2hytt^ cos 
C x *+h z Y^ p (l + cos 2110-^^3 C 1 -c° s 2 iii') + C2+2«Y^p)sin 2^ 



Xy^vV?/. fl-cos 2iiO+Y^ M sin 2^ 
-Jiyl-Y/^ sin 20- Y^p ( 1- cos 20) 


s i n 20 + Y^ , ' M (l + cos 2v ) 
AV-WW./1 +cos 20)-Y^: 1 sin 



(147) 


This is a two-degree-of-freedom problem again, but now with periodic 
coefficients (./period ,) because of the two-bladed roto^the solution for 
the stability is obtained by Floquet theory techniques. Tyj 
the stability boundaries on the K x and K„ plane are shown in figure 9, for 

v/i» = i, h . 0.3, v = 4 , i x »i- i x /mm - y 

divergence and flutter boundaries are much as for the W - 3 whirl 
problem (fig. 8). The periodicity of the system also introduces the possi 1 
ity of a divergence-like instability that occurs with a frequency 1/rev in 
addition to occurring on the real axis. (Such an instability may occur at any 
multiple of one-half the fundamental frequency of the system, which is 2/rev 
for the present case with a period of it.) Such instability regiopp 
the solution presented in figure 9. One region occurs for low stiffness, 
where the natural frequencies of the pylon (u ? and m y ) are near zero, that 
region is buried in the classical flutter region though. A 1/rev /eigenvalue 
also appears where u x or is near 1/rev, since that is w icre - 1 

would be expected to be near 1/rev. This new kind of mstabili y g 

proprotor and pylon system, introduced by the periodic coefficient., 
how special the two-bladed rotor is . 


AIRCRAFT STABILITY DERIVATIVES 


This chapter examines the contributions of the proprotor to the aircraft 
stability derivatives, and hence the influence of the proprotor on the entire 
airrraft The aircraft rigid-body motions are assumed to occur at low tre- 
q^y and^re “SSSttrf to the rotor shaft without modification by the wing 
motion. The expressions obtained previously for the response of the flappi g 
rotor to low-frequency shaft motion are used to estimate the rotor contnb 
tions to the forces and moments on the aircraft. This means, of course, that 
the influence of the rotor lag and the wing motion is being neglected. The 
rotor contributions are compared with the usual airplane contributions 
stability derivatives, from the tail and wing in particular, to evaluate the 
relative importance of the rotor terms. Nomenclature particular to the 
craft stability derivative analysis is given at the end of the chapter. 

While the assumption that the motion occurs at low frequencies for the 
aircraft rigid-body degrees of freedom is valid, it is not quite correct to 
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r ar{, „ 

“^Lf ™ ;'„r HtatiVC aSSCSS "' ent ° f ih * ”*”' ^ey Z very "° 

t,,o ^zrjzzi rrsrrr^ x f £r i ° ; , 

- &s {%±> €FH 

moments and the ri s kl-|,odv motions of the aft, „ ft „ITt 't ™ , forccs ;,nd 

1?Z£ f Sl n Foilow? STSTw* r ° t0r '-—-'and shaft «i«” 
center of gkJtV^T «S" “ “«»* “ *• 
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F, 


N.y* 


/V, 


7 / 


2 ( _< ? - -j /7 np-n ^ 

2(A/„ + hU) 


2 h 


c + hr - 




T sgn Yj 


(148) 


in figure S^ssumos^speci fi^ r two ? ont ™ rotatin S rotors. The rotor mode] 
the rear) for' t o o l potion of rotat ion (clockwise viewed from 

i- . 1 U1 l ,K »otot placed on the right wine as in fl anro i n tv 
direction of rotation is accounted for hv n f Y Y? 1°' ™ opposite 

3?JW$ ?jt r; SfSss 0 vr he 

"™° tho 

gravity is then rigid-body degrees of freedom of the center of 


Xt 


= - 7a - 


i. 


V 


p sgn Si 


i/p 


*P 


7 6 


£y 

- — r sgn Q 


a x = 


-r 


a y = <? 
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axes . 

_ the* aircraft derivatives to which there are 
Now examine some o . rnnsider the pitch moment due to 

contributions from the proprotors First, consider the pitc 

angle of attack, C m ; from its definition, 

a M,, + hH 


AC, 


'm„ 


Wy/4 S W d W _ _J_ 

. - Xr 


3a 

2ttckz 
^w a w L 
aa 


y(S.,/2A)c 


IC^/oa + fr(2C#/ca)”| 

*P J 

' 2 C y/aa + fr(2C g /ga)~ | 


W 


(150) 




rotor contribution in terms o h w - some tric parameters to account for 
involves the ratio of the r th | rot0 r and airplane coefficients. 

Z tl f 7o7^TCZ directly related^ tta.lrcr.ft ^tatic^in » 

posit ioiFof^the rtlck- fH^uiral poiM (as a fraction of 3„) . The change 
in the aircraft static margin due to the rotors is 


oa 


Ah ~ 

n V{S w /2A)c w C L 


2 Cm /oa + h(2C H /oa ) 

— — y~ 


'l 


Xp 


(151) 


the'rotors^ longitudinal stability. 

The vertical force due to angle of attack, C Za , is 

r 2 Cu/oa> 

W 4 H i H 


the static margin and h n is moved forward by 


A C, 


^a 


3 F z /qS. 
3a 


oa 


vs w X P V{S W /2A) 


^C H /oa\ 


(152) 


The contribution from 
C„ = -C L ; therefore 


the rest of the aircraft, primarily from the wing. 


is 


AC, 


oa 


C. 


V(S w /2A)C La 


^2 C H /oa^ 


(153) 
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A positive value means an increment in 
contribution. 


the same direction as the wing 


Tbe drag force due to the longitudinal veloci 


ty ’ C *u> is 


AC. 


x u ~ 


Zu 


4 T 


2aa 


VS, 


W p V( SJ2A) 


C T /aa 

k P 


2 aa 


V (S W /2A) 


( T \ ) ( 154 ) 


Because of the rest of the aircraft, C r = - 2 C n = a ■ * 

the rotors is x u D '‘-'w ari( l the influence of 


AC. 


x,. 


la a 


C. 


X , 


Vf/A 


(~ T \) 


(155) 


where a positive val 


^ a positive value means a contribution in tu 
aircraft drag (a damping term) Usinv the h 6 dlrectl0n as the 

lent radius for 7\ yields " 8 tH approxi matlon based on the equiva- 


A C. 


X, 


C '• 


X. 


- 2aa cos $ 
~ WTa 4 


contribution increases^S^ vvhile the aircraft drag 
with 7. The rotor thrust Lc to inHow ' T° T COntrib ^n decreases 
contribution is always in the same di£t& “ Srf^g^? ^ 


The pitch moment due to pitching velocity, C m is 


A C, 


■m - 


3N / qS c 

y ^ w w 

9 qr? (2V 
w 


2 a a 


V(S /2A)a 2 
w w 


2C^/aa + h(2Cjj/oa) 
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( 156 ) 


“ here « 

m - Th “ * he ‘™ « co°n”i ZZllTon ^ 


s+z 


tH 

~2A~ 


aa 

Va+ 


2C M /aa + h{2C }1 !aa) 


also be expressed in terms of the maneuve^ 15 y equ ; red - This coefficient may 
where M is the aircraft mass A pZTtTve AC^T 'r™*? ^ = , 

the stick-fixed maneuver margin, moves V fo SSJ? the . rotor then decreases? 

-re horizontal tail to maintain the sanS margin 


y 


( 157 ) 
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The rolling moment due to roll rate, C j, p » 1S 


oa 


AC C 




*P l J2V 2V{SJ2A) 


?2C H /aa 

X-n 


_8 V° a ' 

I 2 “2 

w 


(158) 


The contribution due to the wing is approximately - ' C L a /& ’ 

3 {CqJocl) / 3 d ,, — 


AC 


1 

IR = 


4aa 


2 V oa , ±_ Q . 


(159) 


'l-p V(S u /2A)C£ a 


; c. 


2 X 


zj 


The ^ term is due to the change in the rotor ^^*7/6, 

« th(; confribution irllwayfL the sLe direction ns the »in S ter», and 
increases with V. 


The yawing moment due to sideslip, 

F 


hCyic 


* N z /qS w l w 


V(S w f 2A)l. 


~2Cm^/ oa. - h{2C y/oa] ”| 

y p 


(160) 




The contribution of the vertical tail is C n $ = V?yV S iAj* so the tail 
required to counter the rotor contribution 


S v l v oa 


2A Va 


' v 


'2Cpi/oa - h{2Cy/oa) 


yi 
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where a positive value means that more tail is required. 
The yaw moment due to yawing rate, 


2oa 


ACV 


dN z /qS Jw M 

^ 72 ^ ‘ V(S w /2A)l/ 


- 2 C M /oa - h{2C Y /oa ) V 


2 s 7 
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The contribution of the vertical tail is C nr « -2 <^,v7W and, with 
3 (Cy/cra) /3^p = 


5 J, , 2 
v y 

24 


aa 


7a. 


2C Mx - h{2Cy/oa) ^ 

' d 2 X 

x 


(163) 
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The thrust coefficient ^ e is S negative re s e a et T r * qu j red to cou "ter the rotor. 

Trlrlfl the r0t ° r thrUSt d ^ in ^ contributes^o ^ ’ 


The rolling moment due to sideslip, c t is 


= 


3jV /qS Z 
x [ w w 

36 


oa 


2V(S u /2A) 


(sgn (2)1- 


2C y/ oa 

~Tp 


(164) 


«*«• r is the wing 

the rotors is r up) ' ™ en 'he equtvalent dihedral due to 


AT = 




V(S w /2A)C La 
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where a positive value means that more tail is retired. 
The yawing moment due to rolling rate, C„ p , ts 
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ROTOR FORCE AND MOMENT DERIVATIVES AT V/QB = 1, FROM THE LOW-FREQUENCY RESULTS 

(EQS. (102) AND (103)) 





cantilever rotor. These low-frequency derivatives, obtained with equations 
C102') and (103), were evaluated using the parameters in table 1; only vg, Kp, 
h y and V are required for the rotor derivatives, the others are aircraft 
parameters. The flap frequency is the primary parameter, responsible for the 
major differences between the derivatives for the gimballed and cantilever 
rotors; Kp is also important for the cross-derivatives of the gimballed rotor. 

The rotor forces and moments due to the hub inplane velocity xp or yp 
are basically produced by the coefficients // y and The derivative -H/xp, 

always positive and generally larger for the cantilever blade, gives positive 
damping then. The corresponding pivot moment derivative -(My + hH)/xp is also 
always positive. The cross-derivative -H/yp is always negative for the canti 
lever rotor, while it is proportional to Kp for the gimballed rotor and there- 
fore has its sign. The pivot moment derivative - (M u + hH)/yp is also 
proportional to Kp for the gimballed rotor and has its sign; while, for the 
cantilever rotor, the hub moment dominates the hub force acting on the arm ft, 
so the net pivot moment is always positive. The hub force due to the shaft 
angular velocity -H/a v is large and always negative for the gimballed rotor 
because of the flapping required to process the rotor acting through the nega- 
tive H force (Ho) to give a hub force; the derivative is small for the canti- 
lever rotor. The corresponding pivot moment -(My + hH)/a.y is then also negative 
for the gimballed rotor (negative damping); for fhe cantilever rotor, it is 
always positive because of the hub moment. The cross-derivative H/<x x is nega- 
tive for the cantilever rotor, and has the sign of Kp for the gimballed rotor. 

. * - n ir n n 1 1 Arl 


The moment derivative [My 

1 ■ ■ 4- • , 


ft//)/ also has the sign of Kp for the gimballed 


rotor and is positive for the cantilever rotor because of the hub moment. 

The longitudinal derivatives required (~T\, ~Q\> T rp and are a11 P osltlve 

and increases with V (all equal 0.135 at V = 1, based on just the terms as 

in eqs . (54) ) . 

On the basis of the expressions above and the values given for the lotor 
derivatives, the following influence of the proprotor on the aircraft dynamics 

_ h i — 1 -4- ^ -Jt-i /I c 


is found. The change in uuc uw i..i 

static longitudinal stability, that is 


due 


to the rotor always results in decreased 


forward movement of the neutral 


point h n . Typicaly, h n moves 


forward 10 to 20 percent of the wing mean aero- 


dynamic chord, with the larger value for the cantilever rotors. Like many of 
the rotor effects, this change is roughly proportional to H v /V, which, from 
the equivalent radius approximation, is proportional to sin 4>; hence the rotor 
effect increases with V. 


The rotor increases C ~ somewhat over the wing contribution, that is, 
the rotor contributes to the total lift curve slope of the aircraft by about 
10 to 30 percent (the higher value for cantilever rotors) of the wing contri- 
bution in the same direction, thereby increasing the total aircraft ' ' j • 

This effect is also proportional to H v /V, hence to sin $ (roughly) . lhe 
increase in the magnitude of C* (the coefficient is negative) produces some 
increase in the damping of the aircraft longitudinal short period mode. 


in the same direction as the aircraft drag 
contribution (making C Xu more negative). The rotor contribution is due to the 


The rotor contributes to C Xu 
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rotor thrust damping Ty and hence is an order of magnitude larger than the 
aircraft drag contribution. The result is an increase in the phugoid mode 
damping and also in the sensitivity to longitudinal gusts. The great increase 
in Cx u over the usual aircraft values also may change the motion involved in 
the longitudinal modes. 

The rotor contribution to C m q moves the maneuver point h m , typically 5 to 
7 percent of the MAC. The maneuver margin is decreased for the gimballed rotor 
because of the negative H force and is increased for the cantilever rotor 
because of the hub moment capability. This indicates the need for about 40 per- 
cent more and less horizontal tail effectiveness, respectively. Equivalently, 
because of the rotor, the short period mode damping and frequency are decreased 
for the gimballed rotor and increased for the cantilever rotor. Note that, with 
the high wing configuration of the tilt rotor aircraft, pitching about the cen- 
ter of gravity also introduces a T y contribution to C m , which always increases 
h m (less horizontal tail effectiveness required). ^ 

The rotor contribution to is of the same order and sign as the wing 
contribution. Both the rotor H force and torque terms have the same sign, and 
both (ffy/y and Q^/V) m increase with V approximately as sin <j). The H force con- 
tribution is about the same as that of the torque for the gimballed rotor, but 
larger for the cantilever rotor. As a result, the aircraft roll damping is 
increased significantly, and hence the roll mode time control is decreased to 
typically half the value due to the wing alone. 

The rotor contribution to C n g is always negative and, compared with the 
positive term from the vertical tail, the rotor term is small for the gimballed 
rotor but more significant for the cantilever rotor. An increase in vertical 
tail effectiveness is then required, especially for a cantilever rotor or, 
equivalently, the rotor decreases the aircraft Dutch roll mode frequency. 

The rotor contributes to C Ti ^ by a Ty term and a pivot moment term. To 
counter the rotor contribution (maintain the same net C nr ) , the Ty term always 
requires less tail; the pivot moment term requires more tail for the gimballed 
rotor because of the negative H force, but less tail for a cantilever rotor 
because of the hub moment. The Ty term dominates, though just barely for the 
gimballed rotor, so the net result is a requirement for less vertical tail to 
maintain the same value of Cn r , especially for the cantilever rotor. Equiva- 
lently, the rotor influence means an increase in the Dutch roll damping. 

The rotor contribution to Cio is typically equivalent to about 5° of wing 
dihedral (by the present estimate). The rotor derivative involved is a cross 
term - a vertical force due to a lateral velocity - so it is proportional to 
(Kp for the gimballed rotor), and the contribution to the aircraft deriva- 
tive is proportional to the rotational direction of the rotor ( sgn ft). The 
rotor contribution, in terms of the wing dihedral required to counter it, has 
the same sign as -sgn ft for the cantilever rotor, and the same sign as 
Kp sgn ft for the gimballed rotor. For the rotors considered here, more wing 
dihedral is required for the cantilever ( sgn ft = -1) and less for the gimballed 
rotor ( Kp < 0 and sgn ft = 1); or, equivalently, the rotor produces worse and 
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better spiral mode stability, respectively. Generally, the rotor contributes 
significantly to Cjo, the direction depending on the rotor rotation direction, 
whether a cantilever or articulated rotor is involved, and on pitch/flap cou- 
pling for the articulated rotor. Note that the large tip pylons usually asso- 
ciated with the tilt rotor configuration also contribute significantly to C 

The rotor contribution to Cy is negative and several times the negative 
contribution of the vertical tail. Hence it requires less vertical tail or, 
equivalently, increases the Dutch roll damping. The influence o f Cy Q iS 
usually considered secondary to that of Criy> > however. 

The rotor contribution to C n is an order of magnitude larger than the 
wing term; C np < 0 is adverse yaw; which is the direction of the wing contri- 
bution. Therb is a thrust term and a pivot moment term from the rotor. In 
comparison to the wing term (so that positive ratio means adverse yaw) , the T ^ 
term has the same sign as sgn 52; the pivot moment term has the same sign as 
sgn 52 for a cantilever rotor and the same sign as Kp sgn 51 for a gimballed 
rotor. The pivot moment contribution is small, however, especially for the 
gimballed rotor, and the T ^ term dominates. Hence the net contribution to 
t±Ci v /C lv has the same sign as sgn 51 and, with variations in V, it behaves as 
V sin <f (Tr/VCp). For the rotors used here, the rotor contribution for the 
gimballed rotor is adverse yaw and, for the cantilever rotor, it is favorable 
yaw. Since the rotor results in a coefficient larger by an order of magnitude 
than usual for aircraft, however, it is expected to significantly alter the 
lateral modes of the aircraft. 

The rotor contribution to C^, like that to C n p, is an order of magnitude 
larger than the wing contribution. Relative to the wing term (i.e., hC^/CiJ , 
the torque term has the same sign as -sgn 51, while the hub force term has the 
same sign as -sgn 51 for the cantilever rotor and the same sign as Kp sgn 51 for 
the gimballed rotor. The hub force and torque terms have approximately equal 
magnitudes; hence the sum would be small for an articulated rotor with Kp > 0. 
With the present examples then, the rotor contribution is in the same direction 
as the wing contribution for the cantilever rotor and opposite for the gimbal- 
led rotor. The increase in the magnitude of is expected to significantly 
alter the lateral modes of the aircraft. 

In summary, the influence of the proprotor on the aircraft dynamics is as 
follows. The static longitudinal stability (static margin) is decreased [C m ) ; 
the phugoid damping is increased ( Cx u )> the short period damping is increased 
somewhat by C z (but influenced primarily by C m q)> the roll damping is 
increased the Dutch roll frequency is decreased, especially for the 

cantilever rotor; and the Dutch roll damping is increased, especially for the 
cantilever rotor [Cn , and also Cy^) • Because of C^q, the short-period damp- 
ing and frequency (i'Te., maneuver margin) are decreased for the gimballed rotor 
(due to the negative H force damping) and are increased for the cantilever 
rotor (due to the hub moment capability). The rotor significantly influences 
the spiral mode stability (<?£«), the direction depending on the rotor rotation 
direction and also on vg and Kp. Three of the rotor contributions - 
and C ^ - are lateral/vertical coupling rotor coefficients, hence they are 

proportional to sgn 52 and also depend on Vg and Kp, the rotor coupling param- 
eters (ff*) . Three of the rotor contributions - C Xu> C n and C ^ - are an 
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order of magnitude larger than the usual aircraft derivatives and therefore 
may be expected to significantly alter both the longitudinal and lateral modes 
of motion of the aircraft. 

It is concluded that flapping rotors operating in high inflow contribute 
substantially to the aircraft derivatives. A numerical estimate of the rotor 
contributions to the aircraft stability derivatives, from the expressions given 
here, would probably not be adequate, primarily because rotor lag and wing 
motions were neglected. 

Stability Derivative Nomenclature 




Ny, N z 
u y a 

V> q, r 

^ w 
h 

d w 


H 

a t 


h) 

Oy 

sgn Q 

q 


aircraft longitudinal, lateral, and vertical forces (subscripts 
x, y } z for coefficients) 

aircraft roll, pitch, and yaw moments (subscripts £, m, n for 
coefficients) 

aircraft longitudinal, lateral (sideslip), and vertical (angle 
of attack) velocity perturbations 

aircraft roll, pitch, and yaw angular velocities 

wing span 

mast height (aircraft center of gravity to rotor hub) 
wing chord 
wing area 

aircraft lift curve slope 

aircraft lift coefficient 

horizontal tail area 

horizontal tail arm 

horizontal tail lift curve slope 

vertical tail area 

vertical tail arm 

vertical tail lift curve slope 

rotor rotation direction 

in coefficients, dynamic pressure (1/2) P^ 2 

neutral point 
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Stability Derivative Nomenclature (Concl . ) 


h m maneuver point 

T wing dihedral 

The aircraft stability derivatives are defined as they are introduced. 

SECTION 4: NINE-DEGREE-OF-FREEDOM MODEL FOR A PROPROTOR ON A CANTILEVER WING 


Wing Equations of Motion 


In this chapter, the equations of motion for a proprotor on a cantilever 
wing are derived. The equations of motion and the hub forces and moments were 
obtained previously for the six rotor degrees of freedom (blade collective and 
cyclic flap and lag motion), including the excitation by the six degrees of 
freedom of the shaft motion (eqs. (44) to (48)). The equations of motion for 
the wing elastic bending and torsion motion are now derived. The wing degrees 
of freedom are forced by the rotor hub forces and moments, and each mode of 
wing motion produces a corresponding shaft motion. Thus the two sets of equa- 
tions, for the rotor and for the wing, are combined by substituting the wing 
degrees of freedom for the shaft motion in the rotor equations of motion and 
forces, and then substituting for the rotor forces in the wing equations of 
motion. The result is a coupled set of equations that describes the 
aeroelastic behavior of the proprotor and wing system. 

The proprotor aircraft configuration consists of large diameter flapping 
rotors mounted on the wing tips. Usually, the engine and transmission are also 
mounted in the pylon at the wing tip; hence there is a large mass and inertia 
at the wing tip. The rotor hub forces are transmitted to the wing tip through 
a mast of height h (i.e., forward of the pylon pivot). 

The dynamics specific to the proprotor configuration are of primary 
interest here, that is, the high inflow rotor dynamics coupled with the wing/ 
pylon motion. Hence the model considered is restricted entirely to the air- 
plane cruise mode configuration, with the shaft parallel to the free-stream 
velocity. The model is restricted also to the frequency range most important 
to the coupled wing and rotor motion. Only the lowest frequency wing modes 
are considered, and the elastic motion of the pylon with respect to the wing 
tip is neglected since it usually has a much higher natural frequency than the 
lowest wing modes. The motion considered then is the wing elastic motion with 
the pylon and the rotor shaft rigidly attached to the wing tip. The aircraft 
rigid-body motions are also neglected since they are degrees of freedom of 
low frequency and are not highly coupled with the motions to be examined. 
Neglecting the rigid-body motion places the primary emphasis on the basic rotor 
and wing dynamics, as is desired here, and so may be justified as an appropri- 
ate first step at least. Indeed this rotor and cantilever wing model is very 
useful in establishing the behavior of the proprotor aircraft. Hence a 
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cantilever wing model is used for the fixed system motion, with the 
understanding, of course, that eventually a more complete model including the 
aircraft motion must be used. In addition, the cantilever model corresponds 
to the configuration of many proprotor models that have been tested in wind 
tunnels, including the full-scale tests considered later. 

The model considered consists of a wing and rotor in cruise flight 
configuration, operating in a free-stream velocity V with the shaft always 
parallel to V so that the rotor equilibrium flow is purely axial. The rotor 
operates in high inflow. The wing root is attached to an immovable support 
with cantilever root restraint. The wing motion consists of elastic bending, 

\ ertical and chordwise, and elastic torsion. A pylon with iarge mass and 
moment of inertia is rigidly attached to the wing tip. The rotor is mounted 
on the pylon with the hub forward of the wing elastic axis, with the rotor 
shaft horizontal (parallel to V) . The rotor has three or more blades, with 
first mode flap and lag motion for each blade. The rotor hub forces and 
moments are transmitted through the pylon to the wing tip. 

The wing is assumed to have a high aspect ratio so that strip theory can 
be used for the wing aerodynamics and engineering beam theory for the elastic 
bending; this assumption is well justified for the tilt rotor aircraft, which 
ave, typically, an aspect ratio around 6. Wing sweep, dihedral, and angle of 
attack are considered, but the major effect is that of the wing sweep on the 
position of the effective elastic axis of the wing, hence on the effective 
mast height for the transmission of the rotor hub forces to the wing bending 
and torsion motion. Regardless of the wing sweep, dihedral, or angle of attack, 
the rotor shaft is assumed to be maintained parallel to the free-stream veloc- 
ity in equilibrium trim flight. This assumption is required to avoid periodic 
coefficients in the rotor equations of motion due to an inplane component of 
the trim velocity. 

geometry and motion - The wing geometry is defined by a straight spar 
line that is the locus of the local elastic axis. The wing root is supported 
with cantilever restraint; the rotor shaft is attached rigidly to the wing tip. 
The wing geometry is shown in figure 11. The wing has a constant chord and 

a length y T ^ from root to tip (semispan). The distance along the spar is w „ 
measured from the root. The shaft has length h, the distance the rotor hub is 
forward of the wing tip elastic axis (mast height). The wing spar is basically 
perpendicular to the forward velocity V, but small wing sweep, dihedral, and 
angle of attack are considered. The wing root is attached to a plane defined 
by the forward velocity V and the vertical; then three rotation angles define 
the orientation of the spar with respect to the free-stream velocity: dihedral 

(positive for upward rotation of the wing tip), sweep (positive for 
sweep aft), and angle of attack <$ U2 (positive nose-up). All 3 these angles are 
assumed to be small, an appropriate assumption for proprotor aircraft in cruise 
The rotor shaft must then be rotated by the angles -6 Wl , -6 Uo , and -6 Uo to keep' 
the shaft parallel to the free-stream velocity. z 6 

The rotor (fig. 3) is placed on the tip of the wing (fig. 11) so that the 
wing tip motion is transmitted directly to the rotor shaft. It is also neces- 
sary to account for the rotational direction of the rotors. The aircraft has 
two counterrotating rotors, one on each wing tip, but the rotational direction 
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of the rotor on the right wing (fig. 11) may be either clockwise or counter- 
clockwise. The equations of motion and hub forces for the rotor are derived 
for the rotor model in figure 3, that is, assuming clockwise rotation. This 
model is always used for the rotor, and then the rotational direction of the 
rotor is accounted for by placing the rotor on a right-hand wing (fig. 11) or 
on a left-hand wing as appropriate. The influence of the rotational direction 
is a number of sign changes in the equations of motion, which reflect how the 
hub forces and moments of the standard rotor (fig. 3) excite the right or left 
wing, and how the right or left wing produces motion of the rotor shaft. 


The notation fi indicates the influence of the rotational direction of the 
rotor. The clockwise rotating rotor on the right wing (figs. 3 and 11) is 
denoted by n = 1; the clockwise rotor on the left wing (fig. 3 and the mirror 
image in fig. 11) is denoted by SI = -1: 


n = 


(+1; rotor rotation clockwise on right wing, 
counterclockwise on left wing 

-1; rotor rotation counterclockwise on right wing, 
clockwise on left wing 


(174) 


The wing motion is described by elastic bending and torsion of the wing 
spar; the wing displacement is shown in figure 11. The pylon (and with it the 
shaft) is rigidly attached to the wing tip. The existence of an elastic axis 
of the wing (assumed to be a straight line) means that the wing distortion may 
be described first by elastic torsion of the wing about the local elastic axis, 
without bending the wing; and then by elastic bending of the spar, which 
deflects the elastic axis from the undistorted position without changing the 
torsional deflection. 


A modal description of the wing elastic deformation is used, and only the 
lowest frequency modes are retained. The elastic torsion of the wing results 
in a pitch change of the local wing section (positive nose-up as 

shown in fig. 11). With a modal representation, this motion is written 


'l, 


which is an expansion of 0^, in a series of the mode shapes of the elastic 
torsion motion. The generalized coordinates vp are the degrees of freedom. 

The modal representation is useful because it separates the time and space 
dependence of V Associated with each degree of freedom pp, there is an 
equation of motion with appropriate generalized mass and stiffness, hence a 
natural frequency of each mode. Only the lowest frequency degrees of freedom 
are retained in this study of the basic dynamics; it is sufficient, in fact, 
to consider only one wing torsion mode. Then, 

6 W = C1?5) 
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d e scrib es the motion If the mode shape is normalized to unity at the tip 
WjO = 1, then p gives the nose-up torsion angle in radians' at the tij 

with elas ^ lc bending of the wing results in the deflection of the wine soar 
with components both perpendicular to the wing surface (Vertical nr L! • P 
ben ,ng) and parallel to the wing surface (cho^dwS bend " "‘t de lotion 

Of the spar line normal to the wing surface is - (f u 1 ,w : j Election 

upward deflection. The deflection^, the pllne'tf ihe“iing l s j (TV?* ” 

f and r POSltl d% 0r a rear ” ard deflectiOT - vertical and chorSwiVbiidin* 
of the ScUon iiitt t0 th< \ direct ion of the local principle axes 

the wing spar "hut t[ P ! built ~ in Wlng twist > these a *es are the same all along 

sweep, a. 

modes^ar^retained? ' For^^al'ySs l^suff'"- ^ ^ 

mode each for the z w and x u represcntahons; hence' 01 ™ t0 retal " °" ly °" e ' 


z w = 


4 1 hy Q/y) | 

X W = q 2 0)hyQ/ w ) 


( 176 ) 


where n w is the mode shape of elastic bending of the wing. For the present 

Ep z 2 lisl £ ssE? 

SHS? “ n~r 

represlnt" S ^Mhe^^ for 


The degrees of freedom that represent the wing motion are thus: 

p, wing elastic torsion, positive nose-up (p = at the tip); 

q v wing vertical or beamwise bending, positive upward = z w /y T at the tip); 

q 2 , wing chordwise bending, positive rearward ( q 2 = x w fy T at the tip) 

Associated with these degrees of freedom are mode shape £(» ,) for torsion and 
hyO/y) Tor bending, which are normalized to £ and urn resDectivelv at «.u 

f^i° nSi iw er tbe r mo ^ ion of rotor shaft in terms of the wing degrees of 
freedom. The shaft displacement and rotation fxn w § g f 

point h aft of the hub are required because oAhe MUoJ’.t'thXi^tJi “ 
whtch is specified by the wing degrees of freedom (q, , q„ and p,! ?f ^1 
gleets for the moment the effects of the rotor rotation direction and wine 
sweep, dihedral, and angle of attack, the following shaft k produced 

The wmg torsion deflection p results in shaft pitfh a„. TOe ™ng venjcal 

also'produce^a'slope ^ 

ur tne elastic axis at the tip, it results in shaft roll 
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The wing chordwise bending q z results in longitudinal displacement z p 
iiie wing O / i,r 4 *-v» rtic mrif 1 r» 1 . thfvre IS no fit 


rt he wins tnoruwisc ucnumj mo -- ° ... • . 

of the shaft and also shaft yaw angle a x . With this model , there is no fir. 
order source of shaft lateral displacement y P . _ If the magnitude of the dis- 
placement and rotation at the tip due to the wing degrees of freedom (given by 
the mode shapes) are accounted for, the shaft motion is (for a clockwise ro . 
ing rotor on the right wing, figures o and 11, with - 

“x = -<7? V (y T J 




J «3 


■1 W 'V-L w' 

a y = = p 

a z = 

Xp = = q \ PTi 


(177) 


UP 


= 0 


- 


W 


- q ?VTw 


Note that the wing bending motion produces coupling of the longitudinal and the 
lateral/vertical groups of the rotor, q y and q 2> giving both longitudinal 
motion of the shaft (s p and a 3 ) and lateral/vertical motion (x p and a ). The 
coupling is not strong however, and it is found from the behavior of the sys- 
tem that wing chordwise bending q „ is basically a longitudinal motion, and 
vertical bending and torsion (^ and p) belong with the lateral/vertical group. 

Consider the clockwise rotating rotor on the left wing, that is, fl = -!• 
the only change is in the direction of the shaft angle due to the slope of the 
elastic' axis at the tip during bending (consider the mirror image in fig; U. 
including the definitions of the bending and torsion deflection) . rt this 
change in sign is incorporated by use of the 9 notation, 






= (y T J ) 


(178) 


and the rest of the shaft motion in equations (177) is unchanged. If the wing 
motion produced a lateral shaft displacement y p , that, too, would change sign 
with the direction of rotor rotation. 

Consider the effect of wing sweep, dihedral, and angle of attack. The 
wing tip displacement and rotations, along with the wing motion, are de me 
with respect to the wing spar and the section principal axes - which are 
rotated bv 6,, 6,, OJ and with respect to the wind axes. Hence the wing P 

motion will hive a slightly different decomposition into the shaft motion 
which remains in the wind axes. For example, vertical bending q^ produces, 
addition to vertical displacement x p and shaft roll cm as given previously, 
some shaft pitch cm due to the wing sweep some lateral displacemen y p 

to the wing dihedral 6., , and some shaft yaw and axial displacement z p due 

to angle of attack S^. 1 After a consideration of the complete set of wing and 
shaft motions, the result is 
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a ;r = -q 2 ^x'(y T J + p^wi - ^x'(y Tw ) 
a y = p + - q l ^w 3 ^ w '(y Tw ) 

“ 3 = - P ^ 3 + ^ 2 ^ 2 %'^) 

X P = ^ ■ ^Vr a 
i/p = 'W^ - <? 2 fi S% u 

2 p = ~ q l y T w ~ q \ & w 2 }Jt w 

Note that the effect of the rotational direction of the rotor Q 
change the signs of c^, a 3 , and y p . 
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is only to 


than 7116 inf * uence °f the wing sweep (and dihedral) requires more attention 

elastic ^xi^r 6 ' C °" sider . the “"swept wing, represented by a straight unswept 
lastic axis line with cantilever restraint at the root. This structure has P 

Ji^rh characteristics: elastic torsion at an inboard wing section results in 

placement n of S the 0atb ° ard sect i° ns > but produces no vertical or chordwise dis- 
ann^H t • StlC aX1S fr ° m itS undistorted position; and a force 

PP . d 1 the V ' ang tlp at the elastic axis results in bending of the wing but 
due tTlhTs force 0 ” ^ tl0n " ince . there is no torsion moment about any section 
and torsion motions. ® 15 n0 6laStic C ° Uplirg of the wia S bend ing 


restr^i nf he ng . ls now swept ’ that behavior would be maintained if the root 

shaft ^' a a fu Pt ‘ In that CaS6j the description developed for the 
shaft motion produced by the wing torsion and bending would be correct includ 

s rof thft 5 ° swe H e v nd dihedrai - However * tMs is ^ way s; e ;r 

mgs Cof the type used for proprotor aircraft) are built. The wings are 
usuaiiy built with a center box structure in the fuselage, where thf spars are 
nswept, and only the wing structure outside the fuselage has swept spars. 

fuselage ^ seve ^ al points > where the wing box is tied to the 

useiage. (The wing used m the full-scale, wind-tunnel tests considered 

In ,!L WaS 3 j° bUllt thlS way,) 0116 a PProach to treating such a structure is 
use a good structural dynamics analysis to calculate the coupled bending and 

restraint t°and ’TV"" PylOT ' *“*"«"* influence oTJhl roof® 

estraint and sweep. Such an approach is useful if available and, in fact it 

pro a y necessary if an accurate representation of a specific design is 

r™tatio?th» a " al rr iS r n desired 1™. 

Sin(e this (eno^ a eS .° y ' ele, " e,>ts *>« fundamental to the behavior. 

rather JhantEeL. r 3 ? 3neral rumination of proprotor dynamics 

sign of a specific vehicle (with a swept wine) such a 
representation is adequate. P gJ ’ Cn a 


The model used to represent 
except for a bend at span station 


a swept wing has a straight elastic axis line 
y w = yB w ’ where the wing sweep and dihedral 
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are entered. The root restraint and the wing inboard of y B axe unswept, while 
the outboard section is swept. Hie mast height h is measured from the wing tip 
elastic axis to the rotor hub. The following behavior is expected of such a 
model typical of swept wing behavior: elastic torsion of the inboard sections 

produce s^a pitch change at which then rotates the entire outboard section 

about the extension of the inboard spar line, thereby producing displace 
of the outboard elastic axis from its undistorted position, and a force (verti- 
cal for the sweep effect) at the tip elastic axis produces a torsional moment 
about the wing sections inboard of y B thereby producing torsion motion of the 
wing rather than just bending as for ¥he unswept case. The torsion and bending 
motions of the wing are then elastically coupled. The first eff ecu in partic- 
ular, means that torsion of the wing now results m the displacement of the 
wing spar, which is produced only by bending for the unswept wing. Specifi- 
cally, the torsion motion produces now displacements at the tip - e ^ c< p ° e 
shaft - as well as rotation of the shaft. An advantage of this model is that 
the shaft motion due to the wing degrees of freedom can be obtained by simple 
geometric considerations, and the wing equations of motion can be obtame y 
a simple extension of the methods (based on section force and moment equilib- 
rium) used for an unswept wing. Both are the result of retaining the 
sentation of the wing by an elastic axis line. The wing motion is defined then 
by elastic torsion about the local elastic axis (producing local elastic pitch 
changes given by the mode shape ^ , and elastic axis deflections due to the 
sweep), followed by elastic bending of the spar (producing local spar disp ace 
ment given by the mode shape n w , from the displacement due to torsion) . The 
effect of the bent elastic axis is to couple the influence of torsion an 
bending; with sweep, the tip is displaced due to torsion, and also the torsi 
mode is excited by forces at the tip elastic axis. 

A swept wing is characterized by an effective elastic axis for the 
vertical bending of the tip: at some point on the shaft or its extension the 

application of a vertical force results in purely vertical displacement of the 
shaft, with no rotation. Without sweep of the wing, this point would be at 
the wing tip elastic axis; but with sweep, a force there produces a nose-down 
pitch motion of the shaft also (for aft sweep). Hence the effective elastic 
axis is some distance forward of the wing tip elastic axis. This feature is 
included in the model used here. Generally, the effective elastic axis ( or 
the tip) lies between the actual wing tip elastic axis and the extension o 
the inboard (unswept) spar line - the actual position depending on the degrees 
of root restraint and sweep and other structural details. 

The model used here is, in fact, only conceptual. It is used because it 
allows an elementary derivation of the equations of motion, including the most 
important features of the swept wing behavior. The bend in the spar is no a 
feature that can be determined accurately from the geometry of the wing. Hence 
u R , the spanwise location of the spar bend, is just a parameter that charac- 
terizes the influence of the wing sweep on the coupled bending and torsion 
motions of the wing, which is determined by the details of the wing struct ura 
construction, geometry, and root restraint. This parameter (it will be found 
that what is really required is i w (y B J) can be identified by matching the 
behavior of the model to that known from a better model or from actual exp 
ment. The most important feature of the swept wing is the effective elastic 
axis position at the wing tip. If that is known from a good structural 
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analysis of the wing, or measured experimentally, then the corresponding value 
? yB w wh * ch , glves .that effective elastic axis position for the model used 
ere may be determined. Then that single parameter completely determines all 
the characteristics of the swept wing in this model. 

a Z 111 ® inf | uen f of the bent elastic axis model on the rotor shaft motion 
due to the wing degrees of freedom can now be determined from simple geometric 

tiCnrf r ?h tl0nS - "‘ ng . pi J Chd " fl « tlon « rotates the entire outboard p or- 
tion of the wing about the inboard spar line? so torsion of the wing produces 

at Wing tip ' vertica Hy due to sweep of the wing and longi- 
tudinally due to dihedral. The rotation angle is d w (y R ) = pg r VR ) and the 

: rm - ] L Pip 1S s y? w S W3 f to first order for small 6 W and 

in W the W shaft vert^* ~ 1 5 VTw&Wi for dihedral. The increments 

in the shaft vertical and longitudinal displacement are then 


AX P = -PW y Bj y TJ» 3 
w w ° 

A3 P = ~ P ^w ^ y B, 


(180) 


'W 3 / 

The complete equations for the shaft motion due to the wing degrees of 


freedom are 
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The length k EA = h - is the distance from the effective tip elastic 

a(is d?sn?Lle»°l ) Ub V" *1 ■ 'f rIy ' for ‘ Wt > the vertical elastic 

and/ f!r I! For convenience, n w ’ u is written for V (y T ) 

h) * ZO ^ ^ " 


f jl in g Equations of Motion - The equations of motion for the wing degrees of 
freedom, lowest mode elastic bending and torsion, are 

{I qw + m PyT w y<i\ + C q q x + K q x + S w p = M q + M q 

1 l aero 1 1 


1 vo tor 


%li 2 * C q jt 2 * K qi q 2 - S w I p = U ♦ Hq _ 

z z aero z rotor 

+ % y P + c fP + K pP+ ««?, - Vu 2 5 2 ■ Xp aevo * «P : 


-rotor 


( 182 ) 
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for the vertical bending, chordwise bending, and torsion motion, respectively. 
Each equation has inertia, structural damping, and structural spring terms, 
forced by wing aerodynamic forces and by the rotor hub forces and moments act- 
ing on the wing tip. The generalized inertia of the wing bending mode is 

x = rn w n/ dy w , and for the torsion mode, I pw = j W d V w > where 

Hw 0 
m,s is the mass per unit length of the wing and I 0 is the wing section moment 
of inertial. To these wing inertias are added the pylon inertia terms: m p is 

the pylon mass (without the rotor) and Ip x and I Py are, respectively, the. 
pylon yaw and pitch moments of inertia about the wing tip effective elastic 
axis. The inertia coupling of the bending and torsion of the wing is due to 
the offset of the pylon center of gravity: S w = rn P y Tl /-PEA » where Z PEA 1S the 
distance the pylon center of gravity (without the rotor) is ahead of the wing 
tip effective elastic axis. The structural spring terms are written 


K, 


r yT w 

'l "I 


EI ^ 


C K 




r y ?W 2 

- J d y w 




r w 2 

K P = I 

•/n 


where EI ZZ is the section beamwise modulus/inertia product and, similarly, for 
the chordwise and torsion elastic restraint. These expressions are not used 
here, however, since a very accurate estimate of the bending and torsion mode 
shapes would be required to evaluate them. 

The derivation of these modal equations follows the standard methods of 
aeroelastic analysis. For typical proprotor configurations, the jjylon mass is 
so large that it dominates the wing inertias, that is, < ' < ™PyTw . . 

I « Jp . Hence the inertia is primarily that of the pylon and rotor, with 
the wing Contributing only the elastic restraint of the motion. This is 
fortunate in that calculating the wing inertias requires an accurate estimate 
of the mode shapes, while the pylon mass and moments of inertia are well- 
defined characteristics that are easily determined. The wing structural spring 
constants (K a , 3 K q „ and K p ) are best determined by adjusting their values so 
that the predicted frequencies of the modes match the frequencies measured 
experimentally. By this procedure, the wing inertias and structural parameters 
are determined from the characteristics most well defined and easily measured: 
pylon inertia and natural frequencies of the modes. It is particularly impor- 
tant to match the measured frequencies well since they have the most important 
influence on the dynamic behavior of the system. This procedure is satisfac- 
tory here, but, for the preliminary design of an actual vehicle, a good 
structural analysis method for predicting the generalized masses and stiff- 
nessess of the wing and pylon modes is necessary. It is evident that (for the 
lowest wing modes at least) the wing mode shapes have a secondary influence on 
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the equations of motion. A very rough approximation to and £ is 

satisfactory to estimate the wing contributions I q and I p to the inertias 

since they are dominated by the pylon contributions; the structural spring 
constants and other parameters such as n y '( y T ) are determined by matching 

them to the measured characteristics of the wing (or those calculated by a 
more accurate method) rather than using their definitions in terms of the wing 
mode shapes (which are not really available). 

The wing structural damping constants (C p ^ , Cq 3 and C p ) are determined, 

as the spring constants are, by matching the theoretical results to the meas- 
ured^ characteristics . By definition, the damping constants are given by 
r ~ A 3s with the appropriate inertia and spring constant, where g s is the meas- 
ured structural damping coefficient of the wing (twice the fraction of critical 
damping). The wing structural damping is typically 1 or 2 percent of the 
critical damping. 


The wing motion is excited by the rotor forces and moments acting at the 
hub. For the moment,. if the effects oi r the rotor rotation direction and of 

the wing sweep, dihedral, and angle of attack are neglected, the rotor forcing 
terms are 
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Mp = M 

votov 


+ hH - C* y m Tq 
PTT w 3 
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f -w 

where C p ^ = J d yw^T w S 2 / 3 * The excitation of wing vertical bending 

^1 is due to rotor torque and vertical force, with the effectiveness of the 
former determined by the slope of the mode shape at the tip n y ' (yj< ) (written 

as q w ' in eqs . (183)) and for the latter, by its displacement ) = y T . 

Similarly, the excitation of chordwise bending q is due to the pivot yaw 
moment ( M x - hY ) and the thrust force. The wing torsion motion is excited by 
p, the pivot pitch moment (My + hH) , and by the trim thrust. The thrust term 
results because wing vertical bending q elevates the rotor trim thrust above 
e m oard sections and so gives an arm about which the trim thrust produces 
a torsion moment; the constant involved, C* , was evaluated using the 

approximation £ y = y w /y and n y = ,y w 2 /y^ for the wing modes. 


Introducing the influence of tlie rotor rotation direction, that is, 
putting the rotor in figure 3 on a j'eft-hand wing, simply changes the signs of 
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the moment contributions to the wing bending excitation. If the convention _ 
for Q (given previously) is used to indicate these sign changes, the result is 


Ak 


1 




rotor 


+ Ij rp h 

1 w 
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M - -fin 1 - hY) - y? 1 

rotor 


and M-n remains unchanged. 

r rotor 

If wing sweep, dihedral, and angle of attack are considered, the first 
influence is a slightly different decomposition of the rotor forces, which are 
defined with respect to the shaft (wind) axes, into the wing tip axes for 
determining the excitation of the wing bending and torsion. With the 
appropriate coordinate rotations, the result is 

M q, = [Q ~ QS w 3 ( M y + hH) ~ S CMx ' hY)] ' 

ro tor 

+ yy> ( H - Y - ? ? ) 

W L c. 


M = fin,' [~M X + hY + Q6 W (H u + hH) - 

4 O u/ 1 £- 


ro tor 


+ y T 0 -t - 

W it. 
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M = M u + hH + Q + (M a . - hY) 

rotor 
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C * 

pq 


y, r T{-q 

-u 
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q 2 ) 


However, these expressions are for a wing that is swept along its entire length, 
including the root restraint. The model considered here has the inboard por- 
tion unswept, and only the wing outboard of = y B ^ has sweep and dihedral. 

With this model, there are also the additiona., terms. 
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1 rotor 
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2 rotor 
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rotor 
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3 1 1 w I 

The terms in equations (186) (except for last two in the torsion moment) simply 
account for the decrease in the effect of the and 6 W3 components of the 
rotor forces and moments in exciting the wing when not all of the wing has 
sweep and dihedral. These terms combine with the similar terms in equations 
(185) to produce a reduced net influence. For example, sweep of the wing 
produces a component 6 Ws M y of the rotor pitch moment that tends to produce 
vertical bending motion, that is, it contributes to M q . The correspond- 

rotor 

ing terms in equations (185) and (186) combine to give 


* «w 3 y (y Bw )My = -SwjVG/^) - 

which illustrates the decrease in effectiveness due to y B . The other terms 

are similar, a straightforward decrease in the effectiveness of the effect of 

the sweep and dihedral because they occur only on the outboard portion of the 
wing. 


The remaining terms in equations (186), the torsion moment due to H and T, 
are the effective elastic axis influence. The torsion moments due to the hub 
vertical force H (from eqs. (185) and (186)) combine to give 

m Protor = ft - ]#• From this, it follows that a vertical force 

on the shaft at h = ^ w (y Bw )S u ^y T ^ produces no torsion, and so, by definition, 

that is the location of the effective elastic axis at the wing tip. Now 6 ,, y 
is the distance of the wing tip behind (for aft sweep) the root elastic 3 

axis, and S y Q/ gy ) is a constant less than 1 (of order y Bw /y Th} ); thus the 

effective elastic axis is moved forward of the wing tip elastic axis, a frac- 
*"^ on toward the root spar line. The distance from the wing tip 

effective elastic axis to the rotor hub is then 


h EA = h ~ ( 18? ) 
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where h is the actual distance from the wing tip spar to the hub. The 
corresponding vertical displacement of the effective elastic axis, due to the 
dihedral, is 


Z EA ~ 


( 188 ) 


w 


(Zero geometric vertical displacement of the shaft from the wing tip has been 
assumed. ) 


The parameter K (j/d ) can be evaluated from a knowledge of the effective 


w 


elastic axis position. If Ah is the distance the effective elastic axis is 
ahead of the wing tip spar (A h = h - Ti^) , then 


Ah 
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The other parameters required can be estimated in terms of K w (y B ) as 
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Only a single parameter remains to be determined for this model, £ (j yg ). It 


W 


is determined by matching to the correct effective elastic axis position, hence 
the most important feature is correct and this model should provide an adequate 
estimate for the other influences of sweep and dihedral. 


The complete expressions for the rotor forcing terms in the wing 
equations of motion are then: 
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where n u ' is written for nw' (i/y ) and for 
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The wing equations of motion are normalized by dividing by (N/ 
this normalization is denoted by superscript *, as in 7 /(N/2')I b = I* and 

similarly for the other inertia, damping, and spring constants. The rotor 
forcing terms then take the usual rotor coefficient form. The wing aerodynamic 
terms (presented in the next section) are normalized by dividing by 
y(N/2)I^ = (i rr/2)aa , which is also denoted by superscript *. This introduces a 

factor y for the wing aerodynamic terms, as for the rotor aerodynamics. The 
equations of motion for the wing vertical bending, chordwise bending, and 
torsion modes are then: 


J* + m/ 

H w 


0 

T q 

* %*£ * p ■* 

- 

<7 * 

b w 


~^W*^W 2 

\ 


0 0 

Vh\ r 

K q 


0 

0 


% 


Cp * J 


V 


W Id 2 


W 


L P J 

y 



K * 

$2 


V/ \S^ T Jd2C T lod) -C* q y T Jy2C T /aa)& w ^ 


0 

0 


K p *. J 




+ Y 


-On/Ofiy (1 - £,) 


W 


1 

-1 


w 


- Qr] w S 

-fin ' 

w 

-«««. Cl - UJ 


w 


- 2 C M /a a' 

y 

2C Mx /oa 


+ Y 


y T w ~ f26 w 3 ?J C 1 - 

_<Sw 2 y T w + ‘ V 04 *,** 1 - 5 W ) 


E4 


'W * VV 1 ' ? » 2ri 

* i/j^d - c/) 

a y f 1 - v* 




1 

K> 

1 



+ y 

- 2 W 

Cg/aa + y 

CM 

1 


-“V 1 - V- 


- 2z ea -J 


C T /aa 


(192) 


An exception to the normalization is the definition of m p *: 
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m p * = 


It also follows that S w * = S w / (N/2)I b = m F * zp M / y T ^. The Lock number y appears 

as a factor of all the aerodynamic terms (the rotor forces and moments have 
inertia terms, too, but always with a factor y” 1 )- This single parameter 
accounts for the relative influence of the aerodynamic and inertia forces; 
specifically, it is the only parameter that varies with air density p, all 
other constants being the ratio of inertias (that is, of course, the reason or 
the normalization by (N/2)I b ). The spring and damping constants are 
K * = K q /(N/2)I b and C*^ = C q J{N/2)I b (and similarly for chordwise bending 

and torsion) . Since these constants have already been made dimensionless, in 
terms of the dimensional values the normalized spring and damping rates are 
K* = K q /(N/2)I b V 2 and C* = C q / {N/2)I b \a\ . The effect of a variation m 

rotor rotational speed is to change the relative spring and damping rates _ of 
the wing, because the wing natural frequencies are really constant dimensional 
values (cycles/sec). The wing frequencies then change with 51 when expressed 
as dimensionless quantities (per rev). This is, in fact, the only influence 
of 51 on the equations of motion, besides a possible change in the rotor natural 
frequencies v R and v r with 51. In the inertias, as in the spring and damping 
constants, it is not necessary to introduce the intermediate^steps of making the 
quantities dimensionless (using p, 51, and R ') ; for example, I Fp - Tp^/(/l//2) b 

is correct with dimensional quantities on the right-hand side. The rotor 
radius R enters only in normalizing lengths, such as y T ^ and and, of 

course, the air density p appears only in y- 



where F z is the vertical aerodynamic force on the wing section (lift); F^, 

the chordwise force (profile and induced drag); and M w , the aerodynamic moment 

about the local elastic axis. The section forces F z an d F x are defined with 

10 10 

respect to the section principle axes, not with respect to the free stream. 

The integrals of these section forces over the span, weighted by the appropri- 
ate mode shape, give the net forces that excite the wing degrees of freedom. 
The velocity seen by the section has perturbations that result from the wing 
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degrees of freedom and from aerodynamic gusts; any interference with the rotor 
is neglected. From the velocity perturbation, the perturbations of the section 
forces are found (following a procedure similar to that used for the rotor aero- 
dynamics), and hence the wing aerodynamic coefficients (the trim terms are 
dropped) : 
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(Eqs. (195) continued on next page.) 
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where C T and Cn are the aircraft trim lift and drag (profile and induced) 

u o u o 

coefficients; CV and C D are their derivatives with respect to angle of attack, 
ot ct 

The section moment characteristics are given by x^ y the distance the aerody- 
namic center is behind the elastic axis, and c ^ the nose-up moment coeffi- 
cient about the aerodynamic center. The coefficients can be corrected for 
unsteady airfoil effects (by use of a lift deficiency function) and for the 
increased dynamic pressure in the rotor slipstream; such corrections would be 
small, however, and probably would not be consistent with the accuracy of the 
methods used to derive the coefficients. The constant dy^ = y^f {-noa) 

accounts for the difference in the normalization of the wing and rotor 
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coefficients. The dependence on forward velocity V is shown explicitly, except 
that for a fixed gross weight; Cr and Cj) are proportional to 7“ 2 while Cr 

o ot 

and are constant (at high speed at least, where the induced drag is small). 

The constants e n and f n are integrals of the wing mode shapes, which account 
for the way the motion produces forces on the wing: 
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The constants were evaluated using the approximate mode shapes r)^ 5 

a ? d 5 y w^ T w ( which are reasonably close with the large mass on the wing 


The derivation of the coefficients follows the standard techniques of 
strip theory in aeroelasticity, similar, for example, to the derivation of the 
contributions of the wing to the aircraft stability derivatives. The most 
important wing aerodynamic coefficients are the vertical bending forces due to 
direct angle-of-attack changes: and c ^ an £ e an gle 

of attack due to wing vertical bending, torsion, and vertical gusts is 
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The lift perturbation is 6F Zw = 6L = {l/2)V 2 c w C L J>a, and the moment exciting 
vertical bending is 
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If one normalizes by dividing by 0/2 )cra, coefficients C* , °q l qi anc C q 

are identified (as given in eqs. (195)). The remaining coefficients are 
derived in a similar fashion. 


Equations of motion for proprotor and wing- All the elements are now 
available to construct the equations of motion for the proprotor and cantilever 
wing system: the rotor equations of motion (eqs. (44) and (47)), the rotor hub 

forces and moments (eqs. (45), (46), and (48)), the shaft motion due to the 
wing degrees of freedom (eqs. (181)), and the wing equations of motion (eqs. 
(192)). It is only necessary to perform the matrix multiplications required 
because of the substitutions. The result is a set of linear ordinary differ- 
ential equations for the nine degrees of freedom: 

8 cyclic flap (longitudinal tip path plane tilt) 

10 

8 cyclic flap (lateral tip path plane tilt) 

10 

r cyclic lag (lateral rotor center-of-gravity offset) 

10 

x, cyclic lag (longitudinal rotor center-of-gravity offset) 

10 

8q collective flap (coning) 

C 0 collective lag (or rotor speed perturbation) 

wing vertical (or beamwise) bending 
<72 wing chordwise bending 

p wing torsion 

with the inputs 
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01 cyclic pitch (lateral control plane tilt) 

0 19 cyclic pitch (longitudinal control plane tilt) 
0 Q collective pitch 

Uq longitudinal gust 

g lateral gust 

(j 

vertical gust 


The matrix multiplications could be left to the machine, but it is better done 
by hand since there is considerable cancellation of terms in the coefficients 
of the equations of motion. The final expressions for the coefficients remain 
rather complex, however, because of the many (small) terms that involve the 
wing sweep, dihedral, and incidence angles. 


For the autorotation case, the collective lag degree of freedom becomes 

the rotor speed perturbation in the rotor model used here by setting 

I* = I* =1 and v r = 0 (as discussed earlier). Setting the collective lag 
C 0 ct 

natural frequency to zero assures that no torque is transmitted to the wing, 
but it is probably better to directly drop the Cg/aa forcing of the wing 
motion from equations (192) (primarily forcing of ^ ). 


Structural damping of the rotor blades is also added to this model; 
however, since the blade flap and lag damping are high already because of the 
high inflow aerodynamic forces, the low structural damping of the rotors con- 
sidered here is not very important to the dynamics. A term is added 


to the rotating lag equation of motion, and I^g Q 


- 1 ) 


l/2 


is added to the 


rotating flap equation (the structural damping does not act on the centrifugal 
spring term in vg) . The transformation to the nonrotating degrees of freedom 
and equations of motion follows as usual. Different damping coefficients are 
allowed for the rotor cyclic and collective modes of the rotor, specifically 
to account for the collective lag mode in autorotation (which must have zero 
structural damping as well as zero spring) and the coning mode of the gimballed 
rotor. The structural damping parameter g s is twice the fraction of critical 
damping, which is typically 0.5 to 1 percent for the cantilever rotor blades 
considered here. 


Simplified equations- The coefficients of the nine equations of motion 
for this proprotor model are simplified considerably if all the effects of 
sweep, dihedral, and angle of attack are neglected. Since these angles are 
always small for the model considered here, they generally contribute only 
small corrections to the coefficients. An exception, however, is the effective 
clastic axis shaft at the wing tip due to the sweep, which is an important 
aspect of the wing structural dynamics; this may be included by use of, for 
the mast height h y the distance from the hub to the effective elastic axis 
rather than to the actual wing tip spar. The simplified version of the equa- 
tions of motion then is obtained by neglecting all effects of , and & w 
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except that the effect of the effective elastic axis is included if h 7 , is 
used for h wherever it appears. In addition, the calculation of the rotor 
coefficients is simplified by considering only the 0 % terms (i.e., eqs. (54)); 

this primarily limits the results to below the critical tip Mach number and is, 
in fact, the form used almost exclusively for the rotor aerodynamics in the 
results presented here. Similarly, only the wing aerodynamic forces due to 
^L a are retained. This simplified model is not usually used here, and for the 
design and analysis of actual vehicles these simplifications of the structural 
and aerodynamic features would probably not be satisfactory. It will be shown, 
however, that the simplified model incorporates the essential features of the 
high inflow proprotor dynamics and so may be useful in further studies. 


Dropping the , and terms from the shaft motion expressions 

(eqs. (181)) yields 12 3 



and the wing equations of motion (eqs. (192)) reduce to 





+ 


0 

0 


0 





K* 

Ql 

0 


'- C pq y ?w <2y2 C T /aa) 
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The simplification is quite substantial. If one substitutes for the shaft 
motion into the rotor equations and forces, and then for the rotor forces into 
the wing equations, the matrix multiplications are easily carried out. The 
result is a set of nine equations of motion of the form: 


A 2 x + A\X + AqX = Bv (198) 


where the degree-of-freedom vector is 



and the control vector is 



The coefficient matrices are given below for the simplified equations. For 
convenience, in these matrices, q is written for Qqy' G/fy ) > U ^ ^ or 

h EA , and V for V + v; all C 7 and Cq terms are dropped and the structural 
damping of the blade is neglected; superscripts * on the inertias and spring 
constants are dropped. 
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Air Resonance 


The incorporation of the wing vertical bending motion into the model 
introduces the possibility of a mechanical instability involving that degree 
of freedom and the rotor blade inplane motion. Since instability involves the 
elastic bending of the airframe, it may occur in flight; hence it is termed 
air resonance. The basic mechanism involved, however, is identical with that 
in classical helicopter ground resonance (involving inplane hub motion, 
typically due to vibration of the helicopter on its landing gear while on the 
ground). The analysis and conclusions from ground resonance are thus directly 
applicable to air resonance also. 

Certain conclusions immediately follow then from the ground resonance 
analysis (ref. 34). Air resonance involves a coincidence of the wing vertical 
bending frequency (for the proprotor configuration) and the frequency (in the 
fixed system) of the lower lag mode ? - 1 . Any resonance with the upper lag 
mode (£ + 1) are always stable. Furthermore, the resonances with the lcwer 
lag mode are also stable if the rotor is stiff inplane, that is, if > 1/rev. 
Air resonance instability is possible only with a soft inplane rotor, that is, 
with v < 1/rev. 

Air resonance instability then occurs at a resonance of the frequencies 
and 1 - of the wing and rotor modes, with v p < 1/rev. Such conditions 

must be avoided, for example, by use of a wing that is stiff enough so that any 
resonance occurs at a rotor speed much higher than the normal operating 
rotational speed. Alternatively, it is possible to stabilize any resonances 
that occur by including sufficient damping - structural or aerodynamic - in 
the wing and rotor motions. 


High inflow operation of the rotor results in an aerodynamic damping of 
the lag motion, Q which is of order 1 in high inflow, compared to the order 
of the inflow squared in low inflow. Thus increasing forward speed in the air- 
plane configuration greatly increases the lag damping and eventually will 
stabilize the air resonance motion. Forward speed also contributes to stabil- 
ity by increasing the wing aerodynamic damping. Therefore, above a certain 
speed, no instability occurs even at a coincidence of the wing and lag fre- 
quencies; it is desirable, of course, that this speed be as low as possible, 
at least below the aircraft stall speed. 


A simple model will yield an estimate of the damping required, so of the 
forward speed required, to stabilize the air resonance motion. Consider the 
rotor lag and wing vertical bending degrees of freedom (C 1C -, and q ) 

retain only the direct damping of the wing and lag motion as the only aerody- 
namic influence. The homogeneous equation of motion, in Laplace form, reduces, 
for this set of degrees of freedom, to 


T c *(e 2 +C e+v 2 -l) 
-T/(2s + ( 7 c ) 

0 


I ? *(2s + C ? ) 

J^ ,t (s 2 +C^s+v^ 2 -l) 

S r*^T s2 
S w 



(199) 
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where the definitions of the coefficients follow by a comparison with the 
complete set (eqs. (198)). In particular, has been introduced for the 

total lag damping, M* C for the total win 8 vertical bending damping, and u>^ 

for the natural frequency (per rev) of the wing vertical bending motion. The 
characteristic equation is then 

I * 2 [(s 2 + C 8 + V 2 - l ) 2 + ( 2 s + C ) 2 ]M* (S 2 + c S + to* ) 

S b? ‘I II 

- I *(s z + C s + v 2 - l)S* 2 y 2 s 4 = 0 (200) 

? C c ^ 1 W 

The only coupling of the rotor and wing degrees of freedom is the inertia 
coupling that results from S'* 2 ; it follows then, as in ground resonance, that 
the damping required must also be of order 5* 2 . Examine the resonant case, 
where = 1 - v r (v < 1). Assume that the system is exactly on the stabil- 

ity boundary, so that s = too is a solution of the characteristic equation, 

with a) real. Let co 2 be an order 5* 2 distance from w 2 . Then, to the lowest 

2 C 
order in S * , the characteristic equation is 

V 22 W \ 2 * 1 - 


,2 2 4 

■v"c vv 


- 1 - « 0 


( 201 ) 


Since the resonant case is being considered, for to^ = 1 - v^, then 

- 1 • “?! = 2(V C - 11 = ‘ 2 "<7i 

v c 2 + 1 - = 2 \ 

and equation (201) becomes 

Vb w bV“?, 2 \ - = 0 

For dynamic stability, it is required then that 




S* 2 yl to 2 
2v 


( 202 ) 


that is, that the product of the wing and rotor damping be above a certain 
critical value. 
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To estimate the damping available, 
the wing: 


Mi C n = -yC* . 


consider the aerodynamic damping of 

i Vr; 

= y c T - 

lj a 5 T\aa 


and the aerodynamic lag damping: 


I *C = yQ‘ = y 

c e r s r 


The contributions of the structural damping of the wing and blade are negligible 
at the forward speeds involved. The boundary is usually at a low enough V that 
the small V approximation to Q • is accurate still; therefore, 

]/2 

I *C = yQ- S y — 

C ? S 4 

The requirement on forward speed V for air resonance stability is then: 


V 3 > 


S* 2 / 2v 


Y 2 (1/20)C7 <3 I mo 

^ W 


(203) 


which gives a critical V only a few percent below that for Q- = (V sin and 

so the low inflow approximation for Q £ is within the general accuracy of the 
result. This result is, of course, in dimensionless form, so the velocity is 
really the inflow ratio V/QR and the frequencies are all dimensionless values 
(per rev) . 


I’or the type of soft inplane rotor considered here, (per rev) does not 
vary much with ft, at least near and above the normal operating rotational 
speed. Consequently, the value of (per rev) at resonance with 1 - v is 

also independent of ft. Then the entire right-hand side of equation (203) for 
the ciitical 1 is a constant; it depends on the rotor (v^) and on some geo- 
metric properties of the rotor and wing. It does not, however, depend on ft or 
on (dimensional) because it is in dimensionless form. 


The criterion for stabilizing the air resonance motion is then a 
requirement for a fixed value of the ratio V/QR for a given rotor (soft 
inplane, with (per rev) independent of ft). The rotor rotational speed for 
the instability is given by the criterion of resonance of the q and ^ - 1 
modes; therefore, 1 


ft 


?l 


?! 


OJ/7 / 1 — V 

q \f c 


(204) 
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where in equation (204) ii»^ is dimensional and v^, per rev. Thus the rotor 

rotational speed (RPM) for resonance is directly proprotional to the wing 
vertical bending frequency (Hz). The forward speed required for stability at 
the resonance is given in terms of a critical inflow ratio: 


V = 


OR 


n 


critical 


(205) 


Thus the aircraft velocity (knots) required for stability is directly propor- 
tional to the rotor rotational speed (rpm) . Then increasing the wing vertical 
bending stiffness increases the ft at which air resonance occurs, but it also 
increases the V required for stability. Both effects are linear with w (Hz) 
i lo 

or for a fixed (per rev); on the V - ft plane, the locus of the 


highest air resonance instability speed is a constant V/QR line, that is, a 
straight line through the origin. The increase in ft is favorable for air 
resonance occurring above the normal operating range, but the corresponding 
increase in V is unfavorable. 


SECTION 5: RESULTS OF THE THEORY AND COMPARISON WITH FULL-SCALE TESTS 


Proprotor Dynamic Characteristics 


In this chapter (and the following two chapters), the dynamic 
characteristics of a proprotor operating at high inflow on the tip of a canti- 
lever wing are examined. The investigation uses the results of the theory 
developed previously. The cases examined are based on actual proprotor 
designs - two full-scale proprotors that were tested recently in the Ames 40- 
by 80-Foot Wind Tunnel on a dynamic test stand consisting of a cantilever wing 
simulating the full-scale aircraft wing stiffnesses. One rotor was designed 
and constructed by the Bell Helicopter Company, and the other by the Boeing 
Vertol Company. The Bell and Boeing rotors are shown in the configuration for 
the dynamic tests in the 40- by 80-Foot Wind Tunnel in figures 12 and 13, 
respectively. The two rotors differ primarily, so far as their dynamic char- 
acteristics are concerned, in the placement of the rotating natural frequencies 
of the blade flap and lag motions. The Bell rotor has a gimballed hub and 
stiff inplane cantilever blade attachment to the hub, hence vg = 1 (nearly, for 
it does have a weak hub spring) and > 1 ; it also incorporates positive 
pitch/flap coupling, Kp < 0 or 63 < 0, to increase the blade flap/lag stability. 
The Boeing rotor has a cantilever or hingeless hub with soft inplane blade 
attachment, hence vg > 1 and < 1 . The flap frequency is large even for a 
hingeless rotor because of the operation at lower rotor speed in airplane mode 
(the hover value of Vg is more typical of a hingeless rotor helicopter). The 
different placement of the blade frequencies, at opposing extremes of the 
range of choices, results in quite different dynamic characteristics for the 
two aircraft. A description of the rotors and the full-scale test results are 
given in references 25 to 28. For convenience, in the following discussion. 
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the gimballed, stiff-inplane rotor is referred to as the Bell rotor; the 
hingeless, soft-inplane rotor is referred to as the Boeing rotor. 

The dynamic characteristics of the rotors are examined in the following 
two chapters. Of primary interest are the eigenvalues - the frequency and 
damping of the roots of the coupled wing/rotor modes of motion. The basic 
parameters for this investigation are forward speed 7, rotor rotational speed 
ft, and inflow ratio 7/ ft/?. A velocity sweep, varying 7 (and so 7/ft/?) at con- 
stant ft is of interest since it is the way the rotor actually operates in air- 
plane mode. The rotor speed ft determines the relative values (i.e., 
dimensionless, per rev) of the wing natural frequencies, and so has a fundamen- 
tal influence on the dynamics. Varying ft also changes the blade natural fre- 
quencies, especially for the cantilever rotors considered here. A rotor speed 
sweep at constant 7 also varies the inflow ratio 7/12/?. The inflow ratio 7/ft/? 
is the primary parameter for the rotor aerodynamics and the wing aerodynamics 
as well. Varying 7/ft/? may also change the blade frequencies because of the 
change in the rotor collective pitch angle with the inflow ratio. 

Several elements in the theoretical model are examined to determine their 
influence on the proprotor dynamics: blade lag degrees of freedom, wing aero- 

dynamics, rotor speed perturbation degree of freedom, the complete expressions 
for the rotor aerodynamics, and a simplified theoretical model. The influence 
of the rotor lag motion can be examined by studying the effect of dropping the 
and C\ s degrees of freedom from the complete nine-degree-of- freedom model. 
The influence of the wing aerodynamics can be examined by setting all wing 
aerodynamic coefficients to zero and by comparing the results with those from 
the model that includes the wing aerodynamic forces. 

The rotor is usually considered to be operating in autorotation for the 
theoretical results presented here because this state is found to be least 
stable and also it is the configuration in which the full-scale dynamic tests 
were conducted. In autorotation, i 0 is the rotor speed perturbation degree of 
freedom, which is achieved by setting = 0 (as discussed previously); also, 

no rotor torque is delivered to the wing tip. The other extreme examined, 
considered the powered operation case, is when the huh operates at constant 
angular velocity (ft) with no perturbation; then c 0 is the collective lag 

degree of freedom, which is achieved by setting v f = v , and the rotor torque 

0 

perturbations are transmitted to the wing. This representation of the powered 
case is, of course, the limit of a perfect governor, but it provides an 
indication of the real powered state dynamics of the proprotor. 

The rotor aerodynamic coefficients are usually calculated for the present 
results using only the e 0 terms (eqs. (54)). This approximation simplifies 

the calculations considerably since then the coefficients depend only on V/QR, 
and it provides an adequate representation of the rotor aerodynamics for the 
range of inflow considered here, at least so long as the critical tip Mach 
number is not exceeded. The validity of the approximation can be assessed by 
comparisons with the dynamics that result when the complete expressions for 
the rotor aerodynamic coefficients are used (eqs. (50)). Since only a check 
on the influence of the complete expressions is desired here, representative 
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analytical expressions are used for the airfoil characteristics rather than 
the characteristics of the actual rotor blade sections (as described in eqs. 

(55) and (56)), To evaluate the complete coefficients, the angle-of-attack 
distribution over the blade span in the trim operating state is required, hence 
the rotor collective pitch is required. The determination of the collective 
pitch requires a performance analysis to find 0o.75 for a given Cj> or Cq; a 
simpler procedure - the use of an approximate collective pitch angle based on 
the inflow angle at 75 percent radius (plus about 1° to account for the mean 
operating section angle of attack) - is also evaluated. The evaluation of the 
complete rotor aerodynamic coefficients is also the only time the blade twist 
is required in the calculations (with the exception that when only the 

terms are used, the blade is assumed to have the twist required to maintain 
unstalled flow over the entire span). 

The effects of the wing sweep, dihedral, and angle of attack considerably 
complicate the coefficients of the equations of motion. Therefore, a simpler 
theory is evaluated in which all the 6^ , 6^ , and 5^ terms are dropped. The 

shift of the effective elastic axis at the wing tip is included by use of hgA 
(the distance from the hub to the effective elastic axis at the wing tip) for 
the mast height h (the physical shaft length in the unswept wing model). The 
simplified theory also uses only the o 0 terms for the rotor and wing aerody- 

namic coefficients. Of course, that is the usual approximation used here for 
the rotor coefficients, and it includes the primary effects of the wing coef- 
ficients. When the option of using the complete expressions for the rotor 
coef f icieint s is not included, the calculations are simplified considerably. 

The two full-scale rotors were also tested on a quarter stiffness wing, 
which had natural frequencies half those of the full stiffness wing. There- 
fore, by operating at half normal rotor rotational speed, the dimensionless 
wing frequencies (per rev) could be maintained and a given F/ft/? achieve^ with 
half the normal forward velocity F. With the quarter stiffness wing, it was 
possible to simulate operation at forward speeds twice the maximum capability 
of the wind tunnel, at least so far as the wing natural frequencies (maintained 
at the same per rev values) and the rotor aerodynamics (primarily a function 
of F/ft/?) were concerned. The rotor stiffness was, however, not correspondingly 
scaled down; therefore, the increase (per rev) in the blade flap and lag 
natural frequencies where the rotor was slowed down to half normal ft violated 
the simulation. For the gimballed rotor, at least the flapping mode (which 
has only centrifugal stiffening, except for the weak hub spring) was simulated 
on the quarter stiffness wing; but the lag frequency and both the flap and lag 
frequencies for the cantilever rotor were considerably different. Consequently, 
the results of the quarter stiffness wing tests are not presented here as 
simulating operation at twice the actual speed; the results are studied since 
they do provide additional data for investigating the dynamic characteristics 
of the proprotor. 

The theoretical studies concentrate on the eigenvalues as primary 
indicators of the dynamic behavior. However, the eigenvectors and also the 
frequency response of the system are examined at a typical operating condition 
(F/ft/? = 0.7). Finally, the results of the theory developed here are compared 
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with experimental data from the full-scale tests and also with results from 
other theories. 

With the nine-degree-of-freedom model of the present theory, there are 
eighteen roots or eigenvalues. The eigenvalues for each mode are usually a 
complex conjugate pair, which are then described by a frequency go and damping 
ratio e: 


a) = Im\ 

? = -Re\/ | X | 

The damping ratio C is the fraction of critical damping for the mode. The 
coincidence of the notation for the damping ratio and the rotor lag angle is 
unfortunate, but both are well established and the context should always 
clarify which is meant. 

Nine eigenvectors or modes correspond to the eigenvalues. The modes are 
identifiable by their frequency (which is near the uncoupled frequency, approx- 
imately the appropriate natural frequency as indicated in parentheses below) 
and by the participation of the degrees of freedom in the eigenvector. The 
modes are labelled here as follows: 


i 

<? 2 

P 

3 - 1 
e - i 

3 

3 + 1 
C + 1 

C 


wing vertical bending (gg^ ) 
wing chordwise bending (u> n ) 

h 2 

wing torsion (cjp) 
low-frequency flap (v - 1) 
low-frequency lag (v^ - 1) 
coning (v g ) 

high-frequency flap (vg + 1) 
high-frequency lag (v^ + 1) 
collective lag (v c ) 


Each mode involves motion of all nine degrees of freedom, of course (q^ 9 q 2 > 
p, V C 1 C* **1S’ and The low “ and hi ^ h - fre q u ency modes were 
discussed earlier. 


Description of the full-scale proprotors. - The parameters that describe 
the two full-scale proprotors are given in table 3 and in figures 14 to 17. 
The wing described in table 3 is for the configuration tested in the 40- by 
80-Foot Wind Tunnel. The rotors are flight-worthy designs. (For further 
descriptions of the full-scale rotors and aircraft, see references 25 to 28.) 
The parameters given in table 3 are those required for the present theory. 
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TABLK 3 . - DESCRIPTION OF THE BELL AND BOEING FULL-SCALE PROPROTORS 
AS TESTED IN THE AMES 40- BY 80-FOOT WIND TUNNEL 


otor 

Type 

Number of blades, N 
Radius, R 
Lock number, y 
Solidity, a 

Pitch/flap feedback, Kp 
Lift curve slope, a 
Rotor rotation direction, 

sgn fi 

Tip speed, OR 
Rotational speed, fi 


Blade frequencies (per rev) 


Gimballed, stiff 
inplane 
3 

3.82 m (12.5 ft) 
3.83 
0.089 
-0.268 
5.7 

+ 1 

183 m/ sec 
(600 ft/sec) 

458 RPM 
7.63 Hz 
48.0 rad/sec 

[l * 0.0355(^) 2 ] ' 

(fig. 16(a)) 
1.85 

figure 16(b) 


Boeing 


Cantilever, sofl 
inplane 
3 

3.97 m (13 ft) 
4.04 
0.115 
0 


160 m/sec 
(525 ft/sec) 

386 RPM 
6.43 Hz 
40.4 rad/sec 


figure 17(a) 


figure 17(b) 


Blade inertias 


0 for autorotation, for powered 


142 kg-m 2 
(105 slug-ft 2 ) 
1.000 


1.0 


.955 



TABLE 3.- DESCRIPTION OF THE BELL AND BOEING FULL-SCALE PROPROTOR 


AS TESTED IN THE MIES 40- BY 80- FOOT WIND FUNNEL 



Bell 

Blade inertias S^* 

1 .OSS 

% 

1.212 

Mi * 

6.160 

'r 


Except in autorotation 
where 

r* 

J* 

c o a 

% 


Blade structural damping 


Ps g 

0.1 percent 

g Ho 

9 S 

.5 percent 

.5 percent 

9,1 

. 5 percent 


1 

1 

0 


-Jig 

Semi span, yy 

V.llOrd , C*y Q 

Mast height, 

u 

rl 


Pylon center 

• rorti as 


n EA 

of gravity, zp m 


r* 

-"P, 


X 


1 1 




I* 

P 

L W 
C* 

pq 


1.333 

.413 

. 261 
.342 
.050 

76.9 

1.086 

1.206 

4.03 


.0141 

.667 


Continued 

Boeing 

1.092 
1 . 286 

4.344 


0.5 percent 
.5 percent 
.5 percent 
. 5 percent 


1 . 281 

.398 

.276 

.354 

.145 

110.93 

1 . 323 
2.926 

2.822 

. 009 

.667 


TABLE 3.- DESCRIPTION OF THE BELL AND BOETNC FULL-SCALE PROPROTOR 


AS TESTED IN THE AMES 40- BY 80-FOOT WIND TUNNEL - Continued 


Mode shape 


Bel 1 


W 

Du' G/r ) 

M 


0.S35 

1.74 


Sweep, 6^ 

Dihedral, 

Thickness ratio, 

Aspect ratio 

Aerodynamic center, x ^ /o w 

w 

Moment coefficient, e m 

m ac 


-6.5° 

0 

13. S percent 
6 . 6 
-.01 
-.005 


Boeing 


Bull stiffness wing 
Di mens ion a 1 stiffness 




h 

K, 


V 

Dimensional structural damping 


C, 




c r 


Dimensionless damping 


Dimensional ess stiffness 

K* 

K q 

‘2 
v * 

K P 

% 

C ? 2 
r * 

L P 

Typical resulting frequencies 

q 2 

P 


9.20x10 s kg m 2 /sec 2 
(6.793xi0 6 slug ft 2 / sec 2) 
2.50xl0 7 (1 . 840x1 0 7 ) 

1.77x10 s (1.305xl0 6 ) 


9030 kg m 2 /sec 

(6653 slug ft 2 /sec, t = 1 percent) 
27,300 (20,185, 1.8 percent) 

955 (703, 1.5 percent) 


18.72 

50.70 

3.595 

.880 

2.67 

.093 

3.2 Hz 0.42/rev 
5.35 .70 

9.95 1.30 


18.51 

50.14 

3.555 

.732 

2.22 

.077 

2.3 Hz 0.36/rev 
4.0 .62 

9.2 1.43 


TABLE 3.- DESCRIPTION OF THE 
AS TESTED IN THE AMES 40- 

BELL AND BOEING FULL-SCALE PROPROTOR 
BY 80-FOOT WIND TUNNEL - Concluded 


Bell 

Boeing 

Quarter stiffness wing 
Dimensional stiffness 



k q 

2. 28*10° kg m 2 /sec 2 


(1.687*10 6 slug ft 2 /sec 2 ) 

K a 

6.60x10 s (4.863x10 s ) 

K P 

5.30x10 s (3.908xl0 5 ) 

Dimensional damping 




5590 kg 

m 2 /sec 

(2646 slug ft 2 /sec, £ = 0.8 percent) 

% 

4700 (3462, 

0.6 percent) 

Cp 

441 (325, 1 

. 0 percent) 

Dimensionless stiffness 



K* 

4.65 

4.60 

<7l 



K iz 

13.40 

15.25 

K * 
P 

1.077 

1.065 

Dimensionless damping 

c* 

.35 

.29 

q i 



c* 

.458 

.581 

q 2 



Cp* 

.043 

.036 

Typical resulting frequencies 



q i 

q 2 
P ' 

1 . 5 Hz 0. 40/rev 

1.1 Hz 0.55/rev 

2.6 .68 

1.85 .58 

5.65 1.48 

5.55 1.73 

Full-scale aircraft wing 
Typical frequencies 



q 7 

0.45/rev 

0.49/rev 

.78 

.95 

£ 

p 

1.15 

1 

1.22 


142 



The Bell rotor is of gimballed, stiff inplane construction; the Boeing 
rotor is of cantilever, soft inplane construction. The Lock number y is based 
on the representative inertia X^, which is also given in table 3. The Bell 
rotor has -15° of 63, from which Kp = tan 63 = -0.268. The two rotors rotate 
in opposite directions; with sgn ft = 1, the Bell rotor rotates clockwise on 
the right wing, and with sgn ft = -1, the Boeing rotor rotates counterclockwise 
on the right wing. The blade twist and thickness distributions are given in 
figure 14 for the two rotors. The structural properties of the blades are 
given in figure 15. The rotor rotational speed ft and the corresponding tip 
speed ft/?, given in table 3, are the design operating speeds for airplane 
configuration. All theoretical results for frequencies are presented dimen- 
sionless, that is, as a fraction of ft (per rev); table 3 gives the dimensional 
values of ft for reference and to orientate the dimensionless results. 

The blade rotating natural frequencies are given in table 3 and in 
figures 16 and 17 for the Bell and Boeing rotors, respectively. Figure 16(a) 
shows the variation of the Bell flap frequency with rotor speed ft, from the 
expression given in table 3. Since the Bell hub is gimballed, vg is inde- 
pendent of 7/ft/? (collective pitch); the variation shown results from the hub 
spring restraint, which is rather weak at the normal ft. The coning mode 
natural frequency (for which the Bell rotor acts as a cantilever rotor) given 
in table 3 is only approximate, for it should vary with ft and with the inflow 
ratio (with collective) as the lag frequency does. The coning mode, especially 
with this high frequency, does not participate significantly in the proprotor 
and wing dynamics. The lag frequency for the Bell rotor (fig. 16(b)) is 
greater than 1/rev for the stiff inplane blade; the variation with 7/ft/? is 
actually the effect of the collective pitch variation. For normal ft, the Bell 
blade has a lag frequency from about 1.6 to 1.2 over the usual range of inflow 
ratio; for reference, the forward speed is about 250 knots at 7/ft/? =0.7 and 
normal ft. The flap and lag frequencies for the Boeing blade are shown in fig- 
ures 17(a) and 17(b) for this cantilever, soft inplane blade, Vg > 1/rev and 
< 1/rev (for the normal operating ft at least). The frequencies of the 
Boeing rotor do not vary much with inflow ratio (i.e., with collective) because 
the blade has a nearly isotropic shank construction to achieve the soft inplane 
lag frequencies; in fact, a soft inplane rotor for which the lag frequency did 
vary with collective would not be a useful design for the proprotor because of 
the restraints of blade loads and blade stability. For reference, the forward 
speed for the Boeing rotor is about 218 knots at 7/ft/? = 0.7 and normal ft. The 
Bell and Boeing rotors both have a hover rotational speed of about 550 rpm; 
the blade frequencies up to this ft value are given in figures 16 and 17. The 
blade frequencies at half normal operating ft (229 and 193 rpm, respectively, 
for the Bell and Boeing rotors) required for operating on the quarter stiffness 
wing, are also given in figures 16 and 17. 

The blade inertia constants are normalized by the Ifo values given. The 
blade structural damping values are only approximate; the structural damping 
of the blade is in any case negligible compared with the aerodynamic damping 
in high inflow. 

The same wing is used for both rotors. The wing semispan is 5.1 m 
(16.67 ft) from the tunnel floor to the rotor shaft; the difference in the 
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dimensionless values results from 
normalization. The mast height h 
and h^A is f rom the effective 
sweep (6^ < 0) of 6.5° which 


the use 
is from 


tip a 
shift 

tion 

only 


^3 

distance hpA " 
results in the 


rotor radii in the 
spar to the rotor hub, 
elastic axis to the hub. The wing has a forward 
shifts the effective elastic axis at the wing 


of different 
the wing tip 


h = 0.31 
value of 


m aft of the wing tip spar. The elastic axis 
g (yg ) given. The pylon center-of- gravity loca- 


w 


IS 

to 


the distance forward of the wing tip effective 


3 w 


mp * 


ip EA /h "r w - 


The value of 


VO: 


.) 


a structural dynamics analysis of the wing modes; note 
mode shape = y w 2 /yT w gives 


elastic axis; it is used 
given was obtained from 

that the approximate 


0JT U ) = 2. 


The inertia parameters of the wing were obtained from measurements; 
the wing stiffnesses were then obtained by finding the values required to 
match the experimental frequencies. This procedure makes use of the most 
accurately known wing properties (mass and frequencies) and ensures that the 
most important parameters of the wing dynamic characteristics - the 
frequencies - are accurately represented. The result was the dimensional 
stiffnesses given in table 3 for the full stiffness wing and the quarter 
stiffness wing; the dimensional stiffnesses are the same for both rotors since 
the same wing was used. The dimensionless stiffnesses are obtained by dividing 
by (N / 2) Ij-fl ?- . The values given are based on the normal operat ing rotor speed, 
and so must, in general, be multiplied by (tiRp/QR) 2 where ttR 0 = 183 and 
160 m/sec (600 and 525 ft/sec), respectively, for the Bell and Boeing rotors. 

To indicate the meaning of the stiffnesses, table 3 includes typical values 
of the resulting wing frequencies; since these values are for the coupled 
motion - at 100 knots and normal Q (one-half normal Q for the quarter stiff- 
ness wing) - they include the rotor inertia and aerodynamic influences, but 
the latter are small at this low speed. The results for the predicted and 
measured frequencies are always given here in dimensionless form (per rev); for 
reference, the corresponding dimensional values of the wing natural frequencies 
are given in table 3. 


The wing structural damping coefficients were determined from the 
stiffness and the measured values of the damping ratio (with the rotor off and 
no forward speed, so they include no aerodynamic influences) . The wing 
structural damping is not well known, which means there is some uncertainty in 
the predictions of the wing mode damping levels. The values given in table 3 
are based on the measurements with and without the rotor and wing aerodynamic 
influences, and are found to produce reasonable correlation with the data. The 
dimensionless damping values given are based on the normal rotor operating 
speed, so, in general, they must be multiplied by (QR o /0R) . 

The wing used, while having properties representative of the full-scale 
aircraft, actually was somewhat thinner {t w /c w = 13.5 percent as compared with 
the more likely 20-23 percent) and of quite different construction. Wing fre- 
quencies characteristic of the full-scale aircraft (for airplane cruise con- 
figuration and rotor speed) conclude table 3 to indicate the appropriateness 
of the wing used for the full-scale tests. 
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Other theories The results of the theory developed here compared with 
the results of other theories, as well as with the experimental data from the 
u 1-scale tests. The reports of the full-scale tests of the Bell and Boeing 
rotors (refs. 25 and 26) include predictions of the dynamic characteristics 
using theories developed by the companies for their rotors. 


The Bell Helicopter Company uses two theories, a linear model and a 
noniinear model. The Bell linear theory (refs. 15 and 25) consists of a 
C ose -form analysis for calculating eigenvalues and eigenvectors. The model 
as mear blade aerodynamics, but no wing aerodynamics. The degrees of free- 
om consist of flap and lag for each rotor blade and, for the wing and pylon 
there are five modes: vertical bending, chordwise bending, torsion, pylon 

yaw, and pylon pitch. The Bell nonlinear theory (refs. 25 and 29) consists of 
an open-form anaiysis for calculating the time history of the motion by numeri- 
in egration of the equations of motion. The model uses better blade aero- 
ynamics than does the linear theory, and incorporates wing aerodynamics 
(including rotor/wing interference effects). For this application at least 
the same degrees of freedom are used as for the linear theory. The nonlinear 
heory uses normal modes for the airframe (wing and pylon here) degrees of 
freedom, obtained from a NASTRAN calculations. ’ 

The Boeing Vertol Company uses a linear ^heory (refs. 26, 28, and 30) 
that consists of a closed-form analysis for calculating the~eigenvalues and 
eigenvectors. The degrees of freedom consist of flap and lag motion for each 
blade and normal modes for the wing. 


The Gimballed, Stiff- Inplane Rotor 


A 25-ft-diam flight-worthy gimballed, stiff-inplane proprotor, designed 
and constructed by the Bell Helicopter Company, was tested in the 40- by 8 80- 

iTol T7 e l T JUly / 97 °- 1116 configuration for the dynamics test^fig. 

~) consisted of the windmilling rotor mounted on the tip of a cantilever wing 
The rotor was operated in high inflow axial flow. The rotor and wing were 

rh^R 1 !? T reV1 ° USly - llie tCSt results * and also theoretical results from 
the Bell theories, are given in reference 25. Only the data for the case with 
the yaw link m are used here; that configuration had the pylon yaw stiffness 
at the airplane mode value. The theoretical dynamic characteristics of this 
or an wing are discussed, followed by a comparison with the full-scale 
test results and the Bell theoretical results. The aspects of the theory to 
be examined were discussed in the previous section. Y 

™ w | 11 ’ e t P rt dlCteti variation of the eigenvalues of the system with forward 
speed, at the normal airplane mode rotor speed (ft = 458 rpm) , is shown in 
figure 18: the frequency and damping ratio and the root locus. The wing mode 

frequencies show a slight decrease with V. The decrease in the lag frequency 
v with V/nE (really the effect of collective pitch) is apparent in the + i 
mode frequencies. The 0 ± 1 modes show the influence of the negative 6 3 , ^hich 
drops the effective flap frequency below 1/rev. For the rotor operating in 
autorotation, the rotor speed perturbation degree of freedom c is a first- 
order motion (a balance of the inertia and the aerodynamic lag damping)' so its 
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eigenvalue is on the negative real axis (not shown in fig. 18). The time 
constant of the C mode increased greatly from low inflow, where the lag amp- 
ing comes from the o la terms («*) . The damping of the wing chordwise bending 

mode <7, decreases with the damping of the torsion mode p increases. The 
damping of the wing vertical bending mode first increases with V, then it 
decreases. The peak occurs where the frequencies of the C - 1 and q x modes 
cross- hence it is an effect of the coupling of the rotor lag and wing bending 
modes. A similar effect is apparent in the q z damping where the C - 1 and q 2 
frequencies cross, but it occurs at very low inflow so the effect is smal . 

The q , mode becomes unstable at 495 knots (V/QR = 1.39); and the mode, at 
600 knots (V/QR = 1.69); these are the characteristic high inflow instabilities 
of the proprotor and wing configuration. Of course, the helic *J- **P_ M f* " um " 
ber ( M max = M t i p [ 1 + (k/O/?) 2 ] l/2 ) is unity at about 550 knots (V/QR - 1.55), 
certainly an upper limit on the validity of the theory. 

The influence of the rotor lag motion on the system stability is shown in 
figure 19 which compares the damping of the wing modes with and without the 
and Cl<7 degrees of freedom in the theory (g Q is the rotor speed perturba- 
tion for this autorotation case, so it must be retained); figure 19(c) shows 
the complete root locus without Z IC and ^5 in comparison with figure 18(c). 

The rotor lag motion has a very important influence on the wing modes, espe 
ciallv the q damping. Without the resonance of the lower-frequency lag mode 
with wing vertical bending, the damping does not show a peak at the resonance, 
rather it continues to increase with V until the high inflow effects appear at 
about 450 knots (V/QR = 1.27); then the damping drops off very quickly. 5o 
the rotor lag motion has the following influence on the stability: when the 

C - 1 frequency is greater than the q l frequency (roughly, when 
v > u + l/rev), the resonance increases the q l mode damping, while, *then 

t he g -1 frequency is less than the q l frequency, it decreases the damping. 

The result is that the rotor lag motion significantly reduces the wing vertical 
bending mode damping at the higher speeds. The speed at which the q 1 mode 
becomes unstable is not changed much, however, which indicates that the lg 
inflow instability mechanism is not greatly influenced by the lag motion. The 
reduction in damping at high speed is then more importantly accompanied by a 
great reduction in the rate at which the damping decreases, which is very 
beneficial The rotor lag motion thus makes the high inflow instability less 
severe The flap and lag modes in high inflow are highly coupled by the aero- 
dynamics. Eliminating the s ic and z is motions therefore greatly influences 
the behavior of the 6 ± 1 mode, as shown in the root loci of figures 18 and 19. 
The rotor flap motion is expected to be important to the proprotor and wing 
dynamics. The above comparison shows that the rotor lag motion can be equa y 
important. The influence of the lag motion is a combination of the high inflow 
aerodynamic forces and inertia coupling with the shaft motion. 


The influence of the rotor speed perturbation degree of freedom and the 
wing aerodynamics on the system stability is shown in figure 20. The basic 
case is autorotation operation, including the wing aerodynamics. Eliminating 
the wing aerodynamic forces decreases the q 1 and p mode damping, but has on y 
a small influence on the q 2 damping, which indicates that the C L q wing damping 

is the most important effect. Powered operation (including the wing 
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aerodynamics again) is considerably stabilizing for all the wing modes, over 
the autorotation case. The powered model considers the hub rotating at a 

TT? S ? e f d > so the ? mode becomes the elastic motion of the blades about 
the hub with spring restraint v The powered operation has little influence 
1 and C ± 1 modes, or on the wing mode frequencies. 

influence of the complete expressions for the rotor aerodynamic 
coefficients, that is, including the a t , a,, a 

^ U ( 


as the C£ 


terms, is shown in figure 21 


n , and c-j terms as well 

a x M 

The better rotor aerodynamic model 


reduces the predicted stability of the wing modes, for both autorotation and 
were cases. T e details of the coupling of the high-frequency rotor modes 
P + * and ? + 1 are also changed somewhat. The complete rotor aerodynamic 
coe lcients were calculated both by use of the correct collective pitch from 
a performance analysis (the col lective pitch required for Cn = 0 for autorota- 
ion or for the C T needed in equilibrium cruise for powered flight), and bv 
use of an approximate collective value based on the inflow at 75 percent radius 
f 0.75 ta u (K/M’)/(3/4) + 1.25°). The performance calculation is very sen- 
sitive to the collective pitch used, but figure 21 shows that the dynamics 
behavior is not; the approximate collective used is, in fact, within 1 or 2° 
o the correct value for both autorotation and powered flight at high speed 
and so evidently is an adequate representation of angle-of-attack distribution 
lhe complete expressions for the aerodynamic coefficients give somewhat dif- 
ferent numerical values compared with those obtained when only the terms 

are used (fig. 5), but the general behavior remains the same. An exception is 
when the drag divergence critical Mach number is exceeded. The helical tip 
Mach number exceeds the critical Mach number for the blade section character- 
istics used {M ovit = 0.9) at about 475 knots (V/tiR = 1.33); it exceeds the 
sonic value (M = 1) at about 550 knots (V/QR = 1.55). These points (fig. 21) 
are limits to the validity of the theory, but the main effects of the better 
blade aerodynamic model occur below these limits. It is concluded that using 
y the terms m the rotor aerodynamic coefficients is satisfactory for 

studying the basic behavior, and, in fact, is quite accurate so long as V is 
neither too small (low inflow) nor too large (stall and compressibility). For 
is example, the range in which the expressions are adequate is approxi- 

mately 7 = 25 to 350 or 400 knots (7/fiff = 0.1 to 1.0 or 1.1). When one pre- 
dicts the characteristics of an actual aircraft, however, especially the 
lgh-speed stability, the best available rotor blade aerodynamic model should 
be used, which probably means tabular data for the lift and drag coefficients 
as a function of angle of attack and Mach number for the blade sections used. 

Figure 22 shows the influence of using the simplified theoretical model 
on the predicted system stability. The effect is that of eliminating the wing 
sweep terms (except for the effective elastic axis shift, which is retained 
rou g EAj> the basic model already uses only the terms in the rotor 

aerodynamics and has no angle of attack or dihedral. The effect of the better 
blade aerodynamics was discussed previously; the effect of dihedral is 
expected to be similar to that for sweep; and there is little influence of 
angle of attack generally (either experimentally or theoretically, for small 
ang es at least). Therefore, the simpler theoretical model is quite 

A-5032 
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satisfactory for studying basic proprotor dynamics, giving the same general 
characteristics as the more involved model. For the design of an actual air- 
craft, however, a good structural analysis of the wing and pylon motion should 
be used. 

Figure 23 shows the behavior of the system dynamics during a rotor 
rotational speed sweep at 185 knots. The decrease in the wing mode frequencies 
is almost exactly proportional to £2 _1 , that is, the dimensional frequencies are 
nearly constant during the Q variation. The lag frequency decreases with £2 
faster than the (per rev) wing frequencies do. The C - 1 mode again shows a 
frequency resonance with the q ^ mode with increased damping when the £ - 1 
frequency is higher, and decreased damping when it is lower than the q-^ fre- 
quency. Some of the damping variation probably results from the high inflow 
influence. At low £1, a resonance of the 8+1 and p modes occur, which is 
apparent in both the frequency and damping of these two eigenvalues. 

The dynamic characteristics of the Bell rotor on the quarter stiffness 
wing, at half normal operating rotor speed (£1 = 229 rpm) , are shown in fig- 
ure 24, including a comparison with the full stiffness wing results (plotted 
vs. V/9R) . The frequencies of the modes are given in figure 24(a) (except for 
the 8, 6+1, and ? + 1 modes), and the great increase in the lag frequency 
that results from slowing the rotor is evident (see also fig. 16(b)). The 
wing frequencies are fairly well matched between the quarter- and full- 
stiffness wings. However, because of the difference in lag frequencies, the 
damping for the wing modes is not well simulated on the quarter-stiffness wing 
(figs. 24(b) and (c)), especially for the q l mode, which, for the full- 
stiffness wing, encounters a resonance with the t, - 1 mode. The influence of 
the rotor lag motion may be removed from the full-stiffness wing theory by 
eliminating the and ^5 degrees of freedom and, indeed, the q damping on 
the quarter-stiffness wing correlates well with that case. With the increased 
lag frequency on the quarter-stiffness wing, the p mode (instead of the q j 
mode) encounters a £ - 1 mode resonance, with a corresponding influence on the 
torsion damping. , 

Figure 25 shows the eigenvalues and eigenvectors for the Bell rotor at 
the typical cruise condition V/VH = 0.7, £2 = 458 rpm, V = 249 knots. This fig- 
ure is a time vector representation of the modes, so the eigenvector set for a 
given mode rotates counterclockwise at w = Im\ and decreases exponentially at 
a rate given by Re\. The projection of each vector on the real axis gives the 
participation of the degrees of freedom in the motion during the damped oscil- 
lation of the system in that mode. The degrees of freedom not shown for a 
given mode have a magnitude negligible on the scale used (i.e., less than 
about 5 percent of the maximum). The autorotation and powered cases show 
little difference except for the motion, of course, and in the wing mode 
eigenvalues. The rotor degrees of freedom participate significantly in the 
wing modes. The 6 ± 1 and C ± 1 modes show the coupling of the flap and lag 
motions due to the high inflow aerodynamics, but little coupling with the wing 
motion or with the collective rotor degrees of freedom. If B^p leads 6^5 in 
the time vector representation, the flap mode is progressive (the tip path 
plane wobbles in the same direction as the rotor rotation) and, similarly, for 
the lag modes. With the stiff-inplane rotor (v r > 1) and negative 63 (so the 
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effective vg < 1 ), then the g - 1 , g + 1 , and £ + 1 modes are progressive and 
the c, - 1 mode is regressive, as expected. 

The frequency response of the Bell rotor to each of the six input 
quantities is shown in figures 26 and 27 for autorotation and powered flight, 
respectively. The magnitude of the response of each degree of freedom to the 
input is shown; the rotor is operating at V/W = 0.7, ft = 458 rpm, and V = 249 
knots (the same as for the eigenvectors in fig. 25). The frequency response 
of the system is a good indicator of the dynamics involved, particularly the 
peaks in the response that occur at the resonant frequencies if the degree of 
freedom can be excited by that input. The frequencies of the eigenvalues are 
also shown (lower right) to identify the resonances. The wing vertical bending 
resonance (q ) is most important for the cyclic inputs (olq, g q, Q\q, and 0 ^), 
and the chordwise bending resonance ( q 2 ), for the collective inputs ( uq and 6 0 ). 
There are also significant resonances with the upper- frequency rotor modes 
(B + 1, C + 1). The degrees of freedom usually show significant excitation at 
the higher frequencies, especially near resonance with the wing modes, even if 
there is small or negligible steady-state response. The major differences 
between the powered and autorotation cases are the steady-state response 
(especially for the collective inputs), which carries into the low frequencies, 
too, and the response of the motion. 

The response shown at very low frequencies in figures 26 and 27 indicates 
the static response of the system to the six inputs. The system generally 
separates into a longitudinal or collective group (variables g Q and £ 0 and 
inputs Uq and 0 O ) and a lateral/vertical or cyclic group (variables $\q, B 15 , 
t>\C> and and inputs olq, g q, 0^, and 0jg). The wing variables (q^, q 2 , p) 
couple the two groups, but are excited most by the cyclic group. In autorota- 
tion, the static response of the cyclic rotor variables to the cyclic inputs 
is of order 1 for the flap motion and of order 0.1 for the lag motion; their 
response to the collective inputs is negligible. The static response of the 
collective rotor variables to the cyclic inputs is negligible; the response of 
3 0 bo 0 O is small, and to Uq, it is negligible. The response of Cq bo the 
collective inputs is of order 1 ; = of course, as discussed earlier 

(eq. ( 86 )). The static response of the wing variables to the collective inputs 
is negligible; the response of q , and p to the cyclic inputs is of order 0.05 
and the response of q 2 is of order 0.005. For powered flight, there is negli- 
gible effect on the response to the cyclic inputs compared to autorotation, 
but the response to collective inputs (0 Oj uq ) increases significantly. The 
static response of the cyclic flap motion to the collective inputs is then of 
order 0.05, the response of the cyclic lag motion is of order 0.005, the 
response of q j is of order 0 . 1 , the response of q 2 is of order 0.05, and the 
response of p is of order 0.01. The static response of the collective vari- 
ables (gg and £ 0 ) to the collective inputs in powered flight is of order 0 . 2 . 

Consider a comparison of the predicted dynamic characteristics for the 
Bell rotor with experimental results from the full-scale tests in the 40- by 
80-Foot Wind Tunnel and with the results of the Bell theories. Full-scale 
experimental data are available for the frequency and damping ratio of the 
wing modes. The data are limited by the tunnel maximum speed (about 200 knots) 
and by the use of an experimental technique that gave primarily only the damp- 
ing ratio for the wing vertical bending mode. The data were obtained by use 
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of an aerodynamic shaker vane on the wing tip (evident in fig. 12; the same 
technique was used for the Boeing rotor, fig. 13). The vane was oscillated to 
excite the wing motion desired; when sufficient amplitude was obtained, the 
vane was stopped and the system frequency and damping were determined from the 
subsequent decaying transient motion. Tills configuration is best suited for 
excitation of the wing vertical bending mode (<? ) . 

Figure 28 shows the variation of the system stability with velocity at 
the normal operating rotor speed (ft = 458 rpm) , in terms of the frequency and 
damping ratio for the wing modes. The results of the present theory are com- 
pared with the experimental data from the full-scale test, and with the 
results of the Bell linear and nonlinear theories. Reasonable correlation 
with both experiment and the Bell theories is shown. The good correlation of 
the frequencies predicted by use of the present theory with the experimental 
data (fig. 28(a)) follows because the wing stiffnesses were chosen specifically 
to match the measured frequencies (at around 100 knots). The difference 
between the predicted damping levels of the Bell linear and nonlinear theories 
is largely due to the neglect of the wing aerodynamic forces in the former. 

Figure 29 shows the variation with rotor speed ft of the wing vertical 
bending mode damping for the Bell rotor at V - 185, 162, and 150 knots. Rea- 
sonable correlation is shown with both the experiment and the Bell theories. 

For V - 162 and 150 knots, the predictions from the Bell theories are available 
only at normal operating rotor speed (ft = 458 rpm, from fig. 28(b)). 

Figure 50 shows the variation of the system stability with forward speed 
for the Bell rotor on the quarter-stiffness wing, at half normal operating 
rotor speed (ft = 229 rpm). During the full-scale test of this configuration, 
the available collective pitch was limited to the value reached at about 155 
knots (at ft = 229 rpm). Since the rotor was operated in autorotation, the 
collective pitch and inflow ratio V/QR were directly correlated. The maximum 
value of the inflow ratio was reached at 155 knots, where 7/ ft?? = 0.875. Above 
this speed, the collective was constant, and the inflow ratio was fixed at 
about 7/ ft?? = 0.840. The increase in velocity above 155 knots was accompanied 
by an increase in the rotor speed ft to keep the inflow ratio at the constant 
value demanded by the collective limit. The theoretical predictions include 
the actual rotor speed. The predicted frequency and damping with the rotor 
speed maintained at a constant value (ft = 229 rpm) are also shown in figure 50. 
The true values of the inflow ratio 7/ft/? for the experimental points above 
155 knots are shown in figure 50. Reasonable correlation is shown with both 
the experiment and the Bell theories. The decrease in the frequencies at high 
speed is produced mainly by the increasing ft. The increasing ft at high speed 
due to the collective limit significantly reduces the wing vertical bending 
and torsion damping, primarily because of the decrease in the effective (i.e., 
per rev) wing frequencies. 

The variation with rotor speed of the wing vertical bending and torsion 
damping for the Bell rotor on the quarter stiffness wing is shown in figure 51 
for 7 = 150 and 170 knots. Reasonable correlation is shown with experiment 
and the Bell theories. 
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Figure 32 shows the rotor flapping due to shaft angle of attack. The 
correlation of experiment and theory is shown in figure 32(a). The present 
theory predicts fairly well the magnitude of the flapping due to shaft angle 
of attack and also the longitudinal flapping However, the theory under- 

estimates the lateral flapping 315 by about a factor of 1/2. The results of 
the Bell linear theory are almost identical with the results of the present 
theory. The Bell nonlinear theory, however, predicts well the lateral flapping 
&15 also, as shown in figure 32(a). As discussed in reference 25, the better 
prediction of Bjp with the nonlinear theory is probably due to the inclusion of 
the influence of the wing-induced velocity on the rotor motion. Further evi- 
dence for that conclusion is the single point in figure 32(a) for which the 
present theory adequately predicts the lateral flapping 3i5- That point is 
from the powered test, which was not conducted on a wing. The lateral flapping 
3^5 is small compared to the longitudinal flapping Bi q, so the present theory 
does predict the magnitude of the flapping well. The azimuthal phase predic- 
tion has the same order of error as does 3i5, however. Figure 32(b) shows the 
predicted and experimental variation of the flapping with inflow ratio V/QR. 

The theoretical results are for a velocity sweep at normal rotor speed (ft = 

458 rpm) , while the experimental results include limited variation of ft as 
well as V, and the flagged points are even for the quarter-stiffness wing. 

Yet the flapping correlates well with the single parameter V/QR, indicating 
that the primary influence is the rotor aerodynamic forces. The underpredic- 
tion of fi 5 is again observed; the single point that agrees with the theory is 
the powered test point. 

Figure 33 shows the variation of the wing vertical bending ( q j) damping 
with 7/ft/?, during velocity sweeps for the Bell rotor on the ful 1 -stiffness and 
quarter-stiffness wings. The full-scale experimental data show a definite 
trend to higher damping levels with the full-stiffness wing than with the 
quarter-stiffness wing, and this trend correlates well with the present theory. 
The difference in damping at the same inflow ratio results from the lag 
motion. Figure 33(b) shows the frequencies of the £ - 1, q^ y and p modes for 
the full-stiffness and quarter-stiffness wings. The full -stiffness wing has a 
resonance of the £ - 1 and q modes that increases the q j damping below the 
resonance and decreases it above, and produces the peak m the damping 
observed in figure 33(a). Slowing the rotor on the quarter-stiffness wing 
greatly increases the lag frequency and removes it from resonance with q ^ 
(instead there is a resonance with the p mode, as shown in figure 33(b) and 
discussed earlier) . Another way to remove the influence of the rotor lag 
motion - in the theory - is to simply drop the and Z\s degrees of freedom 
from the full-stiffness wing case. Then the predicted wing vertical bending 
damping is almost identical to that for the quarter-stiffness wing (fig. 33(a)). 


The Hingeless, Soft-Inplane Rotor 


A 26-ft-diam, flight-worthy, hingeless, soft-inplane proprotor, designed 
and constructed by the Boeing Vertol Company, was tested in the 40- by 80-Foot 
Wind Tunnel in August 1972. The configuration for the dynamics test (fig. 13) 
consisted of the windmilling rotor mounted on the tip of a cantilever wing. 
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with the rotor operating in high inflow axial flow. The rotor and wing were 
described previously. The full-scale test data for the quarter-stiffness wing 
runs, and the theoretical results from the Boeing theory are from reference 26. 
The theoretical dynamic characteristics of this rotor and wing are discussed, 
followed by a comparison with the full-scale test results and the Boeing 
theoretical results . 

The major changes in the dynamic behavior compared with that of the Bell 
rotor are due to the different placement of the blade frequencies. The Boeing 
rotor has 1 - of the same order as the wing vertical bending frequency (as 
did the Bell rotor), but the soft-inplane rotor with < 1 introduces the 
possibility of an air resonance instability, that is, a mechanical instability 
that results from the resonance of the £ - 1 and modes. This instability 
will occur at a definite Q (for resonance) which, in this case, is above the 
normal operating rotor speed and at low forward speed. At high enough V/UR, 
the lag damping Q £ becomes large enough to stabilize the resonance. An analyt- 
ical discussion of the air resonance instability was given earlier. Besides 
introducing the possibility of an instability, at high SI and low V , the lag 
motion of the sofc-inplane rotor generally decreases the wing vertical bending 
mode stability. 

The Boeing rotor has cantilever blades with vg sufficiently above 1/rev so 
that vg - 1 is very close to the wing vertical bending mode frequency. Hence 
the 3 - 1 mode takes on many of the characteristics of the q mode, especially 
at high VftiR* In fact, it is usually the 3-1 mode that becomes unstable at 
high inflow rather than the q mode. By the time the 3-1 root enters the 
right half plane, the mode has however assumed the character of a wing vertical 
bending mode (this behavior is discussed further in terms of the eigenvectors 
of the two modes). Thus the high inflow instability mechanism is the same as 
observed already for the Bell rotor. 

The predicted variation of the eigenvalues of the system with forward 
speed, at the normal airplane mode rotor speed (Q ~ 386 rpm) , is shown in fig- 
ure 34: frequency and damping ratio and root locus. The flap frequency is 

greater than 1/rev and the coupled frequency of the 3-1 modes increases some- 
what with the inflow ratio. The lag frequency is less than 1/rev and the 
coupled frequency decreases with the inflow ratio. Since < 1, the £ - 1 
mode frequency increases. The proximity of the 3-1 and £ - 1 mode fre- 
quencies to the q^ y and even the q 2> frequencies is apparent in figure 34(a). 

The q damping is quite low at low speeds and has a minimum around 200 knots 
because of the influence of the £ - 1 mode, that is, the air resonance behav- 
ior. The qi damping increases at high 7, but there is considerable coupling 
of the 3-1 and q modes (as indicated by the frequencies and the eigenvec- 
tors). The 3-1 mode damping decreases very quickly at high speed and, by the 
time the root crosses into the right half plane at V = 480 knots (V/ttR = 1.54), 
the mode is really a wing vertical bending instability, that is, the high inflow 
inflow proprotor and wing instability. This change in the character of the 
3-1 and q^ modes is shown in figure 34(d), which presents the eigenvectors at 
V = 250, 400, and 500 knots. At low speed, the eigenvector on the left is 
clearly identifiable as the q^ mode, and the eigenvector on the right as the 
3-1 mode, based both on the frequency of the root and on the participation of 
the degrees of freedom in the eigenvector. As forward speed increases, the 
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wing vertical bending motion decreases in one mode and increases in the other. 
The mode that is originally the rotor low-frequency flap mode (B - 1) becomes 
unstable just before 500 knots, and by that time this mode has assumed the 
character of the primary wing vertical bending mode. Note that the wing 
vertical bending motion is characterized not simply by the q 1 degree of free- 
dom, but also by the motion of ^ £ 0 , and p associated with the mode. 

The wing chord ( q 2 ) mode damping decreases with speed until an instability 
is encountered at V = 510 knots (7/121? = 1.64). This is an air resonance 
instability, as indicated by the coincidence of the ? - 1 and q 2 mode fre- 
quencies at this speed (fig. 34(a)). Wing chord bending produces a lateral 
motion at the rotor hub forward of the wing tip, which couples with the rotor 
lag motion. An air resonance instability can occur at even high speed with 
the wing chord mode because the wing aerodynamic damping of that mode remains 
small. The q 2 instability occurs at a slightly higher speed than the B - l/q 1 
instability, so, in some cases, it may be the critical boundary. 

The wing torsion (p) mode couples with the rotor coning (B) mode in this 
case, mainly due to simply a coincidence of the damping and frequencies of 
the two modes. These modes have fundamentally different character (Bq is in 
the longitudinal group of variables and p, in the lateral/vertical group) and 
do not really want to couple. The roots try to cross on the root locus plane 
(fig. 34(c)) and instead exchange roles; the coupling is significant only in 
a narrow region near 300 knots; elsewhere, the roots are clearly distinguish- 
able. While this coupling does not have great physical significance, it is 
discussed because a slight change in the parameters or in the model may elimi- 
nate the coupling. For comparison with such cases, it is most convenient to 
plot the damping (fig. 34(b)) as if the root loci really did cross, that is, 
by joining the corresponding p and B pieces. This practice is followed in the 
comparisons that follow. 

The influence of the rotor lag motion on the system stability is shown in 
figure 35, which compares the damping of the wing modes with and without the 
Cl C anc * Ci s degrees of freedom in the theory. The rotor lag motion has a 
large and important influence on the wing modes. The rotor lag motion sub- 
stantially decreases the stability of the q j and B - 1 modes. The B - 1 mode 
behavior remains the same when the lag degrees of freedom are eliminated, but 
the sharp damping decrease (and instability) occurs at a speed about 250 to 
300 knots higher, beyond the scale used in figure 35(a). The low damping of 
the mode around 200 knots is shown to be due to coupling with the £ - 1 
mode, that is, air resonance behavior. The rotor lag motion also decreases the 
q mode stability at high 7, another air resonance effect. The lag motion 
stabilizes the p mode, but that is not really needed. The complete root locus 
is shown in figure 35(c), which is to be compared with the root locus that 
includes the lag motion (fig. 34(c)). 

The influence of the rotor speed perturbation degree of freedom and the 
wing aerodynamics on the system stability is shown in figure 36. The basic 
case is autorotation operation, including the wing aerodynamics. The rotor 
speed perturbation degree of freedom generally decreases the stability, that 
is, powered operation is more stable, especially for the mode where the air 
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resonance behavior is much less noticeable. The wing aerodynamic forces 
generally increase the stability. 

The influence of the complete expressions for the rotor aerodynamic 
coefficients is shown in figure 37(a) (for clarity, only the autorotation case 
is shown for the 6 - 1 and p modes). Generally, the use of the better blade 
section aerodynamics decreases the predicted stability of the wing modes. The 
correct collective was obtained from a performance analysis for powered and 
autorotation operation of the Boeing rotor. The approximate collective used, 
based on the inflow angle at 75 percent radius, was 00-75 = tan " 1 (7/12/?)/ (3/4) + 
1.0°. The helical tip Mach number reaches the blade critical Mach number 
(M = 0.9) at about 500 knots (7/12/? = 1.61) and reaches the sonic value at 
about 580 knots ( T //f IR = 1.86). The conclusions are the same as for the Bell 
rotor: the basic behavior of the system is described well with only the a £ 

terms in the rotor aerodynamic coefficients, but the complete expressions 
should be used to obtain correct predictions for actual vehicles, particularly 
for the high-speed stability boundaries. The use of the approximate collective 
does not influence the dynamics much, although it is, of course, not 
satisfactory for performance calculations. 

Figure 38 shows the influence of the use of the simplified theoretical 
model on the predicted system stability. As for the Bell rotor, it is con- 
cluded that the simpler model is satisfactory for studying the basic behavior, 
but for the design of a particular vehicle, the best available model should be 
used. 


Figure 39 shows the variation of the system eigenvalues with rotor speed 
12 for the Boeing rotor at 50 knots. At this low speed, the C - 1 and q ^ fre- 
quency resonance around 530 rpm (fig- 39(a)) results in an air resonance insta- 
bility in the < 7 . mode (fig. 39(b)). At resonance, there is a corresponding 
increase in the 5 - 1 damping. The resonance and corresponding instability 
occur above the normal rotor operating speed (12 = 386 rpm) even with the wing 
used for the wind-tunnel test, which was softer in bending than the full-scale 
design. The general decrease in the 5 ± 1 mode damping with 12 results from 
the low lag damping at low inflow. Figure 40 shows the variation of the 
eigenvalues with 12 for the Boeing rotor at 192 knots. The t, - 1 and q j 
resonance again occurs at about 12 = 500 rpm, but this speed is sufficient to 
stabilize the air resonance motion. Figure 41 summarizes the Boeing rotor air 
resonance behavior: the wing vertical bending mode (< 7 ^) damping variation with 

rotor speed 12 for V = 50 to 192 knots. The stabilizing influence of the for- 
ward speed is shown. An earlier section derived an estimate for the V value 
required to stabilize the air resonance motion. In this case, resonance occurs 
with the wing vertical bending frequency of about 0.28/rev, and of about 
0.8/rev; with the other parameters required from table 3, equation (203) 
gives 7/127? > 0.268 for the stability requirement. At this speed, resonance 
occurs at 12 = 500 rpm, so the velocity requirement is F > 108 knots. The use 
of the equivalent radius approximation for the rotor lag damping Q £ gives 
instead 7/12 R > 0.285 or V > 114 knots, which is only about 6 percent higher. 

The estimate compares well with the calculated boundary of about 120 knots 
(fig. 41), better, in fact, than is reasonable to expect from the simple model 
used for the air resonance estimate. 


154 



The dynamic characteristics of the Boeing rotor on the quarter-stiffness 
wing, at half normal operating rotor speed (Q = 193 rpm) are shown in figure 42. 
The frequencies of some of the corresponding modes on the ful 1 -stiffness wing 
are also shown in figure 42(a) (plotted vs. V/fiR; the 6+1, £ + 1, and 6 fre- 
quencies are not shown for the full-stiffness wing). The wing mode frequencies 
are well matched to the full stiffness wing values (per rev), but slowing the 
rotor increases both the flap and lag frequencies of the blade considerably, 
the flap frequency to near 2/rev and the lag frequency to near 1/rev, as com- 
pared with about Vg = 1.35 and = 0.75 at normal Q (see also fig. 17). The 
lag frequency moves nearer 1/rev and thus the c - 1 mode frequency is lower for 
the quarter-stiffness wing (besides the influence on the dynamics, the lag fre- 
quency near 1/rev also means large vibration and blade loads). With the rotor 
frequencies so different, the system damping shown in figure 42(b) has much 
different behavior than for the full-stiffness wing (compare with fig. 34(b)), 
especially for the q. and 6-1 modes. Figure 43 shows the variation of the 
system stability with rotor speed ft, for the Boeing rotor on the quarter- 
stiffness wing at 80 knots. Air resonance effects are evident in q at about 
400 rpm and in q^> at about 500 rpm. The peak in the q damping at 225 rpm is 
due to the coupling with the 6-1 mode. 

Figure 44 shows the eigenvalues and eigenvectors for the Boeing rotor at 
the typical cruise condition of V/QR - 0.7, Q, = 386 rpm, V = 218 knots. With 
this soft-inplane (v^ < 1) and cantilever (vg > 1) rotor, the r, - 1, 6 + 1, 
and £ + 1 modes are progressive and the 6-1 mode is regressive as expected. 

The frequency response of the Boeing rotor to each of the six input 
quantities is shown in figures 45 and 46 for autorotation and powered flights, 
respectively. The magnitude of the response of each degree of freedom to the 
input is shown; the rotor is operating at V/QR - 0.7, Q = 386 rpm, and V = 218 
knots (the same as for the eigenvectors in fig. 44). The steady-state (low- 
frequency) response, compared with that of the Bell rotor, shows only the fol- 
lowing differences: with the hub moment capability of the cantilever rotor 
(vg < 1), the flap motion with respect to the shaft is reduced, and the wing 
motion is increased. There is increased lag motion because of the softer 
blade inplane restraint (lower v^) and there is a change in the azimuthal 
phasing of the cyclic rotor response (e.g., $\c and Big) to the cyclic inputs 
(e.g., 0ig and 0^g) because of the change in rotor frequencies. 

Consider a comparison of the predicted dynamic characteristics for the 
Boeing rotor with experimental results from the full-scale tests in the 40- 
by 80-Foot Wind Tunnel and with the results of the Boeing theory. Full-scale 
experimental data are available for the frequency and damping of the wing 
vertical bending mode, obtained by the same shaker vane excitation technique 
as used with the Bell rotor. 

Figure 47 shows the variation of the system stability with velocity at 
the normal operating rotor speed (ft = 386 rpm) in terms of the frequency and 
damping ratio for the wing modes. Reasonable correlation of the present theory 
with both experiment and the Boeing theory is shown. However, data are 
available only for wing vertical bending mode damping. 
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Figure 48 shows the variation of the wing vertical bending mode damping 
for the Boeing rotor with rotor speed ft at 7 = 50 to 192 knots. These runs 
were conducted to investigate the air resonance behavior of this proprotor and 
wing configuration. Reasonable correlation is shown with both experiment and 
the Boeing theory, except at the higher tunnel speeds. There the data show 
considerable scatter because the tunnel turbulence made analysis of the 
transient motion difficult. 

Figure 49 shows the variation of the system stability with forward speed 
for the Boeing rotor on the quarter-stiffness wing, at half normal operating 
rotor speed (ft = 195 rpm) . Figure 50 shows the variation with rotor speed ft 
for the Boeing rotor on the quarter-stiffness wing at V = 80 knots. Reasonable 
correlation is shown with both experiment and the Boeing theory. The ft sweep 
also show's the air resonance behavior in both theory and experimental data. 


SECTION 6 : COMPARISONS WITH OTHER INVESTIGATIONS 


In this chapter, the present theory and the results obtained are 
compared with the published work of other authors; of primary interest here 
are the theoretical models developed in the literature. 

Hall (ref. 8 ) discussed the role of the negative H force damping on the 
high inflow proprotor behavior, reviewed the problems found in the XV-3 flight 
tests, and reviewed the results of the 1962 test of the XV-3 in the 40- by 80- 
Foot Wind Tunnel. He presented an investigation of the influence of various 
parameters on the stability of the rotor and pylon, particularly forward speed, 
pylon pitch and yaw spring rate, and pitch/flap coupling ( 63 ); this investiga- 
tion used the full-scale XV-3 test results, model tests that simulated the 
XV-3 configuration, and analysis results from a theory presented in the paper. 
Hall derived the equations of motion for a two-bladed rotor on a pylon; N = 2 
was chosen because the analysis was to support the XV-3 investigation. His 
model then had three degrees of freedom: flap angle B (teetering), pylon 

pitch <p x , and pylon yaw <py. The present results (eqs. (146) for N = 2) agree 
with Hall's equations, with the following correspondence of notation: 

Present Notation Hall 

~ a y 

a <f> 

M F (1/2)/% 

H H sin 4> 

Y H cos ip 

Hall considered only the case of Vg = 1, that is, no hub spring restraint of 
the teetering blade; therefore, no hub moment due to flapping is transmitted 
to the pylon motion in his model. The aerodynamic forces Mg and H were 
expressed in terms of integrals of the blade section forces F z and F x over the 
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span, which agree with the present results except that the radial drag force 
was neglected. Hall did not, however, expand the rotor aerodynamic forces in 
terms of the perturbed rotor and pylon motion because he did not derive a set 
of linear differential equations. Hall solved for the dynamic behavior by 
numerically integrating the equations of motion; hence he found the transient 
motion rather than eigenvalues because of the periodic coefficients for N = 2. 
With this method, the exact, nonlinear aerodynamic forces could be included 
rather than the linearized expansion. 

Gaffey, Yen, and Kvaternik (ref. 11) discussed the proprotor aircraft 
behavior and design considerations in relation to the wing frequencies, gust 
response, and ride quality. The influence of the blade frequencies and pitch/ 
flap coupling on the rotor and the rotor/wing stability were discussed. It 
was shown that a cantilever rotor, that is, Vg > 1 , has greater stability, and 
that vg > 1 reduces flapping significantly but also increases blade loads. 
Expressions were given for the low-frequency response of flapping to shaft 
angle of attack ( Xp/V here) and shaft angular velocity (a^ here) in terms of 
the equivalent radius approximations; the present results (eqs. (96) to (100)) 
agree with their expressions. Experimental and theoretical data were given for 
proprotor/wing stability, flapping, loads, vibration, and gust response. 

Tiller and Nicholson (ref. 13) discussed the stability and control 
considerations involved with proprotor aircraft. They found the following 
influences on the aircraft stability. The proprotors with positive pitch/flap 
coupling and clockwise rotation on the right wing produce, through the negative 
H force, an increased effective dihedral in the effect increasing with 

forward speed. The proprotor negative damping requires a larger horizontal 
tail for the short-period mode frequency and damping; the rotor contribution 
found was on the order of 30 to 40 percent of the stabilizer contribution. 
Similar results were found for the vertical tail requirements (G^g) • The rotor 

thrust damping in yaw (2^) contributes significantly to C nr . The rotor torque 
and hub force during rolling increase by about 30 to 50 percent of 

the wing value. The thrust due to rolling (2^) produces adverse yaw {hC n ^ < 0, 

for this rotor rotational direction) that appreciably alters the Dutch roll 
damping and mode shape. The rotor influence on the aircraft stability deriva- 
tives found here agrees with the results of Tiller and Nicholson. They also 
discussed other features of the proprotor configuration that influence the 
aircraft stability and control: the thick wing, the high roll inertia, the 

nacelle contribution to , the important influence of the interconnect shaft 

stiffness on the lateral derivatives (particularly, C n ^ and C n ) , and the 

control features in helicopter and cruise mode. They point out that the 

influence of the rotors on the lateral dynamics is more complex than that on 

the longitudinal dynamics, but that meeting the requirements is largely a 

matter of enough vertical tail effectiveness. 

Young and Lytwyn (ref. 18) developed a four-degree-of-freedom theoretical 
model (8i£, 615 , a y, and a x ) for studying proprotor dynamics. They found an 
optimum value for the flap stiffness for pylon/rotor stability at about 
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Vg = 1.1. An approximation to this result was obtained by setting to zero the 
term that couples the rotor with the pylon; that is, in the present notation. 


ycv 


2 

3 


y 


i) 



o 


There is then no moment about the pivot due to tip path plane tilt, which 
greatly increases the rotor/pylon stability. This optimum was discussed in a 
previous section and was also the subject of the discussion of reference 18 by 
Wernicke and Gaffey. Young and Lytwyn presented several results for the whirl 
flutter case (a truly rigid propeller on the pylon), which were also discussed 
previously. Young and Lytwyn found the power-on case to be less stable than 
the windmilling case; they were considering, however, the case of a t 0, that 
is, the influence of the a £ terms in the rotor aerodynamic coefficients. The 
present results confirm that the use of the better calculation of the aerody- 
namic coefficients decreases the predicted stability. The really important" 
factor in windmill operation (autorotation) is the rotor speed perturbation 
degree of freedom, which makes the windmilling case much less stable than the 
power-on case. The theoretical model considered by Young and Lytwyn was an 
a-bladed rotor (.V 1 3) on an elastically restrained pylon with pitch and yaw 
degrees of freedom. The blade motion allowed was rigid flapping (ng = r), but 
elastic blade restraint was included so that Vg > 1 was possible. Only the 
lotor tip path plane tilt couples with the pylon motion, so the system reduces 
to four degrees of freedom. The same four-degree-of- freedom model was 
considered here (fig. 1), with the corresponding notation: 

Present Notation Young and Lytwyn 

&C 
8,3 


Although the derivation of the rotor aerodynamic coefficients in the linear 
equations was considerably different from that used here, the final form is 
essentially equivalent to the present result. The corresponding notation for 
the rotor aerodynamic coefficients is 


1 /nr 

16 

io 


-ou 
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Young and Lytwyn evaluate these coefficients assuming constant a and a d over 
the blade span. The set of four equations of motion obtained correspond to 
equations (26) . 

Descriptions of analyses typical of the most sophisticated currently used 
for calculating the dynamic characteristics of tilting proprotor aircraft may 
be found in references 29 and 30. These are, in fact, the most complete 
descriptions available in the literature for the proprotor aircraft analyses. 
Their use lies primarily in the development and support of the design of 
specific aircraft. More elementary models remain valuable for general and 
exploratory investigations of proprotor dynamics. 

Descriptions of the tilting proprotor aircraft, and the design 
considerations involved, may be found in references 10, 15, 22, 27, and 28 . ^ 

The XV- 3 flight test results are described in references 2 and 3, and the XV-3 
tests in the 40- by 80-Foot Wind Tunnel are described in references 1, 8, and 
9. Recent tests of full-scale proprotors in the 40- by 80-Foot Wind Tunnel 
are described in references 14, 15, 25, and 26. Some experimental data from 
small-scale model tests are also available (refs. 11, 15, 16, and 24, for 
example) . 


CONCLUDING REMARKS 


A theoretical model has been developed for a proprotor on a cantilever 
wing, operating in high inflow axial flight, for use in investigations of the 
dynamic characteristics of tilting proprotor aircraft in the cruise configura- 
tion. The equations of motion and hub forces of the rotor were found includ- 
ing the response to general shaft motion. This rotor model was combined with 
the equations of motion for a cantilever wing. In further studies, however, 
the rotor model could easily be combined with a more general vehicle or sup- 
port model, including, for example, the rigid-body degrees of freedom of the 
aircraft. The general behavior of the high inflow rotor has been investigated 
and, in particular, the stability of the proprotor and cantilever wing con- 
figuration. The effects of various elements of the theoretical model were 
examined, and the predictions were compared with experimental data from wind- 
tunnel tests of two full-scale proprotors. 

From the theoretical results for the two full-scale rotors, and 
comparisons with the full-scale, wind-tunnel test data, it is concluded that 
the nine-degree-of-freedom model developed here is a satisfactory representa- 
tion of the fundamental proprotor dynamic behavior. The model consists of 
first mode flap and lag blade motions of a rotor with three or more blades, 
and the lowest frequency wing modes. The limitations of the present theory 
are primarily the structural dynamics models of the rotor blades and the wing 
and the neglected degrees of freedom of the proprotor aircraft system. For 
the rotor, it was assumed that the blade flap and lag motions are not coupled, 
that is, are pure out-of-plane and pure inplane motions, respectively. The 
model neglected the higher bending modes of the blades, and the blade elastic 
torsion degrees of freedom were neglected entirely. For the support, the 
model used only an elementary representation of the structural modes of the 
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wing. The model was limited to the cantilever wing configuration, 
neglecting the aircraft rigid-body degrees of freedom as well as the higher 
frequency modes of the wing and pylon. The present model does incorporate the 
fundamental features of the proprotor aeroelastic system. Hence these limita- 
tions of the model are primarily areas where future work would be profitable, 
rather than restirctions on its current use. 

From a comparison of the behavior of the gimballed, stiff-inplane rotor 
and the hingeless, soft-inplane rotor, it is concluded that the placement of 
the rotor blade natural frequencies of first mode bending - the flap frequency 
Vg and the lag frequency - greatly influences the dynamics of the proprotor 
and wing. Moreover, the rotor lag degrees of freedom was found to have a very 
important role in the proprotor dynamics, for both the soft-inplane 
( V j. < 1/rev) and the stiff-inplane (v > 1/rev) configurations. 

The theoretical model developed here has been established as an adequate 
representation of the basic proprotor and wing dynamics. It will then be a 
useful tool for further studies of the dynamics of tilting proprotor aircraft, 
including more sophisticated topics such as the design of automatic stability’ 
and control systems for the vehicle. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Dec. 26, 1973 
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Velocities seen by rotor blade section and resulting aerodynamic 
forces on the section. 
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Figure 5. Influence of blade section aerodynamics on rotor coefficients for 
the Bell and Boeing rotors, compared with coefficients based on only the 
c „ terms. J 
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Figure 6.- Whirl flutter (two-degree-of- freedom) divergence instability 

boundary. 
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Figure 8.- Typical whirl flutter stability boundaries. 
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Figure 9.- Whirl flutter stability boundaries for a two-bladed rotor. 
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Figure 10.- Axis system and geometry for evaluating proprotor influence on 

aircraft stability derivatives. 





Figure 12.- Bell Helicopter Company full-scale, 25-ft-diam proprotor on 
cantilever wing for dynamic tests in the Ames 40- by 80-Foot Wind 
Tunnel . 
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Figure 13.- Boeing Vertol Company full-scale, 26-ft-diam proprotor on 
cantilever wing for dynamic tests in the Ames 40- by 80-Foot Wind 
Tunnel. 
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(b) Section flapwise modulus/moment product. 

(c) Section chordwise modulus/moment product. 

Figure 15.- Structural characteristics of two full-scale proprotor blades 
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Figure 17.- Blade rotating natural frequencies for Boeing rotor. 
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(a) Frequency of modes. 
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Figure 18.- Bell rotor velocity sweep, Cl = 458 rpm; predicted eigenvalues 
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Figure 18.- Concluded. 
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(c) Root locus. 
Figure 19.- Concluded. 
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C a j Wing vertical bending damping. 

(b) Wing chordwise bending {q^) and torsion (p) damping. 

Figure 20.- Bell rotor velocity sweep, £2 = 458 rpm; comparison of basic 
(autorotation and wing aerodynamics), powered, and no wing 
aerodynamics cases. 
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21.- Bell rotor velocity sweep, ft = 458 rpm; influence of the complete 
expressions for the rotor aerodynamic coefficients. 
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(a) Frequency of the modes. 

(b) Damping ratio of the modes. 

Figure 23.- Bell rotor rpm sweep, V = 185 knots 
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(a) Frequency of the modes. 

(b) Wing vertical bending damping. 

(c) Wing chordwise bending (q 2 ) and torsion (p) damping. 

Figure 24.- Bell rotor velocity sweep, quarter-stiffness wing, 5? = 
compared with the full-stiffness wing results with and without £ 
r degrees of freedom. 
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Figure 25.- Concluded. 










(a) Vertical gust a^, input. 

Figure 27.- Bell rotor in powered operation at V/Q.R =0.7 = 458 rpm 

V = 249 knots), magnitude of response of each degree of freedom to 
input at frequency w. 
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Figure 27.- Continued. 




(d) Lateral cyclic pitch Q 1C input . 
Figure 27.- Continued. 
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Figure 27.- Continued. 
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(a) Frequency of the inodes. 

(b) Wing vertical bending damping (<?,). 

(c) Wing chordwise bending damping (<? 2 ) . 

(d) Wing torsion damping (p). 

Figure 28.- Bell rotor velocity sweep at SI = 458 rpm, comparison with full 
scale experimental data and the Bell theories. 
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(a) Wing vertical bending damping C<7 j D at 185 knots. 

(b) Wing vertical bending damping at 162 knots. 

(c) Wing vertical bending damping (c^) at 150 knots. 

Figure 29.- Bell rotor rpm sweep, comparison with full-scale experimental data 

and Bell theories. 
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Figure 32.- 


Bel 1 rotor flapping due to shaft angle of attack; 
full-scale experimental data and Bell theories 


comparison with 
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(a) Frequency of the modes. 

(b) Damping ratio of the modes. 

Figure 34.- Boeing rotor velocity sweep, 2 = 386 rpn; predicted eigenvalues 
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(a) Frequency of the modes. 

(b) Damping ratio of the modes. 

Figure 34.- Boeing rotor velocity sweep, fi = 386 rpm; predicted eigenvalues 
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(d) Boeing rotor at ft = 386 rpm; variation with velocity of predicted 
eigenvectors for wing vertical bending ( q ) and rotor flap (g - 1) 
modes . 1 


Figure 34.- Concluded. 
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(a) Wing vertical bending (^ 2 ) and rotor flap (3 - 1) damping. 

(b) Wing chordwise bending (<? 2 ) and torsion (p) damping. 

Figure 36.- Boeing rotor velocity sweep, ft = 386 rpm; comparison of basic 
(autorotation and wing aerodynamics), powered, and no wing aerodynamics 
cases . 
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(a) Wing vertical bending (q^) and rotor flap (8 - 1) damping. 

(b) Wing chordwise bending (<? 2 ) and torsion (p) damping. 

Figure 37.- Boeing rotor velocity sweep, = 386 rpm; influence of the complete 
expressions for the rotor aerodynamic coefficients. 
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Simplified theory model 









200 400 600 800 

iiR , fps 


(a) Frequency of the inodes. 

(b) Damping ratio of the modes. 

Figure 40.- Boeing rotor rpm sweep; V = 192 knots. 
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Figure 41.- Boeing rotor air resonance behavior; wing vertical bending damping variation with rotor 

speed at V = 50 to 192 knots. 
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Figure 44.- Boeing rotor roots and modes at V/tiR = 0.7, Q = 386 rpm, V = 218 knots. 
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(b) Powered flight. 
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(a) Vertical gust a input. 

Cr 

Figure 45.- Boeing rotor in autorotation at V/toR = 0.7 (Q = 386 rpm, 
V = 218 knots), magnitude of response of each degree of freedom to 
input at frequency co. 




Figure 45.- Continued. 
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(d) Lateral cyclic pitch 8^ 
Figure 45.- Continued. 
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Figure 45.- Continued. 





o)/X2 

(a) Vertical gust a input. 
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Figure 46.- Boeing rotor in powered operation at V/9R = 0.7 (fi = 386 rpm, 

V = 218 knots) , magnitude of response of each degree of freedom to input 
at frequency w. 
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(a) Frequency of the modes. , 
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Figure 47.- Boeing rotor velocity sweep at ft = 386 rpm; comparison with full 
scale experimental data and the Boeing theory. 
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(a) Wing vertical bending damping (cj'j) at 50 and 60 knots. 

(b) Wing vertical bending damping \q ) at 100 knots. 

(c) Wing vertical bending damping ( q j) at 140 knots. 

(d) Wing vertical bending damping (< 5 ^) at 192 knots. 

Figure 48.- Boeing rotor rpm sweep, comparison with full-scale experimental 

data and the Boeing theory. 
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